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Abstract 

Gamma- Ray Bursts (GRBs), short and intense pulses of low energy 7-rays , have fascinated 
astronomers and astrophysicists since their unexpected discovery in the late sixties. During the 
last decade, several space missions: BATSE (Burst and Transient Source Experiment) on Compton 
Gamma- Ray Observatory, BeppoSAX and now HETE II (High-Energy Transient Explorer), to- 
gether with ground optical, infrared and radio observatories have revolutionized our understanding 
of GRBs showing that they are cosmological, that they are accompanied by long lasting afterglows 
and that they are associated with core collapse Supernovae. At the same time a theoretical un- 
derstanding has emerged in the form of the fireball internal-external shocks model. According to 
this model GRBs are produced when the kinetic energy of an ultra-relativistic flow is dissipated 
in internal collisions. The afterglow arises when the flow is slowed down by shocks with the sur- 
rounding circum-burst matter. This model has numerous successful predictions like the prediction 
of the afterglow itself, the prediction of jet breaks in the afterglow light curve and of an optical 
flash that accompanies the GRBs themselves. In this review I focus on theoretical aspects and on 
physical processes believed to take place in GRBs. 
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I. INTRODUCTION 

Gamma-Ray Bursts (GRBs) are short and intense pulses of soft 7-rays . The bursts last 
from a fraction of a second to several hundred seconds. GRBs arrive from cosmological 
distances from random directions in the sky. The overall observed fluences range from 
10~ 4 ergs/cm 2 to 10 _7 ergs/cm 2 (the lower limit depends, of course, on the characteristic of 
the detectors and not on the bursts themselves). This corresponds to isotropic luminosity 
of 10 51 — 10 52 ergs/sec, making GRBs the most luminous objects in the sky. However, we 
know today that mos t GRBs ar e narrowly beamed and the corresponding energies are "only" 

31Clj | , making them comparable to Supernovae in the total energy 
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around 10 51 ergs 
release. 

The GRBs are followed by afterglow - lower energy, long lasting emission in the X- 
ray, optical and radio. The radio afterglow was observed in some cases several years after 
the bursts. The accurate afterglow positions enabled the identification of host galaxies in 
almost all cases when afterglow was detected and this in turn enabled the determination 
of the corresponding redshifts that range from 0.16 (or possibly even down to 0.0085) to 
4.5. Within the host galaxies there is evidence that (long duration) GRBs arise within star 
forming regions and there is evidence that they follow the star formation rate. 

While not all observed features are understood there is an overall agreement between the 
observations and the fireball model. According to the fireball model GRBs are produced 
when the kinetic energy of an ultra-relativistic flow is dissipated. The GRB itself is produced 
by internal dissipation within the flow while the afterglow is produced via external shocks 
with the circum-burst medium. I will focus in this review on this model. 

The numerous observations of the GRB and the observations of the afterglow constrain 
the fireball model that describes the emitting regions. The evidence on nature of the inner 
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engine that powers the GRB and produces the ultra- relativistic flow is however, indirect. 
The energetic requirements and the time scales suggest that GRB involve the formation of 
the black hole via a catastrophic stellar collapse event or possibly a neutron star merger. 
Additional indirect evidence arises from the requirement of the fireball model of long (several 
dozen seconds) activity of the inner engine. This hints towards an inner engine built on an 
accreting black hole. On the other hand, the evidence of association of GRBs with star 
forming regions indicates that GRBs progenitors are massive stars. Finally, the appearance 
of Supernova bumps in the afterglow light curve (most notably in GRB 030329) suggest 
association with Supernovae and stellar collapse. 

I review here the theory of GRB, focusing as mentioned earlier on the fireball internal- 
external shocks model. I begin in ^Hl with a brief discussion of the observations. I turn 
in £11111 to some generally accepted properties of GRB models - such as the essential ultra- 
relativistic nature of this phenomenon. Before turning to a specific discussion of the fireball 
model I review in §IVI several relativistic effects and in £0 the physical process, such as 
synchrotron emission or particle acceleration in relativistic shocks that are essential ingre- 
dients of this model. In WII I turn to a discussion of the prompt emission and the GRB. 
In Will I discuss modelling the afterglow emission. I consider other related phenomenon 
- such as TeV 7-rays emission, high energy neutrinos, Ultra High energy cosmic rays and 
gravitational radiation in Willi Finally, I turn in £11X1 to examine different 'inner engines' 
and various aspects related to their activity. I conclude with a discussion of open questions 
and observational prospects. 

While writing this review I realized how large is the scope of this field and how difficult 
it is to cover all aspects of this interesting phenomenon. Some important aspects had to be 
left out. I also did not attempt to give a complete historical coverage of the field. I am sure 
that inadvertently I have misse d ma ny i mportant references. I ref er the reader to several 
other recent review papers 

H00Q Q, H H that discuss these and 



other aspects of GRB theory and observations from different points of view. 

II. OBSERVATIONS 

I begin with a short review of the basic observed properties of GRBs. This review is brief 
as a complete review requires a whole paper by itself. I refer the reader to several review 
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papers for a detailed summary of the observations |98l . Ill9l 11771 l416j | . I divide this section to 
three parts. I begin with the prompt emission - the GRB itself. I continue with properties of 
the afterglow. I conclude with a discussion of the rates of GRBs, the location of the bursts 
within their host galaxies and the properties of the host galaxies. 

A. Prompt Emission 

I begin with a discussion of the GRB itself, namely the 7-rays and any lower-energy 
emission that occurs simultaneously with them. This includes X-ray emission that generally 
accompanies the 7-ray emission as a low energy tail. In some cases, called X-ray flashes 
(XRFs), the 7-rays signal is weak and all that we have is this X-ray signal. Prompt 
(operationally defined as the time period when the 7-ray detector detects a signal above 
background) longer-wavelength emission may also occur at the optical and radi o but it is 
harder to detect. However, so far optical flashes was observed In three cases 
simultaneously with the 7-ray emission. 

1 . Spectrum 

The spectrum is non thermal. The energy flux peaks at a few hundred keV and in many 
bursts there is a long high energy tail extending in cases up to GeV. The spectrum varies 
strongly from one burst to another. An excellent phenomenological fit for the spectrum 

n 

was introduced by Band et al. using two power laws joined smoothly at a break energy 
(a - P)E : 

[ (hvf exp(-^) for hu<(&- f3)E ; 

[[(a-^o] M, (M exp(^-a)> for hv > {a - P)E , 

I denote the spectral indices here as a and /3 to distinguish them from the afterglow param- 
eters (a and j3) discussed later. There is no particular theoretical model that predicts this 
spectral shape. Still, this function provides an excellent fit to most of the observed spectra. 
For most observed values of a and (3, vF v oc v 2 N{y) peaks at E p = (a + 2)E . For about 
10% of the bursts the upper slope is larger than -2 and there is no peak for vF„ w ithin the 



observed spectrum. Another group of bursts, NHE bursts, (no high energy) 296| does not 



have a hard component (which is reflected by a very negative value of $). The "typical" 



energy of the observed radiation is E p . E p denned in this way should not be confused with 
the commonly used hardness ratio which is the ratio of photons observed in two BATSE 1 
channels: Channel 3 (100-300keV) counts divided by Channel 2 (50-100keV) counts. The 



break frequen cy and t 



observed flux 
Band et al 



246, 



re peak flux frequencies are lower on average for bursts with lower 



247|. 



present a small catalogue of the 



spe ctra of 52 bright bursts which they 



3231 ] present a larger catalogue with 156 



analyze in terms of the Band function. Preece et al. 
bursts selected for either high flux or fluence. They consider several spectral shape including 
the Band function. 

Fig. [T]shows the d istribution of observed values of the break energy, (a—f3)E , in a sample 
of bright bursts 



323J. Most of the bursts are the range 100 keV < (a-(3)E < 400 keV, with 
a clear maximum in the distribution around (a — (3)Eq ~ 250keV. There are not many soft 
GRB s - that is, GRBs with peak energy in the tens of keV range. However, the discovery 



164] of XRFs - X-ray flashes with similar temporal structure to GRBs but lower typical 



energies - shows that the low peak energy cutoff is not real and it reflects the lower sensitivity 
of BATSE in this range fl8o( . 

Similarly, it is debatable whether there is a real paucity in hard GRBs and there is an 
upper cutoff to the GRB hardness or it just happens that the detection is optimal in this 
(a few hundred keV) band. BATSE triggers, for example, are based mostly on the count 
rate between 50keV and 300keV. BATSE is, therefore, less sensitive to harder bursts that 
emit most of their energy in the MeV range. Using BATSE's observation alone one cannot 
rule out the possibility that there is a population of harder GRBs that e mit equal powe r 
in total energy which are not observed because of this selection effect 
More generally, a harder burst with the same energy as a soft one emits fewer photons. 
Furthermore, the spectrum is generally flat in the high energy range and it decays quickly 
at low energies. Therefore it is intrinsically more difficult to detect a harder burst. A study 
of the SMM (Solar Maximum Mission) data |l58^ suggests that there is a deficiency (by 
at least a factor of 5) of GRBs with hardness above 3MeV, relative to GRBs peaking at 
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lit equal powe r 
I66.l22il.l3i3]. 



1 BATSE is the Burst and Transient Source Experiment on the CGRO (Compton Gamma-Ray Obser- 
vatory), see e.g. http://cossc.gsfc.nasa.gov/batse/ It operated for almost a decade detecting several 
thousand bursts, more than any satell ite before or after it. The BATSE data was published in several 
catalogues see Paciesas et al. |284i |285| for the most recent one 
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~0.5MeV, but this data is consistent with a population of hardness that extends up to 
2MeV. 

Overall the narrowness of the hardness distribution is very puzzling. First, as I stressed 
earlier it is not clear whether it is real and not a result of an observational artifact. If 
it is real then on one hand there is no clear explanation to what is the p hysic al process 



that controls the narrowness of the distribution (see however Guetta et al. |154j ). On the 
other hand cosmological redshift effects must broaden this distribution and it seem likely 
but n ot demon strated yet) that if the GRB distribution extends to z=10 as some suggest 



223] then such a narrow distribution requires an intrinsic correlation between the 



intrinsic hardness of the burst and its redshift, namely that the intrinsic hardness increases 
with the redshift. There is some evidence for such a correlation between E p and the observed 



peak flux 



246 



247(. More recently Amati et al. [5j reported on a correlation between E p 



and the isotropic equivalent energy seen in 12 BeppoSAX bursts that they have analyzed. 
They also report on a correlation between E p and the redshift as, the bursts with higher 
isotropic equivalent energy are typically more distant. These three different correlations are 
consistent with each other if the observed peak flux of bursts is determined by their intrinsic 
luminosity more than by the distance of the bursts. In such a case (because of the larger 
volume at larger distances) the observed more distant bursts are on average brighter than 
nearer ones (see also 311 Cj) . 

Even though the burst hardness distribution shows a single population a plot of the 
lardness vs temporal duration shows that short bursts (see Fig. 0} are typically harder 
79I 19?! . The correlation is significant. Another interesting sub-group of bursts is the NHE 
no h igh energy) bursts - bursts with no hard component that is no emission above 300keV 
296]. This group is characterized by a large negative value of (3, the high energy spectral 
slope. The NHE bursts have luminosities about an order of magnitude lower than regular 
bursts and they exhibit an effectively homogeneous intensity distribution with (V/V max ) = 
0.53 ± 0.029. As I discuss later in 31 1 A 21 most GRB light curves are composed of many 
individual pulses. It is interesting that in many bursts there are NHE pulses combined with 
regular pulses. 

EGRET (The Energetic Gamma Ray Experiment Telescope) the high energy 7-rays de- 
tector on Compton - GRO detected seven GRBs with photon energies ranging from 100 MeV 
to 18 GeV jSi|. In some cases this very high energy emission is delayed more than an hour 
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after the burst 1721 13921 ] . No hi gh-energy cutoff above a few MeV has been observed in any 
GRB spectrum. Recently, 13o| have combined the BATSE (30keV -2Mev) data with the 
EGRET data for 26 bursts. In one of these bursts, GRB 941017 (according to the common 
notation GRBs are numbered by the date), they have discovered a high energy tail that 
extended up to 200 MeV and looked like a different component. This high energy compo- 
nent appeared 10-20 sec after the beginning of the burst and displayed a roughly constant 
flux with a relatively hard spectral slope (F v oc u°) up to 200 sec. At late time (150 after 
the trigger) the very high energy (10-200 MeV) tail contained 50 times more energy than 
the "main" 7-rays energy (30keV-2MeV) band. The TeV detector, Milagrito, discovered (at 
a statistical significance of 1.5e-3 or so, namely at 3a) a TeV signal coincident with GRB 
970417 If true this would correspond to a TeV fluence that exceeds the low energy 7- 

rays fluence. However no further TeV signals were discovered from other 53 bursts observed 



252) . One 



by Milagrito [9] or from several bursts observed by the more sensitive Milagro 
should recall however, that due to the attenuation of the IR background TeV photons could 
not be detected from z > 0.1. Thus even if most GRBs emit TeV photons those photons 
won't be detected on Earth. 



Another puzzle is the low energy tail. Cohen et al. |55j analyze several strong burst s and 



find that their low energy slope is around 1/3 to -1/2. However, Preece et al. |32lL 322] 
suggest that about 1/5 of the bursts have a the low energy power spectrum, a, steeper than 
1/3 (the synchrotron slow cooling low energy slope). A larger fraction is steeper than -1/2 
(the fast cooling synchrotron low energy slope). However, this is not seen in any of the 
HETE spectrum whose low energy resolution is somewhat better. All HETE bursts have 
a low energy spectrum that is within the range 1/3 and -1/2 [16] . As both BATSE and 
HETE use Nal detectors that have a poor low energy resolution |55j , this problem might be 
resolved only when a better low energy spectrometer will be flown. 



2. Temporal Structure 

The duration of the bursts spans five orders, ranging from less than O.Olsec to more than 
lOOsec. Common measures for the duration are T90 (T50) which correspond to the time in 
which 90% (50%) of the counts of the GRB arrives. As I discuss below (see £111 A3|) the 
bursts are divided to long and short bursts according to their T90. Most GRBs are highly 
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variable, showing 100% variations in the flux on a time scale much shorter than the overall 
duration of the burst. Fig |21 depicts the light curve of a typical variable GRB (GRB 920627). 
The variability time scale, St, is determined by the width of the peaks. St is much shorter 
(in some cases by a more than a factor of 10 4 ) than T, the duration of t he burst. Variability 
on a time scale of milliseconds has been observed in some long bursts 25 ll |270(. However, 



only ~ 80% of the bursts show substantial substructure in their light curves. The rest are 
rather smooth, typically with a FRED (Fast Rise Exponenti al De cay) structure. 

Fenimore and Ramirez- Ruiz [9^] (see also Reichart et al. [34 1| ) discovered a correlation 
between the variability and the luminosity of the bursts. This correlation (as well as the 
lag-luminosity relation discussed later) allow us to estimate the luminosity of bursts that do 
not have a known redshift. 

The bursts seem to be composed of individual pulses, with a pulse being the "building 
blocks" of the overall light curve. Individual pulses display a hard to soft ev olution with 
the peak energy decreasing exponentially with the photon fluence 

HHQ- The P ulses 

have the following temporal and spectral features, (i) The light curve of an individual puls e 



is a FRED - fast rise exponential decay - with an average rise to decay ratio of 1 :3 [281]. 
(ii) T he lo w energy emission is delayed compared to the high energy emission 2 {28! | . Norris 
et al. |28o| ] have found that these spectral lags are anti-correlated with the luminosity of the 
bursts: Luminous bursts have long lags. This lag luminosity relation provides another way 
to estimate the luminosity of a burst from its (multi-spectra) light curve, (iii) The pulses' 
low energy light curves are wider compared to the high energy light curves. The width goes 
as ~ .E _0 ' 4 |9_2j. (iv) There is a Width-Symmetry-Intensity correlation. High intensity pulses 
are ( statistically) more symmetric (lower decay to rise ratio) and with shorter spectral lags 
(28! I . (v) There is a Hardness- Intensity correlation. The instantaneous spectral hardness 
of a pulse is correlated to the instantaneous intensity (the pulse become softer during the 
pulse decay) [loT ]. 

Both the pulse widths, St, and the pulse separation, At , have a rather similar log-normal 
distributi ons. However, the pulse separation, distribution, reveals has an exc ess of long 
0. These ,o ng — can be classified as quiescent periods |33o| . relatively 



2 Low/high energy implies the low vs. the high BATSE channels. The four BATSE channels at 20-50keV, 
50-100keV, 100-300keV and > 300keV. 
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long periods of several dozen seconds with no activity. When excluding these qu iescent 



270, £22 . The 



periods both distributions are log-normal with a comparable parameters 

average pulse interval, At = 1.3sec is larger by a factor 1.3 then the average pulse width 

St = lsec. One also finds that the pulse widths are correlated with the preceding interval 
_ I 1 

270]. Ramirez- Ruiz and Fenimore |328| found that the pulses' width does not vary along 
the bursts. 

One can also analyze the temporal behavior using the traditional Fourier transform 
method to analyze. The power density spectra (PDS) of light curves shows a power law 
slope of ~ —5/3 and a sharp break at 1Hz j^ . 

The results described so far are for long bursts, the variability of short (T < 2sec) bursts 
is more difficult to analyze. The duration of these bursts is closer to the limiti ng re solution 
of the detectors. Still most (~ 66%) short bursts are variable with St/T < 0.1 269]. These 
variable bursts are composed of multiple subpulses. 



3. Populations 



Long and Short Bursts The clearest classification of bursts is based on their duration. 
Kouveliotou et al. 198] have shown that GRB can be divided to two distinct gro ups: long 



170 



2QI 



burst with T 90 > 2sec and short bursts with T 90 < 2sec. Note that it was suggested 
that there is a third intermediate class with 2.5sec < Tqo < 7sec. However, it is not clear if 



156|. 



this division to three classes is statistically significant 

An interesting question is whether short bursts could arise from single pea ks of long 
bursts in which the rest of the long burst is hidden by noise. Nakar and Piran [269! have 
shown that in practically all long bursts the second highest peak is comparable in height 
to the first one. Thus, if the highest peak is above the noise so should be the second one. 
Short bursts are a diffe rent entity. This is supported by the observation that short bursts 



are typically harder 



79 



197| . The duration-hardness distribution (see Fig. HJ) shows clearly 
that there are not soft short bursts. 

The spatial distribution of the observed short bursts is clearly different from the dis- 
tribution of the observed long one. A measure of the spatial distribution is the average 
ratio (V/V max ) = {(C/C mm )" 3/ ' 2 ), where C is the count rate and C min is the minimal rate 
required for triggering. In a uniform Eucleadian sample this ratio equals 0.5 regardless of 
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the luminosity function. One of the first signs of a cosm olog ical origin of GRBs was the 



deviation of this value from 0.5 for t 
BATSE short bursts sample 



186 



248, 



re B ATSE sample 



3oJ 



256|. The the (V/V max ) of the 



is significantly higher than (V/V max ) of the long 



bursts sample. Note that more recently 380] suggested that the two values are s imila r and 



the distribution of long and short bursts is similar. However, Guetta and Piran [152j finds 
{V/V max )i ong = 0.282 and (V/V max )i ong = 0.390 (I discuss this point further in AH C 3|) . This 
implies that the population of observed short bursts is nearer on average than the popula- 
tion of the observed long ones. This is not necessarily a statement on the location of short 
vs. the location of long bursts. Instead it simply reflects the fact that it is more difficult to 
detect a short burst. For a short burst one has to trigger on a shorter (and hence noisier) 
window the detector (specifically BATSE that triggers on 64 ms for short bursts and on 1 
sec for long ones) is less sensitive to short bursts. I discuss later, in flll C 31 the question of 
rates of long vs. short bursts. 

So far afterglow was detected only from long bursts. It is not clear whether this is an 
observational artifact or a real feature. However, ther e was no X-ray afterglow observed for 
the only well localized short hard burst: GRB020531 Chandra observations show an 

intensity weaker by at least a factor of 100-300 than the intensity of the X-ray afterglow 
from long bursts at a similar time [46]. Afterglow was not observed in other wavelength as 



well 



121| 



As identification of hosts and redshifts depend on the detection of afterglow this implies 
that nothing is known about the distribution, progenitors, environment etc.. of short burst. 
These bursts are still waiting for their afterglow revolution. 

X-ray Flashes (XRFs) are X-ray burs ts with a similar temporal structure to GRBs 



but lower typical energies. Heise et al. 164j ] discovered these flushes by comparing GRBM 
(GRB Monitor) with sensitivity above 40 keV and WFC (Wide Field Camera) triggering 
on BeppoSAX 3 . In 39 cases the WFCs were triggered without GRBM triggering implying 
that these flashed do not have any hard component and most of their flux is in X-ray . 
The duration of 17 of these transients (out of the 39 transients), denoted X-ray flashes 
(XRFs), is comparable to the duration of the X-ray emission accompanying GRBs. The 
peak fluxes of the XRFs are similar to the X-ray fluxes observed during GRBs in the WFCs 



3 see http://www.asdc.asi.it/bepposax/ for information on BeppoSAX and its the different instruments. 
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10 8 ergs/sec/cm 2 ) but their pe ak en ergy is clearly below 40 keV. These finding confirmed 



the detection of Strohmayer et al. 



398] of 7 GRBs with E p < lOkeV and 5 additional GRBs 



with E p < 50keV in the GINGA data. 

Barraud et al. analyze 35 bursts detected on HETE II 4 . They find that XRFs lie 
on the extension of all the relevant GRB distributions. Namely there is a continuity from 
GRBs to XRFs. Detailed searches in the BATSE data revealed that some of these bursts 



have also been detect ed b y BATSE [189]. Using the complete search in 90% of the WFC 
data available, Heise jl63| find that the observed frequency of XRFs is approximately half 
of the GRB-frequency: In 6 years of BeppoSAX observations they have observed 32 XRFs 
above a threshold peak-luminosity of 5 x 10~ 9 erg/s/cm 2 in the 2-25 keV range compared 
with 54 GRBs (all GRBs above BATSE thresholds are observed if in the field of view). 

By now Soderberg et al. \38^ discovered optical afterglow from XRF 020903 and they 
suggest that the burst was at z = 0.25. They also suggest a hint of an underlying SN 
signal (see §IIC4J) peaking between 7-24 days after the initial XRF trigger. Afterglow was 
discovered from XRF 030723as well |l0o| . 



4- Polarization 

Recently, Coburn and Boggs |54j reported on a detection of a very high (80%±20%) linear 
polarization during the prompt 7-ray emission of GRB 021206. This burst was extremely 
powerful. The observed fluence of GRB 021206 was 1.6 • 10~ 4 ergs/cm 2 at the energy range 



of 25-100Kev 



175 



176]. This puts GRB 021206 as one of the most powerful bursts, and the 
most powerful one (a factor of 2-3 above GRB990123) after correcting for the fact that it 
was observed only in a narrow band (compared to the wide BATSE band of 20-2000keV). 
Coburn and Boggs |54j analyzed the data recorded by the Reuven Ramaty High Energy 
Solar Spectroscopic Imager (RHESSI). The polarization is measured in this detector by the 
angular dependence of the number detection of simultaneous pairs of events that are most 
likely caused by a scattering of the detected 7-rays within the detector. The data data 
analysis is based on 12 data points which are collected over 5sec. Each of these points is 
a sum of several independent observations taken at different times. Thus the data is some 



HETE II is a dedicated GRB satellite that aims at locating quickly bursts with high positional accuracy. 



See http:/ /space. mit.edu/HETE/ for a description of HETE II and its instruments 
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kind of convolution of the polarization over the whole duration of the burst. 

Coburn and Boggs [54j test two hypothesis. First they test the null hypothesis of no 
polarization. This hypothesis is rejected at a confidence level of 5.7a. Second they estimate 
the modulation factor assuming a constant polarization during the whole burst. The best 
fit to the data is achieved with II = (80 ± 20)%. However, Coburn and Boggs find that 
the probability that \ 2 is greater than the value obtained with this fit is 5%, namely the 
model of constant polarization is consistent with the analysis and observations only at the 
5% level. 

Rutledge and Fox [.SSoj ] reanalyzed this data and pointed out several inconsistencies within 



the methodology of Coburn and Boggs |54|. Their upper limit on the polarization (based 
on the same data) is ~ 4%. In their rebuttle Boggs and Coburn 39] point out that the 
strong upper limit (obtained by Rutledge and Fox [359] is inconsistent with the low S/N 
estimated by these authors. However, they do not p rovide a clear answer to the criticism 
of the methodology raised by Rutledge and Fox 359]. This leaves the situation, concerning 



the prompt polarization from this burst highly uncertain. 



5. Prompt Optical Flashes 

The Robotic telescope ROTSE (Robotic Optical Transient Search Experiment) detected a 
9th magnitude optical flash that was concurrent with the GRB emission from GRB 990123 
The six snapshots begun 40sec after the trigger and lasted until three minutes after 
the burst. The second snapshot that took place 60sec after the trigger recorded a 9th 
magnitude flash. While the six snapshots do not provide a "light curve" it is clear that the 
peak optical flux does not coincide with the peak 7-rays emission that takes place around 
the first ROTSE snapshot. This suggests t hat the optical flux is not the "low energy tail" 
of the 7-rays emission. Recently, Fox et al. 10l| reported on a detection of 15.45 magnitude 
optical signal from GRB 021004 193 sec after the trigger. This is just 93 seconds after 



2161 ] reported on 



the 100 sec long burst stopped being active. Shortly afterwards Li et al. 
a detection of 14.67 magnitude optical signal from GRB 021211 105 sec after the trigger. 
Finally, Price et al. |324j detected a 12th magnitude prompt flash, albeit this is more than 
1.5 hours after the trigger. Similar prompt signal was not observed f rom any other burst in 
spite of extensive searches that provided upper limits. Kehoe et al. searched 5 bright 
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bursts and found single-image upper limits ranging from 13th to 14th magnitude around 10 
sec after the initial burst detection and from 14 to 15.8 magnitudes one hour later. These 
upper limits are consistent with the two recent detections which are around 15th mag. The 
recent events of rapid detection suggest that we should expect many more such discoveries 
in the near future. 



6. The GRB- Afterglow Transition - Observations 

There is no direct correlation between the 7-ray fluxes and the X-ray (or optical) afterglow 
fluxes. The extrapolation of the X-ray afterglow fluxes backwards generally does not fit the 7- 
ray fluxes. Instead they fit the late prompt X-ray signal. These results are in nice agreement 
with the predictions of the Internal - External shocks scenario in which the two phenomena 
are produced by different effects and one should not expect a simple extrapolation to work. 

The expected GRB afterglow transition have been observed in several cases. The first 
observation took place (but was not reported until much latter) already in 1992 Bep- 
poSAX data shows a rather sharp transition in the hardness that takes place several dozen 
seconds after the beginning of the bursts. This transition is see n clearly in the different 



energy bands light curves of GRB990123 and in GRB980923 



129]. Connaughton 



have 



averaged the light curves of many GRBs and discovered long and soft tails: the early X- 

re a fterglow can be 



ray afterglow. Additional evidence for the transition from the GRB to t 
observed in the observations of the different spectrum within the GRB 

B. The Afterglow 

Until 1997 there were no known counterparts to GRBs in other wavelengths. On Feb 
28 1997 the Italian-Dutch satellite BeppoSAX detected X-ray afterglow from GRB 970228 



415 1 



IotI . The exact position given by BeppoSAX le d to the discovery of optical afterglow 
Radio afterglow was detected in GRB 970508 jloJ |. By now more than forty X-ray after- 



glows have been observed (see http://www.mpe.mpg.de/~jcg/grb.html for a complete up 



to date tables of well localized GRBs with or without afterglow. Another useful page is: 
http://grad40 .as.utexas.edu/grblog.php). About half of these have optical and radio after- 



glow (see Fig|HJ). The accurate positions given by the afterglow enabled the identification 
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of the host galaxies of many bursts. In twenty or so cases the redshift has been measured. 
The observed redshifts range from 0.16 for GRB 030329 (or 0.0085 for GRB 980425) to 
a record of 4.5 (GRB 000131). Even though the afterglow is a single entity I will follow 
the astronomical wavelength division and I will review here the observational properties of 
X-ray , optical and radio afterglows. 



1. The X-ray afterglow 

The X-ray afterglow is the first and strongest, but shortest signal. In fact it seems to 
begin already while the GRB is going on (see flll A 61 for a discussion of the GRB-afterglow 
transition). The light curve observed several hours after the burst can usually be extrapo- 
lated to the late parts of the prompt emission. 

The X-ray a fterg low fluxes from GRBs have a power law dependence on v and on the 



observed time t 



3161 ] : f u (t) oc v H a with a ~ 1.4 and ft ~ 0.9. The flux distribution, when 



normalized to a fixed hour after the burst has a rather narrow distribution. A cancellation of 
the k corrections and the temporal decay makes this flu x, wh ich is proport iona l to {l-\-zf~ a 



insensitive to the redshift. Using 21 BeppoSAX bursts |316| Piran et al. |310j find that the 
l-10keV flux, 11 hours after the burst is 5 x 10~ 13 ergs/cm -2 sec. The distribution is log- 
normal with <jf x pd 0.43 ± 0.1 (see fig. 03). De Pasquale et al. jz^J find a similar result 
for a larger sample. However, they find that the X-ray afterglow of GRBs with optical 
counterparts is on average 5 times brighter than the X-ray afterglow of dark GRBs (GRBs 
with no detected optical afterglow). The overall energy emitted in the X-ray afterglow is 
generally a few percent of the GRB energy. Berger et al. 26] find that the X-ray luminosity 
is indeed correlated with the opening angle and when taking the beaming correction into 
account they find that L x = fbLx,iso, is approximately constant, with a dispersion of only 
a factor of 2. 

X-ray lines were seen in 7 GRBs: GRB 970508 31^, GRB 970828 444], GRB 990705 
GRB 991216 0, GRB 001025a 428], GRB 000214 Q and GRB 011211 337]. The 



lines were detected using different instruments: BeppoSAX, ASCA (Advanced Satellite for 
Cosmology and Astrophysics) , Chandra and XMM-Newton. The lines were detected around 
10 hours after the burst. The typical luminosity in the lines is around 10 44 — 10 45 ergs/sec, 
corresponding to a total fluence of about 10 49 ergs. Most of the lines are interpreted as 
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emission lines of Fe Ka. However, there are also a radiative- recombination-continuum line 
edge an d Ka lines of lighter elements like Si, S, Ar and Ca (all seen in the afterglow of GRB 



011211 



337|). In one case (GRB 990705, Amati et al. j^) there is a transient absorption 



feature within the prompt X-ray emission, corresponding also to Fe Ka. The statistical 
significance of the detection of these lines is of some concern (2-5 a), and even thought the 
late instruments are much more sensitive tha n the early ones all detections remain at this 
low significance level. Rutledge and Sako |36o| and Sako et al. {36^ expressed concern about 
the statistical analysis of the data showing these lines and claim that none of the observed 
lines is statistically significant. The theoretical implications are far reaching. Not only the 
lines require, in most mo dels, a very large amount of Iron at rest (the lines are quite narrow), 



they most likely require 128j a huge energy supply (> 10 52 ergs), twenty time larger than 
the typical estimated 7-rays energy (~ 5 • 10 50 ergs). 



2. Optical and IR afterglow 

About 50% of well localized GRBs show optical and IR afterglow. The observed optical 
afterglow is typically around 19-20 mag one day after the burst (See fig El- The signal 
decays, initially, as a power law in time, t~ a with a typical value of a m 1.2 and large 
variations around this value. In all cases the observed optical spectrum is also a power law 
. Generally absorption lines are superimposed on this power law. The absorption lines 
correspond to absorption on the way from the source to earth. Typically the highest redshift 
lines are associated with the host galaxy, providing a measurement of the redshift of the 
GRB. In a few cases emission lines, presumably from excited gas along the line of site were 
also observed. 

Technical difficulties led a gap of several hours between the burst and the detection of the 
optical afterglow, which could be found only after an accurate position was available. The 
rapid localization provided by HETE II helped to close this gap and an almost complete 
light curve from 193 sec a fter the trigger (« 93 sec after the end of the burst) is available 



now for GRB021004 



1Q1 1 



Many afterglow light curves show an achromatic break to a st eeper dec line with a ~ 2. 
The classical example of such a break was seen in GRB 990510 159 . 139J ] and it is shown 



here in Fig. It is common to fit the break with the phenomenological formula: F„(t) 
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f*(t/t*)- ai {l - exp[-(t/^)( ai - a2 )](t/^)( ai " a2 )}. This break is commonly interpreted as a 



3 I Ij, H or the viewing 



jet break that allows us to estimate the opening angle of the jet 
angle within the standard jet model |34t| (see flll Dl below). 

The optical light cu rve o f the first detected afterglow (from GRB 970228) could be seen for 



more than half a year . In most cases the afterglow fades faster and cannot be followed 



for more than several weeks. At this stage the afterglow becomes significantly dimer than 
its host galaxy and the light curve reaches a plateau corresponding to the emission of the 
host. 

In a several cases: e.g. GRB 980326 



3J, GRB 970228 



3391 GRB 011121 



38, 



123] 



red bumps are seen at late times (several weeks to a month). These bumps are usually 
interpreted as evidence for a n underlyi ng SN. A most remarkable Supernova signature was 



seen recently in GRB 030329 



167 



395]. This supernova had the same signature as SN98bw 



that was associated with GRB 990425 (see §II04j) . 

Finally, I note that varying polarization at optical wav elengths 



afterglows at the level of a few to ten percent 
are in agreement with rough predictions 
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61 



62 



145 



ras been observed in GRB 



346 



4411 ] . These observations 



126l |365j) of the synchrotron emission model 



provided that there is a deviation from spherical symmetry (see §VFI below). 



3. Dark GRBs 



Only ~ 50% of well-localized GRBs show optical transients (OTs) successive to the 
prompt gamma-ray emission, whereas an X-ray counterpart is present in 90% of cases (see 
Fig. |BJ). Several possible explanations have been suggested for this situation. It is possible 
that late and shallow observations could not detect the OTs in some cases; several authors 
argue that dim and/or rapid de caying transients could bias the determination of th e fraction 



of truly obscure GRBs [24 



116|| . However, recent reanalysis of optical observations 



127 



210 



has shown that GRBs without OT detection (called dark GRBs, FOAs Failed Optical 
Afterglows, or GHOSTs, Gamma ray burst Hiding an Optical Source Transient) have had on 
average weaker optical counterparts, at least 2 magnitudes in the R band, than GRBs with 
OTs. Therefore, they appear to constitute a different class of objects, albeit there could be 
a fraction undetected for bad imaging. 

The nature of dark GRBs is not clear. So far three hypothesis have been put forward to 
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explain the behavior of dark GRBs. First, they are similar to the other bright GRBs, except 
for the fact that their line s of sight pass through large and dusty molecular clouds, that 
caus e hig h absorption 342]. Second, they are more distant than GRBs with OT, at z > 5 



110, 208], so that the Lyman break is redshifted into the optical band. Nevertheless, the 



distan ces o f a few dark GRBs have been determined and they do not imply high redshifts 



81 



,318]. A third possibility is that the optical afterglow of dark GRBs is intrinsically 



much fainter (2-3 ma g be low) than that of other GRBs. 



De Pasquale et al. 



74j find that GRBs with optical transients show a remarkably narrow 



distribution of flux ratios, which corresponds to an average optical-to-x spectral index 0.794± 
0.054. They find that, while 75% of dark GRBs have flux ratio upper limits still consistent 
with those of GRBs with optical transients, the remaining 25% are 4-10 times weaker 
in optical than in X-rays. This result suggests that the afterglows of most dark GRBs are 
intrinsically fainter in all wavelength relative to the afterglows of GRBs with observed optical 
transients. As for the remaining 25% here the spectrum (optical to X-ray ratio) must be 
different than the spectrum of other afterglows with a suppression of the optical band. 



4- Radio afterglow 

Radio afterglow was detected in ~ 50 % of the well localized bursts. Most observations are 
done at about 8 GHz since the detection falls off drastically at higher and lower frequencies. 
The observed peak fluxes are at the level of 2 mJy A turnover is seen around 0.2 mJy 
and the undetected bursts have upper limits of the order of 0.1 mJy. As the localization is 
based on the X-ray afterglow (and as practically all bursts have X-ray afterglow) almost all 
these bursts were detected in X-ray . ~ 80 % of the radio-afterglow bursts have also optical 
afterglow. The rest are optically dark. Similarly ~ 80% of the optically observed afterglow 
have also a radio component (see fig |BJ) • 

Several bursts (GRBs 980329, 990123, 91216, 000926, 001018, 010222, 011030, 011121) 
were detected at around one day. Recent radio observations begin well before that but do 
not get a detection until about 24 hrs after a b urst. The earliest radio detection took place in 



GRB 011030 at about 0.8 days after the burst |40jj. In several cases (GRBs 990123, 990506, 



991216, 980329 and 020405) the afterglow was detected early enough to indicate emission 
from the reverse shock and a transition from the reverse shock to the forward shock. 
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The radio light curv e of GRB 97050 8 (se e fig EJ) depicts early strong fluctuations (of 



order unity) in the flux |l04j . Goodman [1311 ] suggested that these fluctuations arise due to 



scintillations and the decrease (with time) in the amplitude of t he fluctuations arises from 
a transition from strong to weak scintillations. Frail et al. jl04 ] used this to infer the size 
of the emitting region of GRB 970508 at ~ 4 weeks after the burst as ~ 10 17 cm. This 
observations provided the first direct proof of relativistic expansion in GRBs. 

The self-absorbed frequencies fall in the centimeter to meter wave radio regime and hence 
the lower radio emission is within the self-a bsorption part of the spectrum (see WG 31 later). 
In this case the spectrum rises as v 2 \\$>\ . The spectral shape that arises from a the fact 
that the system is optically thick enables us (using similar arguments to those of a simple 
black body emission) to determine the size of the emitting region. In GRB 990508 this has 
lead to ~ 10 17 cm. A comparable estimate to the one derived from scintillations. 

The long-lived nature of the radio afterglow allows for unambiguous calorimetry of the 
blast wave to be made when its expansion has become sub-relativistic and quasi-spherical. 
The light curves evolves on a longer time scale in the radio. Some GRB afterglows have been 
detected years after the burst even after the relativistic-Newtonian transition (see §VII Djl . 
At this stage the expansion is essentially sphe rical and this enables a direct " calorimetric" 
estimate of the total energy within the ejecta |436 1 . 

C. Hosts and Distribution 

1 . Hosts 

By now (early 2004) host galaxies have been observed fo r all but 1 or 2 bursts with optical, 
radio or X-ray afterglow localization with arcsec precision \yi\ . The no- host problem which 
made a lot of noise in the nineties has disappeare d. G RBs are located within host galaxies 



(see Djorgovski et al. 



1771 ] for detailed reviews). While many 



82J, |8J and Hurley et al. 
researchers believe that the GRB host population seem to be representative of the normal 
star-forming field galaxy population at a comparable redshifts, others argue that GRB host 
galaxies are significantly bluer than average and their star formation rate is much higher 
than average. 

The host galaxies are faint with median apparent magnitude R ~ 25. Some faint hosts 
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are at R ~ 29. Down to R ~ 25 the observed distribution is consistent with deep field 
galaxy counts. Jimenez et al. jl82| find that the likelihood of finding a GRB in a galaxy is 
proportional to the galaxy's luminosity. 

The magnitude and redshift distribution of GR B host galaxie s are typical for normal, 



faint field galaxies, as are their morphologies [36, 



84 



168, 



1771 . |282j . While some researchers 



argue that the broad band optical colors of GRB hosts are not dist ingu ishable f rom those of 
normal field galaxies at comparable magnitudes and redshifts |36L 1391 1. others 110 ] asserts 
that t he ho st galaxies are unusually blue and that they are strongly star forming. Le Floc'h 
et al. I213I argues that R-K colors of GRB hosts a re u nusually blue and the hosts may be 



of low metallicity and luminosity. This suggests [212J that the hosts of GRBs might be 
different from the cites of the majority of star forming galaxies tha t are luminous, reddened 
and dust-e nshr ouded infrared starbursts (Elbaz and Cesarsky 



Le Floc'h 



and references therein). 
212j | also suggests that this difference might rise due to an observational bias 



and that GRBs that arise in dust-enshrouded infrared starbursts are dark GRBs whose 
afterglow is not detectable due to obscuration. Whether this is tru or not is very relevant 
to the interesting question to which extend GRBs follow the SFR and to which extend they 
can be use d to determine the SFR at high redshifts 
Totani 40s| . Wijers et al. \i3^ and Paczynski 



287] suggested that GRBs follow the star 
formation rate. As early as 1998 Fruchter et al. jl lo| noted that all four early GRBs with 
spectroscopic identification or deep multicolor broadband imaging of the host (GRB 970228 
GRB 970508, GRB 971214, and GRB 980703) lie in rapidly star-forming galaxies. Within 
the host galaxies the distribution of GRB-host offset follows the light distribution of the 
hosts j3| . The light is roughly proportional to the density of star formation. Spectroscopic 
measurements suggest that GRBs are within Galaxies with a higher SFR . However, this is 
typical for normal field galaxy population at comparable redshifts jl73| . There are some 
intriguing hints, in particular the flux ratios of [Ne III] 3859 to [Oil] 3727 are on average 
a factor of 4 to 5 higher in GRB hosts than in star forming galaxies at low redshifts [84]. 
This may represent an indirect evidence linking GRBs with massive star formation. The 
link between GRBs and massive stars has been stren gthened with the centimeter and sub- 
— d_ of GRB host gaWs Q Q undergoing p r od igioUS StaI formation 
(SFR~ 10 3 M yr _1 ), which remains obscured at optical wavelengths. 

Evidence for a different characteristics of GRB host galaxies arise from the work of Fynbo 
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et al. 



114, 



1151 ] who find that GRB host galaxies "always" show Lyman alpha emission in 



cases where a suitable search has been conducted. This back up the claim for active star 
formation and at most moderate metallicity in GRB hosts. It clearly distinguishes GRB 
hosts from the Lyman break galaxy population, in which only about 1/4 of galaxies show 
strong Lyman alpha. 



2. The Spatial Distribution 

256| was among 



BATSE's discovery that the bursts are distributed uniformly on the sky 
the first indication of the cosmological nature of GRBs. The uniform distribution indicated 
that GRBs are not associated with the Galaxy or with "local" structure in the near Universe. 

Recently there have been several claims that sub-groups of the whol e GRB population 
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230], for example, find 



shows a deviation from a uniform distribution. Meszaros et al. 
that the angular distribution of the intermediate sub-group of burst s (m ore specifically of 
the weak intermediate sub-group) is not random. Magliocchetti et al. |245| reported that the 
two-point angular correlation function of 407 short BATSE GRBs reveal a ~ 2a deviation 
from isotropy on angular scales 2° — 4°. This results is consistent with the possibility that 
observed short GRBs are nearer and the angular correlation is induced by the large scale 
structure correlations on this scale. These claims are important as they could arise only if 
;hese bursts are relatively nearby. Alternatively this indicates repetition of these sources 



245]. Any such deviation would imply that these sub-groups are associated with different 

objects than the main GRB population at least that these subgroup are associated with a 

specific feature, such different viewing angle, 
rn 

Cline et al. [53] studied the shortest GRB population, burst with a typical durations 
several dozen ms. They find that there is a significant angular asymmetry and the (V/V max ) 
distribution provides evidence for a homogeneous sources distribution. They suggest that 
these features are best interpreted as sources of a galactic origin. However, one has to realize 
that there are strong selection effects that are involved in the detection of this particular 
subgroup. 
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3. GRB rates and the isotropic luminosity function 



There have been many attempts to determine the GRB luminosity function and rate 
from the BATSE peak flux distribution. This was done by numerous using different levels 
of statistical sophistication and different physical assumptions on the evolution of the rate 
of GRBs with time and on the shape of the luminosity function. 

Roughly speaking the situation is the following. There are now more than 30 redshift 
measured. The median redshift is z ~ 1 and the redshift range is from 0.16 (or even 
0.0085 if the association of GRB 980425 with SN 98bw should be also considered) to 4.5 
(for GRB 000131). Direct estimates from the sample of GRBs with determined redshifts 
are contaminated by observational biases and are insufficient to determine the rate and 
luminosity function. An alternative approach is to estimates these quantities from the 
BATSE peak flux distribution. However, the observed sample with a known redshifts clearly 
shows that the luminosity function is wide. With a wide luminosity function, the rate of GRB 
is only weakly constraint by the peak flux distribution. The analysis is further complicated 
by the fact that the observed peak luminosity, at a given detector with a given observation 
energy band depends also on the intrinsic spectrum. Hence different assumptions on the 
spectrum yield different results. This situation suggest that there is no point in employing 



sophisticated statistical tools (see however, 224L l308| | for a discussion of these methods) and 



a simple analysis is sufficient to obtain an idea on the relevant parameters. 

I will not att empt to review the various approaches he re. A partial list of calculations 



includes 



56 



90 



169 



224, 



225 



304, 



308, 



379 



380, 



381 



382J. Instead I will just quote results 



of some estimates of the rates and luminosities of GRBs. The simplest ap proach is to fi ; 
(V/V max ), which is the first moment of the peak flux distribution. Schmidt 379, 380 . 138 1 1 



inds using {V/V max ) of the long burst distribution and assuming that the bursts follo w th e 

38o|. 
This 



320j | SFR2, that the pr esent local rate of long observed GRBs is ~ 0.15Gpc~ 3 yr 
Note that this rate from |38o| is smaller by a factor of ten than the earlier rate of 3j 
estimate corresponds to a typical (isotropic) peak luminosity of ~ 10 51 ergs/sec. These are 
the observed rate and the isotr opic peak luminosity. 

Recently Guetta et al. 153 1 have repeated these calculations . They use both the [352! 
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SFR formation rate: 

f 10 - 752 z < 1 
Rgrb{z)= P A , (2) 

^g0.752 pcak ^ > 1. 

and SFR2 from |32o| . Their best fit luminosity function (per logarithmic luminosity interval, 
cHogL) is: 

f (L/L*) a L*/30 <L<L* 
$o(L) = cJ R (3) 

[ (L/L*f L*erg/sec <L< 30L* 

and otherwise with a typical luminosity, L* = 1.1 x 10 51 ergs/sec, a = —0.6 and (3 = —2, 
and c Q is a normalization constant so that the integral over the luminosity function equals 
unity. The corresponding local GRB rate is po — 0.44Gpc _1 yr _1 . There is an uncertainty 
of a factor of ~ 2 in the typical energy, L*, and in the local rate. I will use these numbers 
as the "canonical" values in the rest of this review. 

The observed (BATSE) rate of short GRBs is smaller by a factor of three than the rate 
of long ones. However, this is not the ratio of the real rates as :(i) The BATSE detector 
is less sensitive to short bursts than to long ones; (ii) The true rate depends on the spatial 
distribution of the short bursts. So far no redshift was detected for any short bursts and 
hence this distribution is uncertain. For short bursts we can resort only to estimates based 
on the peak flux distribution. There are indications that (V ~/V max ) of short burst is larger 
(and close to the Eucleadian value of 0.5) than the {V/V max ) value of long ones (which is 
around - 32). This implies that the observed short bursts are nearer to us that the long ones 
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403| possible with all observed short bursts are at z < 0.5. However, Schmidt [380] 
finds for short bursts (V/V max ) = 0.354, which is rather close to the value of long bursts. 
Assuming that short GRBs also follow the SFR he obtains a local rate of 0.075Gpc~ 3 yr _1 - a 
factor of two below the rate of long GRBs! The (isotropic) peak lumi nosit ies are comparable. 
This results differs from a recent calculation of Guetta and Piran |l52| who find for short 
bursts (V/Vmax) — 0.390 and determine from this a local rate of 1.7Gpc" 3 yr _1 which is 
about four times the rate of long bursts. This reflects the fact that the observed short 
GRBs are significantly nearer than the observed long ones. 

These rates and luminosities are assuming that the bursts are isotropic. Beaming reduces 
the actual peak luminosity increases the implied rate by a factor fj~ = 2/9 2 . By now there 
is ev i denc e that GRBs are beamed and moreover the total energy in narrowly distributed 



105 



29 1| . There is also a good evidence that the corrected peak luminosity is much more 
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419]. The corrected peak 



narrowly distributed than the isotropic peak lumi nosity 1153 

luminosity is L pea k(6 2 /2) ~ const. Frail et al. loj | suggest that the true rate is larg er by 
a factor of 500 than the observed isotropic estimated rate. However, Guetta et al. [153] 
repeated this calculation performing a careful average over the luminosity function and find 
that that true rate is only a factor of ~ 75 ± 25 times the isotropically estimate one. Over 
all the true rate is: 33 ± ll/ig 5 Gpc _3 yr _1 . 

With increasing number of GRBs with redshifts it may be possible soon to determine 
the GRB redshift distribution directly from this data. However, it is not clear what are 
the observational biases that influence this data set and one needs a homogenous data set 
in order to per form this calc ulation. Alternatively one can try to determine luminosity 



estimators 95, 28 



3771 1378| from the subset with known redshifts and to obtain, using 



them a redshift distribution for the whole GRB sample. Lloyd-Ronning et al. |223| find 
using the Fenimore and Ramirez-Ruiz sample that this method implies that (i) The 
rate of GRBs increases backwards with time even for z > 10, (ii) The Luminosity of GRBs 
increases with redshift as (1 + z) 1A±0 ' 5 ; (iii) Hardness and luminosity are strongly correlated. 
It is not clear how these features, which clearly depend on the inner engine could depend 
strongly on the redshift. Note that in view of the luminosity-angle relation (see flll Dl below) 
the luminosity depends mostly on the opening angle. An increase of the luminosity with 
redshift would imply that GRBs were more narrowly collimated at earlier times. 



4- Association with Supernovae 

The association of GRBs with star forming regions and the indications that GRBs follow 
;he s tar formation rate suggest that GRBs are related to stellar death, namely to Supernovae 



2871 ] . Additionally there is some direct evidence of association of GRBs with Supernovae. 

GRB 980425 and SN98bw: The first indication of an association between GRB s and 
SNes was found when SN 98bw was discovered within the error box of GRB 980425 120]. 



This was an usual type Ic SN which was much brighter than most SNs. Typical ejection 
velocities in the SN were larger than usual (~ 2 ■ 10 4 km/sec) corresponding to a k ineti c 
energy of 2 — 5 52 ergs, more than ten times than previously known energy of SNes, jl To| . 
Addition ally radio observations suggested a component expanding sub relativistically with 
v ~ 0.3c |200j |. Thus, 1998bw was an unusual type Ic supernovae, significantly more powerful 
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than comparable SNes. This may imply that SNs are associated with more powerful SNes. 
Indeed all other observations of SN signature in GRB afterglow light curves use a SN 98bw 
templates. The accompanying GRB, 980425 was also unusual. GRB 980425 had a smooth 
FRED light curve and no high energy component in its spectrum. Other bursts like this 
exist but they are rare. The redshift of SN98bw was 0.0085 implying an isotropic equivalent 
energy of ~ 10 48 ergs. Weaker by several orders of magnitude than a typical GRB. 

The BeppoSAX Wide Field Cameras had localized GRB980425 with a 8 arcmin radius 
accuracy. In this circle, the BeppoSAX NFI (Narrow Field Instrument) had detected two 
sources, SI and S2. The NFI could associate with each of these 2 sources an error circle 
of 1.5 arcmin radius. The radio and optical position of SN1998bw were consistent only 



le NFI error circle of SI, and was out of the NFI error circle of S2. Therefore, Pian 
identified SI with X-ray emission from SN1998bw, although this was of course 



with t 
et al. 

no proof of association between SN and GRB. It was difficult, based only on the BeppoSAX 
NFI data, to characterize the behavior and variability of S2 and it could not b e exc luded 
that S2 was the afterglow of GRB980425. The XMM observations of March 2002 Q seem 
to have brought us closer to the solution. XMM detects well SI, and its flux is lower than 
in 1998: the SN emission has evidently decreased. Concerning the crucial issue, S2: XMM, 
having a better angular resolution than BeppoSAX NFIs, seems to resolve S2 in a number 
of sources. In other words, S2 seems to be not a single source, but a group of small faint 
sources. Their random variability (typical fluctuations of X-ray sources close to the level 
of the background) may have caused the flickering detected for S2. This demolishes the 
case for the afterglow nature of S2, and strengthens in turn the case for association between 
GRB980425 and SN1998bw. 
Red Bumps: The 



light curves 
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38, 



123 



ate r ed bumps (see flll B 2|) have been discovered in several GRB 
l33£j| . These bumps involve both a brightening (or a flattening) of 
the afterglow as well as a transition to a much redder spectrum. These bumps have been 

n 

generally interpreted as due to an underlining SN [3J3|. In all cases the bumps have been fit 
with a template of SN 1998bw, which was associated with GRB 980425. Esin and Blandford 
88] proposed that these bum ps a re produced by light echoes on surrounding dust (but see 
34o| ). Waxman and Draine 435| purposed another alternative explanation based on dust 
sublimation. 

For most GRBs there is only an upper limit to the magnitude of the bump in the light 
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curve. A comparison of these upper limits (see Fig. ITUj) with the maximal magnitudes of type 
Ibc SNe shows that the faintest GRB-SN non-detection (GRB 010921) only probes the top 
~40th-percentile of local Type Ib/Ic SNe. It is clear that the current GRB-SNe population 
may have only revealed the tip of the iceberg; plausibly, then, SNe could accompany all 
long-duration GRBs. 

GRB 030329 and CN 2003dh: The confirmation of SN 98bw like bump and the 
confirmation of the GRB-SN association was dramatically seen recently 



167, m& in the 



very bright GRB 030329 that is associated with SN 2003dh j5l|. The bump begun to be 
noticed six days after the bursts and the SN 1999bw like spectrum dominated the optical 
light curve at later times (see Fig. ^2 The spectral shapes of 2003dh and 1998bw were 
quite similar, although there are also differences. For example 111 Dl estimated a somewhat 
larger expansion velocity for 2003dh. Additionally the X-ray signal was much brighter (but 
this could be purely afterglow). 

For most researchers in the field this discovery provided the final conclusive link between 
SNe and GRBs (at least with long GRBs). As the SN signature coincides with the GRB 
this observations also provides evidence against a Supranova interpretation, in which the 
GRB arises from a collapse of a Neutron star that takes place sometime after the Supernova 
in which the Neutron star was born - see IIX El . (unless there is a variety of Supranova 
types, some with long delay and others with short delay between the first and the second 
collapses) the spectral shapes of 2003dh and 1998bw were quite similar, although there are 
also differences. For example there is a slightly larger expansion velocity for 2003dh. It is 
interesting that while not as week as GRB 990425, the accompanying GRB 99030329 was 
significantly weaker than average. The implied opening angle reveals that the prompt 7-ray 
energy output, _E 7 , and the X-ray luminosity at 10 hr, Lx, are a factor of ~ 20 and ~ 30, 
respectively, below the average value around which most GRBs are narrowly clustered (see 
HID] below). 

It i s int eresting to compare SN 1999bw and SN 2003dh. Basically, at all epochs Matheson 



et al. |249l | find that the best fit to spectra of 2003dh is given by 1998bw at about the same 
age . The light curve is harder, as the afterglow contribution is significant, but using spectral 
inform ation they find that 2003dh had basically the same light curve as 1998bw. Mazzali 
et al. j25o| model the spectra and find again that it was very similar to 1998bw. They 
find some differences, but some of that might be due to a somewhat different approach to 
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spectral decomposition, which gives somewhat fainter supernova. 

X-ray lines: The appearance of iron X-ray lines (see flll B 1|) has been interpreted as 
additional evidence for SN. One has to be careful with this interpretation as the iron X- 
ray lines are seen as if emitted by matter at very low velocities and at rather large distances. 
This is difficult to achieve if the supernova is simultaneous with the GRB , as t he SN bumps 



imply. This X-ray lines might be consistent with the Supranova model 



4261 ] in which the 



SN takes place month before the G RB. However, in this case ther e won't be a SN bump in 



the 



203 



ight curve! Meszaros and Rees 



239J, Rees and Meszaros |33q and Kumar and Narayan 



suggest alternative interpretations which do not require a Supranova. 



D. Energetics 



Before redshift measurements were available the GRB energy was estimated from the 
BATSE catalogue by fitting an (isotropi c) lu minosity function to the fl ux d istribution (see 
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n 



e.g Cohe n and Piran Guetta et al. |l53| |. Loredo and Wasserman [225j | . Schmidt 
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381] and many others). This lead to a statistical estimate of the luminosity function of 
a distribution of bursts. 

These estimates were revolutionized with the direct determination of the redshift for 
individual bursts. Now the energy could be estimated directly for specific bursts. Given 
an observed 7-ray fluence and the redshift to a burst one can easily estimate the energy 
emitted in 7-rays, E 7jiso assuming that the emission is isotropic (see Bloom et al. for 
a detailed study including k corrections). The inferred energy, E ltiso was the isotropic 
energy, namely the energy assuming that the GRB emission is isotropic in all directions. 
The energy of the first burst with a determined redshift, GRB 970508, was around 10 51 ergs. 
However, as afterglow observations proceeded, alarmingly large values (e.g. 3.4 x 10 54 ergs for 
GRB990123) were measur ed fo r a so . The variance was around three orders of magnitude. 



However, it turned out 



344, 



374j that GRBs are beamed and E liiso would not then be a 
good estimate for the total energy emitted in 7-rays. Instead: E^ = (Q 2 /2)E 1 ^ iso . The angle, 
6, is the effective angle of 7- ray e mission. It can be estimated from the time of the break 
in the afterglow light curve |374j : 



6 = U.lQ{n/E k , so ^y/hl% ys = Omin/E^e^XJays, 



1/6,1/2 



(4) 
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where todays is the break time in days. Ek^ so ,52 is "isotropic equivalent" kinetic energy, 
discussed below, in units of 10 52 ergs, while Ek,e,b2 is the real kinetic energy in the jet i.e: 
Ek,e,h2 — {@ 2 /2)Ek,i SO ,52- One has to be careful which of the two energies one discusses. In 
the following I will usually consider, unless specifically mentioned differently, Ek t i SO ,52, which 
is also related to the energy per unit solid angle as: -Efc,iso,52/4vr. The jet break is observed 



both in the optical and in the radio frequencies. IN 
in the radio differs from that at optical and X-ray 
an additional confi rmation for this interpretation. 



ote that the the observational signature 



159 



374j (see Fig. 125]) and this provides 



Frail et al. [1051 ] esti mated E 1 for 18 bursts, finding typical values around 10 51 ergs (see also 



Panaitescu and Kumar 29l|). Bloom et al. (34 1 find = 1.33 x 10 51 erg and a burst-to- 
burst variance about this value ~ 0.35 dex, a factor of 2.2. This is three orders of magnitude 
smaller than the variance in the isotropic equivalent E T A compilation of the beamed 
energies from 2], is shown in Figs ^] and ED It demonstrates nicely this phenomenon. 
The constancy of E 1 is remarkable, as it involves a product of a factor inferred from the 
GRB observation (the 7-rays flux) with a factor inferred from the afterglow observations 
(the jet opening angle). However, E 7 might not be a good estimate for E to t, the total energy 
emitted by the central engine. First, an unknown conversion efficiency of energy to 7-rays has 



to be considered: E- 



tot 



E^ 



l /2)E li i S0 . Second, the large Lorentz factor during 



the 7-ray emission phas e, ma kes the observed E 1 rather sensitive to angular inhomogeneities 



of the relativistic ejecta 



2061 ] . The recent early observations of the after glow of G RB 021004 
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indicate that indeed a significant angular variability of this kind exists 

The kinetic energy of the flow during the adiabatic afterglow phase, Ef. is yet another 
energy measure that arises. This energy (per unit solid angle) can be estimated from the 
afterglow light curve and spectra. Specifically it is rather closely related to the observed 



afterglow X-ray flux lOoi 201 . 310 ]. As this energy is measured when the Lorentz factor is 
smal ler it is less sensitive than Ej to angular variability. The const ancy of the X-ray flux 
suggest that this energy is also constant. Estimates of E^o \l9 1| show that E 1 « 3Ek t e, 
namely the observed "beamed" GRB ener gy is larger than the estimated "beamed" kinetic 
energy of the afterglow. Frail et al. \l0^ . however, find that E 1 « E^g, namely that the 
two energies are comparable. 

An alternative interpretation to the observed breaks is that we are viewi ng a " universal 



angle dependent, namely, "structured" jet - from different viewing angles 
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347, |44( 



30 



The observed break corresponds in this model to the observing angle 6 and not to the opening 
angle of the jet. This interpretation means that the GRB beams are wide and hence the 
rate of GRBs is smaller than the rate implied by the usual beam ing factor. On the other 



hand it implies that GRBs are more energetic. Guetta et al. jl 531 ] estimate that this factor 
(the ratio of the fixed energy of a "structured" jet relative to the energy of a uniform jet to 
be ~ 7. However they find that the observing angle distribution is somewhat inconsistent 
with the s imple geomet ric o ne that should arise in universal structured jets (see also Nakar 



et al. 



2661 ] . Perna et al. |299[). The energy-angle relation discussed earlier require (see WII II 
below) an angle dependent jet with E{6) oc 9~ 2 . 

Regardless of the nature of the jet (universal structured jet or uniform with a opening 
angle that differs from one burst to another) at late time it becomes non relativistic and 
spherical. With no relativistic beaming every observer detects emission from the whole shell. 
Radio observations at this stage enable us to o btain a direct calorimetric estimate of the 
total kinetic energy of the ejecta at late times 107^ Estimates performed in several cases 
yield a comparable value for the total energy. 

If GRBs are beamed we should expect orphan afterglows (see §VIIK)1 : events in which 
we will miss the GRB but we will observe the late afterglow that is not so beamed. A 
comparison of the rate of orphan afterglows to GRBs will give us a direct estimate of the 
beaming of GRBs (and hence of their energy). Unfortunately there are not even good upper 



limits on the rate of orphan afterglows. Veerswijk 421] consider the observations within 
the Faint Sky Variability Survey (FSVS) carried out on the Wide Field Camerea on teh 
2.5-m Isacc Newton Telescope on La Palma. This survey mapped 23 suare degree down to 
a limiting magnitude of about V=24. They have found one object which faded and was 
not detected after a year. Howeve r, its colors suggest that it was a supernova and not a 



GRB. Similarly, Vanden Berk et al. 



420l | find a single candidate within the Sloan Digital Sky 



Survey. Here the colors were compatible with an afterglow. However, later it was revealed 
that this was a variable AGN and not an orphan afterglow. As I discuss later this limits are 
still far to constrain the current beaming estimates (see §VIIK)1 . 



One ex cept ion is for late radio emission for which there are some limits 



son et al. 
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2m- Levin- 



215J show that the number of orphan radio afterglows associated with GRBs that 



should be detected by a flux-limited radio survey is smaller for a smaller jet opening angle 9. 
This might seen at first sight contrary to expectation as narrower beams imply more GRBs. 



31 



But, on the other hand, with narrower beams each GRB has a lower energy and hence its 
radio afterglow is more difficult to detect. Overall the second factor wins. Using the results 
of FIRST and NVSS surveys they find nine afterglow candidates. If all candidates are asso- 
ciated with GRBs then there is a lower limit on the beaming factor of / b ~ = (# 2 /2) > 13. If 
none are associated with GRBs they find / 6 > 90. This giv e im mediately a corresponding 
upper limit on the average energies of GRBs. Guetta et al. jl53| revise this values in view 
of a recent estimates of the correction to the rate of GRBs to: = 40. 

When considering the energy of GRBs one has to remember the possibility, as some 
models suggest, that an additional energy is emi tted which is not involved in the GRB itself 
or in the afterglow, van Putten and Levinson 4is| . for example, suggest that a powerful 
Newtonian wind collimates the less powerful relativistic one. The "standard jet" model 
also suggests a large amount of energy emitted sideways with a lower energy per solid 
angle and a lower Lorentz factors. It is interesting to note that the calorimetric estimates 
mentioned earlier limit the total amount of energy ejected regardless of the nature of the flow. 
More generally, typically during the afterglow matter moving with a lower Lorentz factor 
emits lower frequencies. Hence by comparing the relative beaming of afterglow emission 
in different wavelength one can estimate the relative beami ng fa ctors, / 6 _1 (i?), at different 



wavelength and hence at different energies. Nakar and Piran |272j use various X-ray searches 
for orphan X-ray afterglow to limit the (hard) X-ray energy to be at most comparable to 
the 7-rays energy. This implies that the total energy of matter moving at a Lorentz factor 
of ~ 40 is at most comparable to the energy of matter moving with a Lorentz factor of a 
few hundred and producing the GRB itself. At present limits on optical orphan afterglow 
are insufficient to set significant limits on matter moving at slower rate, while as mentioned 
earlier radio observations already limit the overall energy output. 

These observations won't limit, of course, the energy emitted in gravitational radiation, 
neutrinos, Cosmic Rays or very high energy photons that may be emitted simultaneously 
by the source and influence the source'e energy budget without influencing the afterglow. 

III. THE GLOBAL PICTURE - GENERALLY ACCEPTED INGREDIENTS 

There are several generally accepted ingredients in practically all current GRB models. 
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Relativistic Motion: Practically all current GRB models involve a relativistic motion 
with a Lorentz factor, T > 100. This is essential to overcome the compactness problem 
(see £11 V Al below) . At first this understanding was based only on theoretical arguments. 
However, now there are direct observational proofs of this concept: It is now gene rally 



accepted that both the radio scintillation 



131] and the lower frequency self-absorption [187] 



provide independent estimates of the size of the afterglow, ~ 10 17 cm, two weeks after the 
burst. These o bserv ations imply that the afterglow has indeed expanded relativistically. 
Sari and Piran 372^ suggested that the optical flash accompanying GRB 990123 provided 
a dir ect evidence for ultra-relativistic motion with T ~ 100. Soderberg and Ramirez- Ruiz 
39o| find a higher value: 1000 ± 100. However, these interpretations are model dependent. 

The relativistic motion implies that we observe blue shifted photons which are signifi- 
cantly softer in the moving rest frame. It also implies that when the objects have a size R the 
observed emission arrives on a typical time scale of R/cT 2 (see ^IVB|) . Relativistic beaming 
also implies that we observe only a small fraction (1/r) of the source. As I discussed earlier 
(see §IIDI and also IIV (^|) this has important implications on our ability to estimate the total 
energy of GRBs. 

While all models are based on ultra-relativistic motion, none explains convincingly (this 
is clearly a subjective statement) how this relativistic motion is attained. There is no 
agreement even on the nature of the relativistic flow. While in some models the energy 
is carried out in the form of kinetic energy of baryonic outflow in others it is a Poynting 
dominated flow or both. 

Dissipation In most models the energy of the relativistic flow is dissipated and this 
provides the energy needed for the GRB and the subsequent afterglow. The dissipation 
is in the form of (collisionless) shocks, possibly via plasma instability. There is a general 
agreement that the afterglow is produced via external shocks with the circumburst matter 
see IVIIjl . There is convincing ev iden ce (see e.g. Fenimore et al. [94j . Piran and Nakar 
31l| . Ramirez- Ruiz and Fenimore \%2^ . Sari and Piran and §VI Al below) that in most 
bursts the dissipation during the GRB phase takes place via internal shocks, namely shocks 
within the rela tivist ic flow itself. Some (see e.g. Dar [70j, Dermer and Mitman [78j, Heinz 



and Begelman 



162], Ruffini et al. [358]) disagree with this statement. 
Synchrotron Radiation: Most models (both of the GRB and the afterglow) are based 
on Synchrotron emission from relativistic electrons accelerated within the shocks. There 
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is a reasonab 
observations 



e agreement between the predictions of the synchrotron model and afterglow 



140 



291 



439] . These are also supported by measurements of linear polarization 



in several optical afterglows (see £111 B 2|) . As for the GRB itself there are various worries 



The main alternative to Synchrotron emission is sy nchrotron-self Compton 
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211 



321 


] anc 


ton 
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383, 384]. The last model requires, of course a 



inverse Compton of external light 
reasonable source of external light. 

Jets and Collimation: Monochromatic breaks appear in many afterglow light curves. 
These br eaks are interp reted as "jet breaks" due to the sideways beaming of the relativistic 



emission 



292, 



344, 



374] (when the Lorentz factor drops below 1/Bq the radiation is beamed 



outside of the original je t reducing the observed flux) and due to the sideways spreading of 



a beamed flow 
structured jet" 



344, 



374 
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. An a lternative interpretation is of a viewing angles of a "universal 

"r~i 



3471 . l44q whose energy varies with the angle. Both interpretations 
suggest that GRBs are beamed. However, they give different estimates of the overall rate 
and the energies of GRBs (see §VIIII below). In either case the energy involved with GRBs is 
smaller than the naively interpreted isotropic energy and the rate is higher than the observed 
rate. 

A (Newborn) Compact Object If one accepts the beaming interpret ation of the 



105, 



2911. It 



breaks in the optical light curve the total energy release in GRBs is ~ 10 51 ergs 
is higher if, as some models suggest, the beaming interpretation is wrong or if a significant 
amount of additional energy (which does not contribute to the GRB or to the afterglow) is 
emitted from the source. This energy, ~ 10 51 ergs, is comparable to the energy released in a 
supernovae. It indicates that the process must involve a compact object. No other known 
source can release so much energy within such a short time scale. The process requires a 
dissipation of ~ 0.1m© within the central engine over a period of a few seconds. The sudden 
appearance of so much matter in the vicinity of the compact object suggest a violent process, 
one that most likely involves the birth of the compact object itself. 

Association with Star Formation and SNe: Afterglow observations, which exist for 
a subset of relatively br ight long bursts, show that GRBs arise within galaxies with a high 
star formation rate (see|83j and flll C 1|) . Within the galaxies the bursts distribution follows 



the light distribution 



This has lead to the understanding that (long) GRB arise from 
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the collapse of massive stars (see flIX D|) . This understanding has been confirmed by the 
appearance of SN bumps in the afterglow light curve (see £111 C 41 earlier) and in particular 
by the associations of SN 1999bw with GRB 980425 and of SN 2003dh with GRB 030329. 

Summary: Based on these generally accepted ideas one can sketch the following generic 
GRB model: GRBs are a rare phenomenon observed within star forming regions and associ- 
ated with the death of massive stars and the birth of compact objects. The 7-rays emission 
arises from internal dissipation within a relativistic flow. This takes place at a distances of 
~ 10 13 — 10 15 cm from the central source that produces the relativistic outflow. Subsequent 
dissipation of the remaining energy due to interaction with the surrounding circumburst 
matter produces the afterglow. The nature of the "inner engine" is not resolved yet, how- 
ever, a the association with SN (like 1998bw and 2003dh) shows that long GRBs involve a 
a collapsing star. Much less is known on the origin of short GRBs. 



IV. RELATIVISTIC EFFECTS 



A. Compactness and relativistic motion 

The first theoretical clu es to the necessity of relativistic motion in GRBs arose from 
the Compactness problem jsfSrj l . The conceptual argument is simple. GRBs show a non 
thermal spectrum with a significant high energy tail (see §11 A 1)1 . On the other hand a 
naive calculation implies that the source is optically thick. The fluctuations on a time scale 
St imply that the source is smaller than cSt. Given an observed flux F, a duration T, and 
an distance d we can estimate the energy E at the source. For a typical photon's energy E~ 
this yields a photon density ~ 4ird 2 F/E^c 3 St 2 . Now, two 7-rays can annihilate and produce 
e + e~ pairs, if the energy in their CM frame is larger than 2m e c 2 . The optical depth for pair 
creation is: 

f e ±cr T 47id 2 F 

* E,c 2 St (5) 
where, f e ± is a numerical factor denoting the average probability that photon will collide 
with another photon whose energy is sufficient for pair creation. For typical valu es and 



306j. This 



cosmological distances, the resulting optical depth is extremely large r e ± ~ 10 15 
is, of course, inconsistent with the non-thermal spectrum. 

The compactness problem can be resolved if the emitting matter is moving relativistically 
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towards the observer. I denote the Lorentz factor of the motion by T. Two corrections appear 
in this case. First, the observed photons are blue shifted and therefore, their energy at the 
source frame is lower by a factor T. Second, the implied size of a source moving towards us 
with a Lorentz factor T is cStT 2 (see £11 V Bl below). The first effect modifies f e ± by a factor 
Y~ 2a where a is the photon's index of the observed 7-rays (namely the number of observed 
photons per unit energy is proportional to E~ a .). The second effect modifies the density 
estimate by a factor T~ 4 and it influences the optical depth as T~ 2 . Together one finds that 
for a ~ 2 one needs T > 100 to obtain an optically thin source. 

The requirement that the source would be optically thin can be used to obtain d irect 
limi t s from sp e cific bursts on the minimal Lorentz factor within those bursts [15I . 91 . 1199L 
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306 



308, 



442]. A complete calculation requires a detailed integration over angular 
integrals and over the energy dependent pair production cross section. The minimal Lorentz 
factor depends also on the maxim al photon energy, E m3u ^, the upper energy cutoff of the 



spectrum. Lithwick and Sari 220| provide a detailed comparison of the different calculations 
and point our various flaws in some of the previous estimates. They find that: 

11 <JTd 2 (m e C 2 ) a+ F /-E'niax\Q!_l r 2a+2/ 1 , ^\2a-2 / R \ 

T - = 180 c 2 5T(a-l) W } 1 {1 + Z) ' (6) 
where the high end of the observed photon flux is given by TE~ a (photons per cm 2 per sec 
per unit photon energy). A lower limit on T is obtained by equating Eq. |UI to unity. 

B. Relativistic time effects 

Consider first a source moving relativistically with a constant velocity along a line towards 
the observer and two photons emitted at R\ and R2. The first photon (emitted at R\) will 
reach the observer at time (R 2 — R\)/v — (R2 — R\)/c before the second photon (emitted 
at R2). For T ^> 1 this equals ~ (R2 — Ri)/2cT 2 . This allows us to associate an "observer 
time" R/2cT 2 with the distance R and for this reason I have associated a scale c5tT~ 2 with 
fluctuations on a time scale 5t in the optical depth equation earlier (see £IIV A|) . This last 
relation should be modified if the source moves a varying velocity (v=v(R)). Now 

f R2 dR 



which reduces to 
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i?/2cr 2 , (8) 



for motion with a constant velocity. The difference between a constant velocity source and 
a decelerating sourc e intr oduces a numerical factor of order eight which is important during 



the afterglow phase [363 ] . 



Consider now a relativistically expanding spherical shell, or at least a shell that is locally 
spherical (on a scale larger than 1/T). Emission from parts of the shell moving at angle 9 
relative to the line of sight to the observer will arrive later with a time delay R(l — cos9)/c. 
For small angles this time delay equals R6 2 /2c. As the radiation is beamed with an effective 
beaming angle ~ 1/T most of the radiation will arrive within a typical angular time scale: 

T ang = R/2cT 2 . (9) 

The combination of time delay and blueshift implies that if the emitted spectrum is a power 
law spectrum with a spectral index a then the observed signal from the ins tanta neous 
emission of a thin shell will decay at late time as a power law with t~( 2-a ) j^j [^Z^]. The 
observed pulse from an instantaneous flash from a thin shell is shown in Fig. EI 

As I discuss later (see WI A|) the similarity between the angular time scale and the radial 
time scale plays a crucial role in GRB models. 



C. Relativistic Beaming and the Patchy Shell Model 

The radiation from a relativistic source is beamed with a typical beaming angle 1/T. 
This implies that if the source that is expanding radially with an ultra-relativistic speed a 
given observer "sees" radiation only from a region that is within T _1 from its line of sight 
to the source. If the radius of the emitting region is R the observer will see radiation from 
a region of size R/Y. Since T is extremely large during the GRB we observe emission only 
from a small fraction of the emitting shell. It is possible, and even likely, that the conditions 
within the small region that we observe will be different from the average ones across the 
shell. This means that the conditions that we infer won't reflect the true average conditions 
within this particular GRB. 

An interesting point related to the internal shocks (discussed later) model in this context 
is the following. According to the internal shocks model individual pulses are obtained 
by collisions between individual shells. Here the inhomogeneity of individual shells could 
be wiped out when the contributions of different hot spots from different shells is added. 
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Alternatively the "inner engine" may produce a consistent angular pattern in which the hot 
spot is in the same position in all shells and in this case averaging won't lead to a cancellation 
of the patchy shell structure. 

Within the internal-external model the GRB is produced by internal shocks in which 
only the relative motion within the flow is dissipated. The bulk Lorentz factor remains 
unchanged. During the afterglow the shell is slowed down by external shocks. As the 
Lorentz factor decreases with time (see Eq. 178)1 we observe a larger and larger fraction of 
the emitting region until T m where 9 is the angular size of the whole emitting region - 
the GRB jet, see §VII HI This has several inevitable implications. If the initial relativistic 
flow is inhomogenous on a small angular scale then different observers looking at the same 
GRB (from different viewing angles) will see different 7-rays light curves. A strong burst 
to one observer might look weak to another one if it is located at an angle larger than 1/T 
from the first. The two observers will see similar conditions later on, during the afterglow, 
as then they will observe the same angular regions. This has the following implications: 
(i) Given that the GRB population originate from some 'typical' distribution we expect 
that fluctuation between different bursts at early time d uring the GRB will be larger than 
fluctuations observed at late times during the afterglow 206J. A direct consequence of this 
behaviour is the appearance of a bias in the observations of GRBs. As we are more likely 
to detect stronger events we will tend to identify bursts in which a 'hot spot' was pointing 
towards us during the GRB phase. If the original GRB shells are inhomogenous this would 
inevitably lead to a bias in the estimates of the GRB emission as compared to the kinetic 
energy during the afterglow, (ii) As the afterglow slows down we observe a larger and larger 
region. The angular structure would produces a variability in the light curve with a typical 
time scale of t, the observed time. These fluctuations will decay later as the Lorentz f actor 



decreases and the observations are averaged over a larger viewing angle. Nakar et al. 



274\ 



have suggested that this is th e source of the early fluctuations in the light curve of GRB 
021004. Nakar and Oren |26?| modelled this process with a numerical simulation. They find 
that the flucutation light curve of GRB 021004 can be nicely fitted by this model and that 
it also explains the correlated fluctuations in the polarization (see also 133| ). 
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V. PHYSICAL PROCESSES 



The observed prompt emission must be generated by energetic particles that have been 
accelerated within the collisionless shocks. The most likely process is synchrotron emission, 
even though there is some evidence that a simple synchrotron spectra does not fit all bursts 



321] (but see however, [16] who finds consistency with the synchrotron model). I consider 



here, the different physical ingredient that determine the emission process: particle accel- 
eration, magnetic field amplification, synchrotron emission and inverse Compton emission 
that could be relevant in some cases. 



A. Relativist ic Shocks 

Shocks involve sharp jumps in the physical conditions. Conservation of mass, energy and 
momentum determine the Hugoniot shock iump conditions across the relativistic shocks for 

n 

the case when the upstream matter is cold (see e.g. Blandford and McKee [30]): 

n 2 = 4Tm (10) 
e 2 = 4Tnim p c 2 

r 2 sh = 2r 2 

where ni 2 ,ei^ are the number density and the energy density (measured in the local rest 
frame) of the matter upstream (region 1) and downstream (region 2). I have assumed that 
the energy density in region 1 is very small compared to the rest mass density. V is the 
Lorentz factor of the fluid just behind the shock and T s h is the Lorentz factor of the shock 
front (both measured in the rest frame of the upstream fluid) . The matter is compressed by 
a factor T across a relativistic shock. The pressure, or the internal energy density behind the 
shock is of order T 2 riim p c 2 . Thus, in the shock's rest frame the relative "thermal" energy 
per particle (downstream) is of the same order of the kinetic energy per particle (ahead of 
the shock) upstream. Put differently the shock converts the 'ordered' kinetic energy to a 
comparable random kinetic energy. In an ultra-relativistic shock the downstream random 
velocities are ultra-relativistic. 

Similar jump conditions can be derived for the Magnetic fields across the shock. The 
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parallel magnetic field (parallel to the shock front) Bu is compressed and amplified: 

sii2 = rsin (n) 

The perpendicular magnetic field B± remains unchanged. 

The energy distribution of the (relativistic) electrons and the magnetic field behind the 
shock are needed to calculate the Synchrotron spectrum. In principle these parameters 
should be determined from the microscopic physical processes that take place in the shocks. 
However, it is difficult to estimate them from first principles. Instead I define two dimension- 
less parameters, e# and e e , that incorporate our ignorance and uncertainties [2881 .1305. 367]. 
It is commonly assumed that these energies are a constant fraction of the internal energy 
behind the shock (see however, Daigne and Mochkovitch [68]). I denote by e e and by the 
ratio between these energies and the total internal energy: 

e e = e e e = AT 2 sh e e nim p c 2 (12) 
e B = B 2 /8n = e B e = AT 2 h e B nim p c 2 

One usually assumes that these factors, e ej B, are constant through out the burst evolution. 
One may even expect that they should be constant from one burst to another (as they reflect 
similar underlying physical processes). However, it seems that a simple model that assumes 
that these parameters are constant during the prompt burst cannot reproduce the observed 
spectrum |68j ]. This leads to explorations of models in which the equipartition parameters 
e £t B depend on the physical conditions within the matter. 

In GRBs, as well as in SNRs the shocks are collisionless. The densities are so small so 
that mean free path of the particles for collisions is larger than the typical size of the system. 
However, one expects that ordered or random magnetic fields or alternatively plasma waves 
will replace in these shocks the role of particle collisions. One can generally use in these 
cases the Larmour radius as a typical shock width. A remarkable feature of the above shock 
jump conditions is that as they arise from general conservation laws they are independent of 
the detailed conditions withi n the shocks and hence are expected to hold within collisionless 



shocks as well. See however 
in GRBs. 



263] for a discussion of the conditions for collisionless shocks 
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B. Particle Acceleration 



It is now generally accepted that Cosmic rays (more specifically the lower energy com- 
oonent below 10 15 eV) are accelerated within shocks in SNRs is the Galaxy (see e.g. Gaisser 



117]). A beautiful demonstration of this effect arises in the observation of synchrotron emis- 
sion from Supernova remnants, which shows X-ray emission from these accelerated particles 
within the shocks. 

The common model for particle shock acceleration is the Diffuse Shock Acceleration 
(DSA) model. According to this model the particles are accelerated when they repeatedly 
cross a shock. Magnetic field irregularities kee p sc attering the particles back so that they 
keep crossing the same shock. The competition 96J] between the average energy gain, Ef/E^ 
per shock crossing cycle (upstream-downstream and back) and the escape probability per 
cycle, P esc leads to a power-law spectrum N(E)dE oc E~~ p dE with 

p = 1 + Zn[l/(1 - P esc )]/ln[(E f /E t )]. (13) 

Note that within the particle acceleration literature this index p is usually denoted as s. 
Our notation follows the common notation within the GRB literature. 

Blandford and Eichler j^J review the theory of DSA in non-relativistic shocks. However, 
in GRBs the shocks are relativistic (mildly relativistic in internal shocks and extremely 
relativistic in external s hocks). Acceleratio n in ultra relativistic shocks have been discussed 



by several groups 



17 



121 



161 



190 



424] In relativistic shocks the considerations are quite 



different from those in non-relativistic ones. Using the relativistic shock jump conditions 
(Eq. ^2 an d kin ematic consi derations one can find (see Achterberg et al. [1] , Gallant and 
Achterberg jl22| . Vietri 422 1) that the energy gain in the first shock crossing is of the order 
T 2 sh . However, subseq uent shock crossing are not as efficient and the energy gain is of order 
unity (Ef/Ei) «2 QQ. 

The deflection process in the upstream region is due to a large scale smooth background 
magnetic field perturbed by MHD fluctuations. A tiny change of the particle's momentum 
in the upstream region is sufficient for the shock to overtake the particle. Within the down- 
stream region the momentum change should have a large angle before the particle overtakes 
the shock and reaches the upstream region. As the shock moves with a sub-relativistic ve- 
locity (~ c/v^3) relative to this frame it is easy for a relativistic particle to overtake the 
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shock. A finite fraction of the particles reach the upstream region. Repeated cycles of this 
type (in each one the particles gain a factor of ~ 2 in energy) lead to a power-law spectrum 
with p « 2.2 — 2.3 (for T s /j ^> 1). Like in non-relativistic shock this result it fairly robust and 
it does not depend on specific assumptions on the scattering process. It was obtained by 
several groups using different approaches, including both numerical simulations and analytic 
considerations. The insensitivity of this result arises, naturally from the logarithmic depen- 
dence in equation and from the fact that both the denominator and the numerato r are 
of order unity. This result agre es ni cely with what was inferr ed fr om GRB spectrum [369] 
or with the afterglow spectrum |29 1| . Ostrowski and Bednarz [283^ point out, however, that 
this result requires highly turbulent conditions downstream of the shock. If the turbulence 
is weaker the resulting energy spectrum could be much steeper. Additionally as internal 
shocks are only mildly relativistic the conditions in these shocks might be different. 

The maximal energy that the shock accelerated particles can be obtained by comparing 
the age of the shock R/c (in the upstream frame) with the duration of an acceleration cycle. 
For a simple magnetic deflection, this later time is just half of the Larmour time, E/Zq e B 
(in the same frame). The combination yields: 

E max « Zq e BR = W 20 eVB 3 R 15 , (14) 

where the values that I have used in the last equality reflect the conditions within the reverse 
external shocks where UHECRs (Ultra High Energy Cosmic Rays) can be accelerated (see 
§VIII CI below). For particle diffusion in a random upstream field (with a diffusion length I) 



one finds that R in the above equation is replaced by y Rl/3. 

The acceleration process has to compete with radiation losses of the accelerated parti- 
cles. Synchrotron losses are inevitable as they occur within the same magnetic field that is 
essential for deflecting the particles. Comparing the energy loss rate with the energy gain 
one obtain a maximal energy of: 

E max « mc 2 (^f « 5 • 10 17 eV(m/m p )r^S-V2. (15) 

The corresponding Lorentz factor is of the order of 10 8 for T s h = 100 and B = 1 Gauss. 
Note that this formula assumes that the acceleration time is the Larmour time and hence 
that the synchrotron cooling time is equal to the Larmour time. Obviously it should be 
modified by a numerical factor which is mostly likely of order unity. 
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C. Synchrotron 



Synchrotron radiation play, most likely, an important role in both the GRB and its 
afterglow. An important feature of synchrotron emission is its polarization (see Fj) . 
Observations of polarization in GRB afterglows and in one case in the prompt emission 
support the idea that synchrotron emission is indeed taking place there (note however that 
IC also produces polarized emission). I review here the basic features of synchrotron emis sion 
focusing on aspects relevant to GRBs. I refer the reader to Rybicki and Lightman |36l| for 
a more detailed discussion. 



1. Frequency and Power 

The typical energy of synchrotron photons as well as the synchrotron cooling time depend 
on the Lorentz factor 7 e of the relativistic electron under consideration and on the strength 
of the magnetic field . If the emitting material moves with a Lorentz factor V the photons 
are blue shifted. The characteristic photon energy in the observer frame is given by: 

{hv s yn)obs = — ^7e r ' ( 16 ) 

m e c 

where q e is the electron's charge. 

The power emitted, in the local frame, by a single electron due to synchrotron radiation 

is: 

4 

Psyn = -VtCUbiI , (17) 

where Ub = B 2 /8n = e#e is the magnetic energy density and o~t is the Thompson cross 
section. The cooling time of the electron in the fluid frame is then ^ e m e (? j P . The observed 
cooling time t syn is shorter by a factor of V: 

tsynhe) = ■ (18) 

Aa T U B lei- 

Substituting the value of 7 e from equation HB] into the cooling rate Eq. [TBI one obtains 
the cooling time scale as a function of the observed photon energy: 



3 2ncm e q e _ 1/2 

Since 7 e does not appear explicitly in this equation t syn at a given observed frequency 
is independent of the electrons' energy distribution within the shock. This is provided, of 
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course, that there are electrons with the required 7 e so that there will be emission in the 
frequency considered. As long as there is such an electron the cooling time is "universal". 
This equation shows a characteristic scaling of t syn (u) oc v~ x l 2 . This is not very different 
from the observed relation 5T oc u~ 0A 93J]. However, it is unlikely that cooling and not a 



physical process determines the temporal profile. 

The cooling time calculated above sets a lower limit to the variability time scale of a 
GRB since the burst cannot possibly contain spikes that are shorter than its cooling time. 
Observations of GRBs typically show asymmetric spikes in the intensity variation, where a 
peak generally has a fast rise and a slower decay. A plausible explanation of this observation 
is that the shock heating of the electrons happens rapidly (though episodically), and that the 
rise time of a spike is related to the heating time. The decay time is then set by the cooling, 
so that the widths of spikes directly measure the cooling time. However, it seems that there 
are problems with this simple explanation. First when plugging reasonable parameters one 
finds that the decay time as implied by this equation is too short. Second, if the cooling time 
is long the shocked region would suffer adiabatic losses and this would reduce the efficiency 
of the process. Thus it is unlikely that the pulse shapes can be explained by Synchrotron 
physics alone. 



2. The Optically thin Synchrotron Spectrum 

The instantaneous synchrotron spectrum of a single relativistic electron with an initial 
energy j e m e c 2 is approximately a power law with F v oc z/ 1 / 3 up to iy syn {^ e ) and an exponential 
decay above it. The peak power occurs at v S yn{le)i where it has the approximate value 

« -3*L = ^VB. (20) 

VsynKle) 3g e 

Note that P v>max does not depend on j e , whereas the position of the peak does. 

If the electron is energetic it will cool rapidly until it will reach j e>c , the Lorentz factor 
of an electron that cools on a hydrodynamic time scale. For a rapidly cooling electron we 
have to consider the time integrated spectrum. For an initial Lorentz factor 7 e : F u oc z/ -1 / 2 

for U syn (j e>c ) < V < V syn {-ie). 

To calculate the overall spectrum due to the electrons one need s to integrate over the 



electron's Lorentz factor distribution. I consider first, following 



375], a power-law distribu- 
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tion a power index p and a minimal Lorentz factor r ) e ,min- This is, of course, the simplest 
distribution and as discussed in WBI this is the expected distribution of shock accelerated 
particles: 

iV( 7e )~ 7e -P for 7e > le>min . (21) 

The condition p > 2 isrequired so that the energy does not diverge at large 7e (Bhattacharya 
Dai and Cheng 



64j consider also distributions with 2 > p > 1 with a maximal energy 
cutoff, see below). The minimum Lorentz factor, 7ejm j n , of the distribution is related to the 
electron's energy density e e and the electron's number density n e as: 

7e,mm = T J = 7^' (22) 

p — 1 n e m e c L p — 1 

The minimal Lorentz factor plays an important role as it characterizes the 'typical' electron's 
Lorentz factor and the corresponding 'typical' synchrotron frequency, v m = v S yn{le,min) ■ 
Interestingly the upper energy cutoff (which essentially exists somewhere) does not play a 
critical role in the spectrum for p > 2. Of course it will lead to a high frequency cutoff of 
the spectrum around v ayn that corresponds to this energy. However, quite generally, this 
happens at the high energy tail far from where the peak flux or the peak energy are emitted. 

A simple modification of the above idea arises if only a fraction, £ej_of the electrons is 
accelerated to high energies and the rest of the electrons remain cold 

0,0- 

If a small 

fraction of electrons shares the energy e e then the typical Lorentz factor would be ^ f T 1 7e,mm, 
where j e ,min is calculated from Eq. |2H above. All the corresponding places where 7e ,mm is 
used should be modified according to this factor. At the same time fewer electrons will be 
radiating. This will introduce a factor £ e that should multiply the total emitted flux. In the 
following discussion I wil l not add this factors into the analysis. Similarly in situations when 
multiple pair are formed \yi\ the electron's energy is shared by a larger number of electron. 
In this case £ e is larger than unity and similar modifications of the spectrum applies. 

The lowest part of the spectrum (strictly speaking the lowest part of the optically thin 
spectrum, as at very low frequencies self absorption sets in, see WC3I below) is always 
the sum of the contributio ns of the tails of all the electron's emission: F u oc i/ 1 / 3 . This is 
typical to synchrotro n [s^jjlsJ . l25sj ] and is independent of the exact shape of the electron's 
distribution. Tavani [40 lL |402] , for example obtain such a low energy spectrum both for a 
Gaussian or for a Gaussian and a high energy power-law tail. The observation of bursts 
(about 1/5 of the bursts) with steeper spectrum at the lower energy part, i.e. below the 
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"synchrotron line of death" |32lL |322] is one of the problems that this model faces. The 
problem is even more severe as in order that the GRB will be radiating efficiently, otherwise 
the efficiency will be very low, it m ust b e in the fast cooling regime and the relevant low 



energy spectrum will be oc v 



-1/2 



55, 



125| . However, as stressed earlier (see §refsec:spec-obs) 



this problem is not seen in any of the HETE spectrum whose low energy tail is always in 
the proper synchrotron range with a slope [lfj and it might be an artifact of the low energy 
resolution of BATSE in this energy range 55| . 

On the other hand the most energetic electrons will always be cooling rapidly (indepen- 
dently of the behavior of the "typical electron"). These electrons emit practically all their 
energy m e c 2 7, at their synchrotron frequency. The number of electrons with Lorentz factors 
~ 7 is oc 7 1_p and their energy oc r y 2 ~ p . As these electrons cool, they deposit most of their 
energy into a frequency range ~ v ay n{l) oc 7 2 and therefore F u oc 7~ p oc v~' p l 2 . Thus the 
uppermost part of the spectrum will satisfy: 

F v = N[j(p)]m e c 2 ^r(u)dj/du oc z/~ p/2 . (23) 

In the intermediate frequency region the spectrum differs between a 'slow cooling' if the 
'typical' electrons with 7 e ,min do not cool on a hydrodynamic time scale and 'fast cooling' if 
they do. The critical parameter that determines if the electrons are cooling fast or slow is 
7e )C , the Lorentz factor of an electron that cools on a hydrodynamic time scale. To estimate 
7e iC compare t syn (Eq. [TS]) with t hyd , the hydrodynamic time scale (in the observer's rest 
frame) : 

3m e c 
4aTU B Tt h yd 

For fast cooling 7 e , m m < 7e,c> while 7 e ,mm > 7e,c for slow cooling. In the following discussion 
two important frequencies play a dominant role: 

= ^synipfe,min) j (25) 
= ^syn vTe,c) • 

These are the synchrotron frequencies of electrons with je, min and with 7 e c . 

Fast cooling (7 e c < j ejm i n ): The typical electron is cooling rapidly hence v c < v m . The 
low frequency spectrum F u oc u 1 ^ extends up to v c . In the intermediate range between, v c 
and v m , we observe the energy of all the cooling electrons. The energy of an electron oc 7, 
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and its typical frequency oc y 2 the flux per unit frequency is oc 7 1 oc v l ' 2 . Overall the 
observed flux, F u , is given by: 

v c < v < v m , (26) 

where v m = v syn {y^ min ),v c = v syn (y e . c ) and F v>max is the observed peak flux. The peak 
flux is at v c F v ^ max = N e P Utmax /4:7rD 2 (where D is the distance to the source and I ignore 
cosmological corrections). The power emitted is simply the power given to the electrons, 
that is e e times the power generated by the shock, dE/dt: 

dE . , 

Pf.„ = e,-. (27) 

The peak energy emitted (which corresponds to the peak of vF v ) is at v m . The resulting 
spectrum is shown in Fig. ESI 

Slow cooling (7 e;C > 7 e>m j„): Now only the high energy tail of the distribution (those 
electrons above 7 e>c ) cools efficiently. The electrons with 7 e ~ 7e )m j n , which form the bulk 
of the population, do not cool. Now f v oc u 1 ^ 3 up to u m , and F v oc v~ p / 2 above u c . In the 
intermediate region between these two frequencies: 

F v = N[(y(p)]P[(y(p)}dy/du oc v~^l 2 \ (28) 

where y[y) is the Lorentz factor for which the synchrotron frequency equals u, N[y] is the 
number of electrons with a Lorentz factor 7 and P[y] the power emitted by an electron with 
7. Overall one finds: 

{iy/^m) ^ F umax , V <C V m i 

(v/v m )- {p - 1)/2 F um v m <v< u c , (29) 

{v c /v m y ip ~ l)l2 (u/u c y p/2 F u>max , v c < v. 

The peak flux is at u m while the peak energy emitted is at u c . The emitted power is 
determined by the ability of the electrons to radiate their energy: 

Pslow -^e Pxyn ipfc.min ) (30) 

where, N e is the number of electrons in the emitting region and P S yn{le,min)-, the synchrotron 
power of an electron with 7 ejm j n , is given by Eq. [T7| 
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Typical spectra corresponding to fast and slow cooling are shown in Fig. [231 The 
light curve depends on the hydrodynamic evolution, which in turn determines the time 
dependence of u m , v c and F, / r r , n/r ,. The spectra presented here are composed of broken 



power laws. Granot and Sari 



144j present a more accurate spectra in which the asymp- 



totic power law segments are connected by smooth curves. They fit the transitions by 



/ub) n ^ 2 ] 1 ' n . The parameter n estimates the smoothness of the transition 



with n ~ 1 for all transitions. 

Fast cooling must take place during the GRB itself: the relativistic shocks must emit 
their energy effectively - otherwise there will be a serious inefficiency problem. Additionally 
the pulse won't be variable if the cooling time is too long. The electrons must cool rapidly 
and release all their energy. It is most likely that during the early stages of an external 
shock (that is within the afterglow phase - provid ed that it arises due to external shocks) 



there will be a tran sition from fast to slow cooling 



Tavani 



401 



187 



259 



261 



430 



431|. 



402| discusses the synchrotron spectrum from a Gaussian electron distribu- 
tion and from a Gaussian electron distribution with a high energy tail. As mentioned earlier 
the Gaussian (thermal) distribution has a typical low frequency z/ 1 / 3 spectrum. However, 
as expected, there is a sharp exponential cutoff at high frequencies. Without a high energy 
tail this spectrum does not fit the observed GRB spectr a of most GRBs (see flll A 1|) . Note, 



however, that it may fit a small subgroup with a NHE 



2961. With an electron distribution 



composed of a Gaussian and an added high energy tail the resulting spectra has the typical 
z/ 1 / 3 component and an additional high energy tail which de pends on the electrons power 
law index. Such a spectra fits several observed GRB spectra 

Another variant is the synchrotron spectrum from a power-law electron distribution with 



1 < p < 2 



28, 64]. In this case there must be a high energy cutoff "f e ,max and the 'typical' 



electron's energy corresponds to this upper cutoff. A possible cutoff can arise from Syn- 
chrotron losses at the energy where the acceleration time equals to the energy loss time (see 
e.g. de Jager et al. |7_3j and the discussion in §VBjl : 



4 x 10 7 B- l/2 . 



(31) 



The resulting "typical" Lorentz factor j em i n differs now from the one given by Eq. E21 
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Bhattacharya 28(, Dai and Cheng 64| find that it is replaced with: 



7e,r. 



P 



p — 1 



I t e L le,Max 



~i V(p-i) 



(32) 



The resulting spectrum is now similar to the one obtained for fast or slow cooling with the 
new critical frequencies v m given by plugging the result of Eq. E21 into Eq. 



3. Synchrotron Self- Absorption 



At low frequencies synchrotron self-absorption may take place. It leads to a steep cutoff 
of the low energy spectrum, either as the commonly known z/ 5 / 2 or as v 2 . To estimate 
the self absorption frequency one needs the optical depth along t he li ne of sight. A simple 
approximation is: a' v ,R/T where a! v , is the absorption coefficient |361 |: 



(p + 2) 

,'2 



The self absorption frequency v a satisfies: a' u ,R/T 



(33) 



1. It can be estimates only once we 



have a model for the hydrodynamics and how do R and 7 vary with time 



142, 



43c 



f-ab sorption frequency depends on the electron distribu- 
z/ 5 / 2 when the synchrotron frequency of the electron 



The spectrum below the the se 
tion. One obtains the well known 
emitting the self absorbed radiation is inside the self absorption range. One obtains v 2 if 
the radiation within the self-absorption frequency range is due to the low energy tail of 
electrons that are radiating effectively at higher energies. For th is la t ter case, which is more 
appropriate for GRB afterglow (for slow cooling with v m < u c ) 

F v oc u 2 [k B T e /(Tm p c 2 )]R 2 , 



184. 



187 



258, 



(34) 



where R is the radius of the radiating shell and the factor k B T e / (Tm p c 2 ) describes the degree 
of electron equipartition in the plasma shock-heated to an internal energy per particle m p c 2 
and moving with Lorentz factor 7. 



The situation is slightly different for a shock heated fast cooling i.e. if v c < v m |143 |. 
In this case we expect the electron's distribution to be inhomogeneous, as electrons near 
the shock did not cool yet but electrons further downstream are cool. This leads to a new 
spectral range v sa <v < v sa > with F v oc z/ 11 / 8 (see Fig. 

Synchrotron self-absorption is probably irrelevant during the GRB itself. Note, however, 
that under extreme conditions the self absorption frequency might be in the low X-ray and 
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this may explain the steep low energy spectra seen in some bursts. These extreme conditions 



are needed in order to make the system optically thick to syn chro tron radiation but keeping 
it optically thin to Thompson scattering and pair creation jl43| |. Self absorption appears 



141 184, 



187, 



430, 



regu larly during the afterglow and is observed typically in radio emission 
439]. The expected fast cooling self-absorbed spectra may arise in the early radio afterglow. 
So far it was not observed. 



D. Inverse Compton 

Inverse Compton (IC) scattering may modify our analysis in several ways. IC can influ- 
ence the spectrum even if the sy stem is optically thin (as it must be) to Compton scattering 



(see e.g. Rybicki and Lightman 3611 ] ). In view of the high energies involved a photon is IC 
scattered only once. After a single IC scattering the photon's energy is so high that in the 
electron's rest frame it is above the Klein- Nishina energy (m e c 2 ~ 0.5Mev), and the decrease 
in the Compton cross section in this energy range makes a second scattering unlikely. Note 
that in some cases (e.g. in forward external shocks) even the first scattering may suffer from 
this problem. The effect of IC depends on the Comptonization parameter Y = 7 2 r e . For 
fast cooling one can show [3671 ] that Y satisfies: 

Y = e e /U B if U e <U B (35) 



Y = yJUjU B if U e > U B , 

where U e and U B are the energy densities of the electron's and of the magnetic field respec- 
tively. IC is unimportant if Y < 1 and in this case it can be ignored. 



If Y > 1, which corresponds to U e > U B (or to e e > e B ) and to Y — yU e /U B , then 
a large fraction of the low energy synchrotron radiation will be up scattered by IC and a 
large fraction of the energy will be emitted via the IC processes. Those IC photons might 
be too energetic, that is their energy may be far beyond the observed energy range. In 
this case IC will not influence the observed spectra directly. However, as IC will take a 
significant fraction of the energy of the cooling electrons it will influence the observations 
in two ways: it will shorten the cooling time (the emitting electrons will be cooled by both 
synchrotron and IC process). Second, assuming that the observed 7-ray photons results from 
synchrotron emission, IC will influence the overall energy budget and reduce the efficiency 
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of the production of the observed radiation. I turn now to each of this cases. 

An IC scattering boosts the energy of the photon by a factor 7^. Typical synchroton 
photon that have been scattered once by IC will be observed at the energy: 

(hvic)obs = — ^7e r - ( 36 ) 

m e c 

The electrons are cooled both by synchrotron and by IC. The latter is more efficient and 
the cooling is enhanced by the Compton parameter Y. The cooling time scale is: 



67rc 3 /Y^/^^ 1/4 m3/ 4 g e 1 /4 

tW = £7/4(^)1/^3/4^ ( 37 ) 

The conditions needed to produce the observed emission using IC are pro bably not ful- 
filled in either external or internal shocks (see however Ghisellini and Celotti jl24j and the 



discussion in WEI below). However even if IC does not produce the observed 7-ray photons 
it still influences the process if Y > 1. First it will add an ultra high energy component 
to the GRB spectrum. This component will typically be at around 7g times the observed 



~ lOOKeV photons, namely at the GeV-TeV range (see e.g. Bottcher and Dermer [41], Vi- 
etri |42.^ and the discussion in §VIII Ajl . This component might have been already observed 
in some GRBs during the early afterglow (see flll A 1|) . Inverse Compton will also speed up 
the cooling of the emitting regions and shorten the cooling time, t syn estimated earlier (Eq. 
IT^|) by a factor of Y . At the same time this also reduces the efficiency (for producing the 
observed 7-rays) by the same factor. 

E. Quasi-Thermal Comptonization 



Ghisellini and Celotti |124j suggested that the prompt GRB emission arises in a quasi- 
thermal Comptonization process. In their model the optical depth within the emitting region 
(of internal shocks) is of order unity leading to a copious pair production. The system is 
optically thick to synchrotron emission. The self-absorbed synchrotron emission is the seed 
for an Inverse Compton emission produced by the pairs. The effective Compton paramter 
in the new system, Y, is: 

kT' kT' 

r = 4r(^- 5 )(l + r)[l + 4( ,)], (38) 

m e c z m e c z 

where T' is the effective temperature of the pair and r is the total cross section for scattering. 

The pairs act as a thermostat controlling the effective temperature within the emitting region 
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to 30-300kev 



399 



400]. The resu lting spectrum from this model is a flat spectrum F v oc u° 



between the hu qn T and kT'T 



124J. The spectrum will evolve rapidly during the burst while 



the pairs are being created and the effective temperature decreases. 



F. Polarization from Relativistically Moving Sources 

Polarization can provide information on both the emission process and on the geometry 
of the emitting regions. Usually the observed polarization is obtained by first integrating the 
Stokes parameters of the radiation emitted by the individual electrons over the electron's 
distribution. This yields the local polarization. Then we integrate over the emitting region 
to obtain the global polarization. In GRBs (both in the prompt emission and in the after- 
glow) the emitting regions move relativistically towards the observed. The implied Lorentz 
transformations play a very important role in the second integration as they change the di- 
rection of propagation of the photons and hence the direction of the local polarization. The 
final results are sometimes surprising and counter intuitive. For example even if the intrinsic 
(local) emission is 100% polarized in the same direction the integration over the emitting 
region would reduce this to 70% polarization. I consider polarization from synchrotron 
emission here, but the results can be easily applied to IC as well. I apply the results derived 
in this section to the possible polarization from the prompt emission and from the afterglow 
in the corresponding sections §VIEI and IVII.TI 

As an example I consider synchrotron emission. Synchrotron emission is polarized with 
and the intrinsic local polarization l evel depends on the spectral index of the energy distri- 
bution of the emitting electrons, p, [36 1| . For typical values (2 < p < 3) it can reach 75%. 
The polarization vector is perpendicular to the magnetic field and, of course, to the direction 
of the emitted radiation. The formalism can be easily adopted also to Inverse Compton for 
which the intrinsic local polarization is higher and could reach 100% when the photons are 
scattered at 90°. 

Consider first a case where the magnetic field is uniform locally (over a regions of angular 
size r _1 ). This could happen, for example, if we have an ordered magnetic field along the 
(f) direction and the observer is more than away from the symmetry axis. This would 
be the case within internal shocks if the magnetic field is dragged from the source or within 
several Poynting flux dominated models. The locally emitted polarization is uniform and 
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is in the plane of the sky and perpendicular to the direction of the magnetic field. In a 
Newtonian system it would combine so that the observed polarization equals the emitted 
one. Howev er, the Lo rentz transformations induce their own signature on the observed 



polarization 



133, 



134]. This is depicted in Fig. It is clear from this figure that the 
polarization vector varies along the observed region (whose angular size is 1/T. Consequently 
the observed global polarization will be smaller than the local polarization. 

The observed stokes parameters are weighted averages of the local stokes parameters 
at different regions of the shell. The instantaneous polarization is calculated using the 
instantaneous observed flux F u (y,T) oc (1 + t/)~( 3+Q \ with a the relevant spectral index 
at this segment, as the weights, where y = (T6) 2 and T is the observer time. The time 
integrated polarization is calculated using the fluences as weights: J °° F v (y,T)dT oc (1 + 

The fluxes depend on how the intensity varies with the magnetic field. For l v oc B°, which 
is relevant for fast cooling 5 (and the prompt GRB), the time integrated stokes parameters 
(note that V = as the polarization is linear) and polarization are given by: 

_>_ _ n I shl ^) J m) 

I SmCh r/o°°(l+l/)- (2+a) « ' ^ } 

and the relative polarization is given by 




n = ^±^, (40) 



where 6 P = <f) + arctan(±^J cot 0) (l^ (see also |228|). For a = 1 Eqs. 123101 yield a 
polarization level of U/U sync h ~ 60%. I.e. 60% of the maximal synchrotron polarization, or 
an overall polarization of ~ 45%. Taking the exact values of a and the de pendence of I v on 



134, 



225 1 ■ 



B for fast cooling and p = 2.5 results in an overall polarization of ~ 50% 

It turns out that on e can get a polarized emission even from random magnetic field 
Gruzinov and Waxman jl48| | and Medvedev and Loeb |255| . This happens if the system 
has non spherical geometry. Consider a two dimensional random magnetic field which is in 
the plane of the shock and assume that the correlation length of this magnetic field is very 
short compared to all other length scales in the system. The Lorentz transformation induce 



5 see however Granot Il33 
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in this case a radial polarization pattern going out from the center (where the velocity of 
the matter is towards the observer and the polarization vanishes). This polarization pattern 
is shown in Fig. It is clear that a simple integration over this pattern will lead to a 
vanishing polarization. 

However, a net polarization can arise in several cases if the overall symmetry is broken. 
Polarization will arise if (see Fig. 1160 : 

• We observe a jet in an angle so that only a part of the jet is within an angle of T -1 . 

• If the emission is nonuniform and there are stronger patches with angular size smaller 
than r~ 1 from which most of the emission arise. 

• We observe a standard jet whose emission is angle dependent and this dependence is 
of the order of 



i i r~i . r~i 



Ghisellini and Lazzati [126], Gruzinov [147| . Sari [3651 ] . Waxman 432] suggested that 
polarization can arise from a jet even if the magnetic field is random. Nakar et al. |275 | 
considered a random magnetic field that remains planner in the plane of the shock (for a 
three dimensional random magnetic field the polarization essentially vanishes). For I v oc B° 
the degree of observed polarization of the emission emitted from a small region at angle y 
is: n(y) /II S y n ch = min(y, 1/y). The overall time integrated stokes parameters are: 



Q 
u 



Cfo 00 P»>, m a + y)- {2+a) ™HyMy){ c ° s(20) \'iy<io 

sin(20) 



i Usynch !?s?nji+v)-^*v*i> ' (41) 

where P' v , m = P' v , m (y, 0) is the emitted power at the synchrotron frequency in the fluid rest 
frame. For a top-hat jet with sharp edges P' v , m is constant for any y and </> within the jet 
and zero otherwise. For a structured jet P' v , m depends on the angle from the jet axis. 

The maximal polarization is observed when one sees the edge of the jet. The probability 
to see the edge of a top-hat jet with sharp edges and an opening angle 6jT ^> 1 is negligible. 
On the other hand a jet with 6jT <§C 1 is not expected. Thus the only physical cases in which 
we can expect a large polarization are 1 < OjY < a few . 

Fig. IT7I depicts the time integrated polarization and the efficiency from sharp edged jets 
with different opening angles as a function of the angle between the jet axis and the line of 
sight, 9 Q b s . The efficiency, e// is defined to be the ratio between the observed fluence at 9 b s 
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and the maximal possible observed fluence at 9 b s = 0. In all these cases the polarization 
is peaked above 40%, however the efficiency decrease sharply as the polarization increase. 
Thus the probability to see high polarization grows when 9j decrease. The probability that 
Oobs is such that the polarization is larger than 30% (-Usynch) while eff > 0.1 is 0.68, 0.41, 
0.2 & 0.08 for 9jT = 0.5,1,2,4 respectively. In reality this probability will be smaller, as 
the chance to observe a burst increases with its observed fl ux. 

These later calculations also apply for IC emission j^, I21J ]. However, in this case the 
intrinsic local polarization is around 100% and hence one can reach a maximal polarization 
of ~ 70%. 

Polarization could also arise if the magnetic field is uniform over random patches within 
a region of size T -1 . Here it is difficult, of c ourse to estimate the total polarization without 
a detailed model of the structure of the jet jl48| . 



VI. THE GRB AND THE PROMPT EMISSION 



I turn now to discussion of the theory of the GRB and the prompt emission. It is 
generally accepted that both the GRB and the afterglow arise due to dissipation of the 
kinetic e nergy of the re lativistic flow. The rela t ivist i c m otion can be dissipated by either 



external 



184, 



234, 



333] or internal shocks 



276. 2 



3341 ] . The first involve slowing down 



by the external medium surrounding the burst. This would be the analogue of a supernova 
remnant in which the ejecta is slowed down by the surrounding ISM. Like in SNRs external 
shocks can dissipate all the kinetic energy of the relativistic flow. On the other hand internal 
shocks are shocks within the flow itself. These take place when faster moving matter takes 
over a slower movi ng sh ell. 

Sari and Piran 13701 have shown that external shocks cannot pro duce variable bursts 

n ri 

(see also Fenimore et al. [94]). By variable I mean here, following 370J that St <C T, 
where T is the overall duration of the burst (e.g. T 90 ) and St i s the duration of a typical 
pulse (see ffll A2|) . As most GRBs are varia ble Sari and Piran |37Cl| concluded that most 
GRBs are produced by internal shocks 



334]. Internal shocks can dissipate only a fraction 
of the kinetic energy. Therefore, they must be accompanied by external shocks that f ollow 
and dissipate the remaining energy. This leads to the internal-external shocks scenario |314| . 
GRBs are produced by internal shocks within a relativistic flow. Subsequent external shocks 
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between the flow and the circum-burst medium produce a smooth long lasting emission - the 
afterglow. Various observations (see £111 A6|) support this picture. I begin with the discussion 
with a comparison of internal vs. external shocks. I review then the prompt emission from 
internal shocks, then the prompt emission from external shocks (which includes contributions 
to the late part of long GRBs and the prompt optical flash). I also discuss the transition 
from the observations of one shock to the other. 



Internal vs. External Shocks 



1. General Considerations 



Consider a "quasi" spherical relativistic emitting shell with a radius R, a width A and 
a Lorentz factor T. This can be a whole spherical shell or a spherical like section of a jet 
whose opening angle 9 is larger than Because of relativistic beaming an observer would 
observe radiation only from a region of angular size ~ r . Consider now photons emitted at 
different points along the shock (see Fig. IT5j) . Photons emitted by matter moving directly 
towards the observer (point A in Fig. ITSj) will arrive first. Photons emitted by matter 
moving at an angle (point D in Fig. ITS)) would arrive after t ang = R/2cT 2 . This is also, 
tft, the time of arrival of photons emitted by matter moving direc tly towards the observer 
but emitted at 2R (point C in Fig. fT%|) . Thus, tn ~ t ang |&u . I37J ] . This coincidence is the 
first part of the argument that rules out external shocks in variable GRBs. 

At a given point particles are continuously accelerated and emit radiation as long as the 
shell with a width A is crossing this point. The photons emitted at the front of this shell 
will reach the observer a time t& = A/c before those emitted from the rear (point B in Fig. 
ITSj) . In fact photons are emitted slightly longer as it takes some time for the accelerated 
electrons to cool. However, for m ost reasonable parameters the cooling time is much shorter 
from the other time scales and I ignore it hereafter. 

The emission from different angular points smoothes the signal on a time scale t ang . If 
t/± < t ang ~ tR the resulting burst will be smooth with a width t ang ~ t^. The second part 



of this argument follows 
(see also Sari and Piran 



rom the hydrodynamics of external shocks. I show later in §VI CI 



3701 ]) that for external shocks A/c < R/cT 2 pa tn ~ t ang and for a 



spreading shell A pa R/cT 2 . Therefore external shocks can produce only smooth bursts! 
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As we find only two time scales and as the emission is smoothed over a time scale t ang , a 
necessary condition for the production of a variable light curve is that t& = A/c > t ang . In 
this case t& would be the duration of the burst and t ang the variability time scale. This can 
be easily satisfied within internal shocks (see Fig lliJj) . Consider an "inner engine" emitting 
a relativistic wind active over a time = A/c (A is the overall width of the flow in the 
observer frame). The source is variable on a scale L/c. Internal shocks will take place at 
R s ~ LT 2 . At this place the angular time and the radial time satisfy: t ang « tf> w L/c. 
Internal shocks continue as long as the source is active, thus the overall observed duration 
T = tA reflects the time that the "inner engine" is active. Note that now t ang ~ L/c < t& 
is trivially satisfied. The observed variability time scale in the light curve, St, reflects the 
variability of the source L/c. While the overall duration of the burst reflects the overall 
duration of the activity of the "inner engine" . 

Numerical simulations 193^ have shown that not only the time scales are preserved but 
the source's temporal behavior is repro duced on an almost one to one basis in the observed 
light curve. This can be explained now [268^ by a simple toy model (see WIB 31 below). 



2. Caveats and Complications 

Clearly the way to get around the previous argument is if t ang < £r. In this case one 
can identify tn with the duration of the burst and t ang as the variability time scale. The 
observed variability would require in this case that: t ang /tn = 5t/T. For this the emitting 
regions must be smaller than R/T . 

One can imagine an inhomogenous external medium which is clumpy on a scale d <C R/T 
(see Fig EH)). Consider such a clump located at an angle 9 ~ to the direction of motion 
of the matter towards the observer. The resulting angular time, which is the difference in 
arrival time between the first and the last photons emitted from this clump would be:~ d/cT. 
Now tam ~ d /cT < tR and it seems that one can get around the argument presented before. 

However, Sari and Piran |37C| have shown that such a configuration would be extremely 
inefficient. This third part of this argument rules out this caveat. The observations limit 
the size of the clumps to d ~ cT5t and the location of the shock to R ~ cTT 2 . The number 
of clumps within the observed angular cone with an opening angle T _1 equals the number 
of pulses which is of the order T/5t. The covering factor of the clumps can be directly 
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estimated in terms of the observed parameters by multiplying the number of clumps (T/ St) 
times their area d 2 = (StT) 2 and dividing by the cross section of the cone (R/T) 2 . The 
resulting covering factor equals St/T <C 1. The efficiency of conversion of kinetic energy to 
7-rays in this scenario is smaller than this covering factor which for a typical variable burst 
could be smaller than 10~ 2 . 

I turn now to several attempts to find a way around this result. I will not discuss here the 
feasibility of the suggested models (namely is it likely that the surrounding matte r will be 
clumpy on the needed length scale [78(, or can an inner engine eject "bullets" 162 1 with an 
angular width of ~ 10~ 2 degrees and what keeps these bullets so small even when they are 
shocked and heated). I examine only the question whether the observed temporal structure 
can arise within these models. 



3. External Shocks on a Clumpy Medium sec: Clumpy 

Dermer and Mitman jzs| claim that the simple efficiency argument of Sari and Piran j3TCl| 
was flawed. They point out that if the direction of motion of a specific blob is almost exactly 
towards the observer the corresponding angular time will be of order d 2 /cR and not d/cT 
used for a "generic" blob. This is narrower by a factor dT/R than the angular time across 
the same blob that is located at a typical angle of T -1 . These special blobs would produce 
strong narrow peaks and will form a small region along a narrow cone with a larger covering 
factor. Dermer and Mitman [3| present a numerical simulation of light curves produced by 
external shocks on a clumpy inhomogeneous medium with 5t/T ~ 10~ 2 and efficiency of up 
to ~ 10%. 

A detailed analysis of the light curve poses, however, several problems for this model. 
While this result is marginal for bursts with d~t/T ~ 10~ 2 with a modulation of 50% it is 



insufficient for bursts with 5t/T ~ 10 3 or i 
time scale of milliseconds has been observed 



the modulation is ~ 100%. Variability on a 
269] in many long GRBs (namely 5t/T can be 



as small as 10~ 4 .). Moreover, in this case one would expect that earlier pulses (that arise 
from blobs along the directi on of motion) would be narrower than latter pulses. This is not 
seen in the observed bursts |328^ . 

Finally the arrival time of individual pulses depends on the position of the emitting clumps 
relative to the observers. Two following pulses would arise from two different clumps that 
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are rather distant from each other. There is no reason why the pulses and intervals should 
be correlated in any w ay. R ecall ( §11 A 2)1 that the duration of a pulse and the subsequent 
interval are correlated 270]. 



4- The Shot- Gun Model 



Heinz and Begelman |162j ] suggested that the "inner engine" operates as a shot-gun 
emitting multiple narrow bullets with an angular size much smaller than T _1 (see Fig l2~T)) . 
These bullets do not spread while propagating and they are slowed down rapidly by an 
external shock with a very dense circumburst matter. The pulses width is given by t ang or 
by the slowing down time. The duration of the burst is determined by the time that the 
"inner engine" emits the bullets. 

This model can produce the observed variability and like in the internal shocks model the 
observed light curve represents also here the temporal activity of the source. However, in 
this model the width of the pulses is determined by the angular time or the hydrodynamic 
time or the cooling time of the shocked material. On the other hand the intervals between 
the pulses depend only on the activity of the inner engine. Again, there is no reason why 
the two distribu tions will be similar and why there should be a correlation between them 
(see gHX2]and 27fi| )- 

5. Relativistic Turbulence 

An inter esting alternative t o shocks as a way to dissipate kinetic energy is within plasma 



turbulence 



226 



387 



J. It has been suggested that in this case the kinetic energy of 
the shock is dissipated downstream to a combination of macroscopic (relativistic) random 
motion of plasma blobs with a Lorentz factor IV Within these blobs the particles have also 
a (relativistic) random velocity with a Lorentz factor T p , such that: T s w TbTp. 

Relativistic turbulence may be the only way to produce variability in a situation that the 
matter is slowed down by the external medium and not by internal interaction. I stress that 
in this case the process is not described by regular shocks and hence some of the previous 
arguments do not hold. Two crucial open questions are i) Whether one can produce the 
observed correlations between pulses and intervals, ii) Why there is no spreading of pulses 
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at later times, as would be expected if the emitting region is slowing down and increasing 
its radius. 



B. Internal Shocks 



1. Hydrodynamics of Internal Shocks 



Internal shocks take place when a faster shell catches a slower one, namely at: 



hit 



cStT 2 = 3 x 10 14 cmr? m <5t 



(42) 



where r 10 o is the typical Lorentz factor in units of 10 2 and St is the time difference between 
the emission of the two shells. I show later that St defined here is roughly equal to the 
observed fluctuations in the light curve of the burst St. Clearly R int < R ext must hold oth- 
erwise internal shocks won't take place. R ext is defined as the location of efficient extraction 
of energy by external shocks (see £IV1 Cj) . If follows from the discussion in £IV1 CI that the 
condition R int < R ex t implies: 



where I is defined by Eq. EHJand it is typically of the order of 10 18 cm, while A is the width 
of the shell and it is of order 10 12 cm. Both conditions set upper limits on T (of the order of 
a few thousands) for internal shocks. If the initial Lorentz factor is too large then internal 
shocks will take place at large radii and external shocks will take place before the internal 
shocks could take place. It is possible that this fact plays an important role in limiting the 
relevant Lorentz factors and hence the range of variability of E p , the peak energy observed 
in GRBs. 

Internal shocks are characterized by a comparable Lorentz factor of order of a few (1 < 
T < 10) reflecting the relative motion of the shells and by comparable densities n in both 
shells. In this case, for an adiabatic index (4/3), the Loretz factor of the shocked region T 
satisfies: 



5r 2 <max( ' Z^A 1 / 4 ) 



(43) 




(44) 



The shocked density h and energy e are: 



n = (4r + 3)n » 4LYi ; e 




(45) 
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Both shocks are mildly relativistic and their strength depends on the relative Lorentz factors 
of the two shells. 



2. The Efficiency of Internal Shocks 

Consider collision between two shells with masses m r and m s that are moving at different 
relativistic velocities: r r ^ Y s 3> 1. The resulting bulk Lorentz factor, T m in an elastic 
collision is: 

(46) 



y m r /r r + m s /r s 

The internal energy, £i nt , in the local frame and Ei nt , in the frame of an external observer, 
of the merged shell: E int = T m £ int , is the difference of the kinetic energies before and after 
the collision: 

E int = m r c 2 (T r - r m ) + m s c 2 (T s - T m ). (47) 
The conversion efficiency of kinetic energy into internal energy is jl93^ : 

e = 1 _ (m r + m s )Y m ^ 
(m r T r + m s T s ) 

As can be expected a conversion of a significant fraction of the initial kinetic energy to 



internal energy requires that the difference in velocities between the shel 



T r 3> T s and that the two masses will be comparable m r ~ m s 6tI . Il93 |. 

Beloborodov |l9| considered internal shocks between shells with a lognormal distribution 
of (r — 1)/ (Tq — 1), where To is the average Lorentz factor. The dimensionless parameter, A, 
measures the width of the distribution. He shows that the efficiency increases and reached 
unity when A is of order unity, that is typical fluctuation in T ar e by a factor of 10 compared 
to the average. Similarly numerical simulations of Guetta et al. [154] show that a significant 
fraction of the wind kinetic energy, on the order of 20%, can be converted to radiation, 
provided the distribution of Lorentz factors within the wind has a large variance and the 
minimum Lorentz factor is greater than w lO 2 ' 5 /^ 9 , where L52 is the (isotropic) wind 
luminosity in units of 10 52 ergs/sec. 

Another problem that involves the efficiency of GRBs is that not all the internal energy 
generated is emitted. This depends further on e e , the fraction of energy given to the electron. 
If this fraction is small and if there is no strong coupling between the electrons and the 
protons the thermal energy of the shocked particles (which is stored in this case mostly 
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s will be significant: 



in the protons) will not be radiated away. Inst ead it will be converted again to kinetic 
energy by adiabatic cooling. Kobayashi and Sari jl95^ consider a more elaborated model in 
which colliding shells that do not emit all their internal energy are reflected from each other, 
causing subsequent collisions and thereby allowing more energy to be emitted. In this case 
more energy is eventually emitted than what would have been emitted if we considered only 
the first collision. They obtain about 60% overall efficiency even if the fraction of energy 
that goes to electrons is small e e = 0.1. This is provided that the shells' Lorentz factor varies 
between 10 and 10 4 . 



3. Light Curves from Internal Shocks 



Both the similarity between the pulse width and the pulse 



s epar ation distribution and 



the correlation between int erval s and the subsequent pulses |270l 13271 ] arise naturally within 
the internal shocks model j^sj]. In this model both the pulse duration and the separation 
between the pulses are determined by the same parameter - the i nterv al between the emitted 
shells. I outline here the main argument (see Nakar and Piran \26^ for details). Consider 
two shells with a separation L. The Lorentz factor of the slower outer shell is Ts = T and 
of the Lorentz factor inner faster shell is Tp = aT (a > 2 for an efficient collision) . Both are 
measured in the observer frame. The shells are ejected at t± and t 2 ~ ti + L/c. The collision 
takes place at a radius R s ~ 2r 2 L (Note that R s does not depend on T 2 ). The emitted 
photons from the collision will reach the observer at time (omitting the photons flight time, 
and assuming transparent shells): 

t Q « ti + R s /(2cT 2 ) « tt + L/c « t 2 . (49) 
The photons from this pulse are observed almost simultaneously with a (hypothetical) pho- 



ton that was emitted from the "inner engine" toget 



explains why various numerical simulations 
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193 



r er w ith the second shell (at t-i). This 
2931 find that for internal shocks the 



observed light curve replicates the temporal activity of the source. 

In order to determine the time between the bursts we should consider multiple collisions. 
It turns out that there are just three types of collisions, (i), (ii) and (iii), that characterize 
the system and all combinations of multiple collisions can be divided to these three types. 
Consider four shells emitted at times U (i = 1,2,3,4) with a separation of the order of L 
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between them. In type (i) there are two collisions - between the first and the second shells 
and between the third and the fourth shells. The first collision will be observed at ti while 
the second one will be observed at t^. Therefore, At ~ £4 — ti ~ 2L/c. A different collision 
scenario (ii) occurs if the second and the first shells collide, and afterward the third shell 
takes over and collide with them (the forth shell does not play any roll in this case). The 
first collision will be observed at £2 while th e sec ond one will be observed at £3. Therefore, 
At ~ £3 — £2 ~ L/c. Numerical simulations |268| show that more then 80% of the efficient 
collisions follows one of these two scenarios ((i) or (ii)). Therefore one can estimate: 



Note that this result is independent of the shells' masses. 

A third type of a multiple collision (iii) arises if the third shell collides first with the 
second shell. Then the merged shell will collide with the first one (again the fourth shell 
does not participate in this scenario). In this case the two pulses merge and will arrive 
almost simultaneously, at the same time with a (hypothetical) photon that would have been 
emitted from the inner engine simultaneously with the third (fastest) shell, t ~ £3. Only a 
20% fraction exhibits this type of collision. 

The pulse width is determined by the angular time (ignoring the cooling time): 5t = 
R s / {2cT 2 s ) where T s is the Lorentz factor of the shocked emitting region. If the shells have 
an equal mass (mi = m 2 ) then T s = yfaT while if they have equal energy (mi = am 2 ) then 
r, « r. Therefore: 



The ratio of the Lorentz factors a, determines the collision's efficiency For efficient collision 
the variations in the shells Lorentz factor (and therefore a) must be large. 

It follows from Eqs. EU] and that for equal energy shells the At-8t similarity and cor- 
relation arises naturally from the reflection of the shells initial separation in both variables. 
However, for equal mass shells St is shorter by a factor of a than At. This shortens the 
pulses relative to the intervals. Additionally, the large variance of a would wipe off the 
At-St correlation. This suggests that equal energy shells are more likely to produce the 
observed light curves. 



At « L/c . 



(50) 




equal mass, 
equal energy. 



(51) 
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C. External Shocks 



1. Hydrodynamics 

Consider the situation when a cold relativistic shell (whose internal energy is negligible 
compared to the rest mass) moves into the cold ISM. Generally, two shocks form: an outgoing 
shock that propagates into the ISM or into the external shell, and a reverse shock that 
propagates into the inner shell, with a contact discontinuity between the shocked material 
(see Fig. |HJ). 

There dual shocks system is divided to four distinct regions (see Fig. the ambient 
matter at rest (denoted by the subscript 1), the shocked ambient matter which has passed 
through the forward shock (subscript 2 or f), the shocked shell material which has passed 
through the reverse shock (subscript 3 or r), and the unshocked material of the shell (sub- 
script 4). The nature of the emitted radiation and the efficiency of the cooling processes 
depend on the conditions in the shocked regions 2 and 3. Both regions have the same energy 
density e. The particle densities n 2 and n 3 are, however, different and hence the effective 
"temperatures," i.e. the mean Lorentz factors of the random motions of the shocked protons 
and electrons, are different. 

Two quantities determine the shocks' structure: T, the Lorentz factor of the motion of the 
inner expanding matter (denoted 4) relative to the outer matter (the ISM or the outer shell 
in the case of internal collisions - denoted 1) , and the ratio between the particle number 
densities in these regions, n±jn\. Initially the density contrast between the spherically 
expanding shell and the ISM is large. Specifically n 4 /ni > T 2 . This happens during the 
early phase of an external shock when the shell is small and dense. This configuration 
is denoted "Newtonian" because the reverse shock is non-relativistic at most (or mildly 
relativistic). In this case all the energy conversion takes place in the forward shock. Only a 
negligible fr actio n of the energy is converted to thermal energy in the reverse shock if it is 
Newtonian [368^ . Let T2 be the Lorentz factor of the motion of the shocked fluid relative to 
the rest frame of the fluid at 1 and let F 3 be the Lorentz factor of the motion of this fluid 
relative to the rest frame of the relativistic shell (4): 

r 2 ^r ; f 3 «l. (52) 
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The particle and energy densities (n, e) in the shocked regions satisfy: 

n 2 ~ 4rni, ; e = e 2 = 4T 2 nim p c 2 ; n 3 = 7n 4 , ; e 3 = e. (53) 

Later, the shell expands and the density ratio decreases (like R~ 2 if the width of the shell 
is constant and like R~ 3 if the shell is spreading) and n^/rii < V 2 (but n^jn\ > 1). In this 
case both the forward and the reverse shocks are relativistic. The shock equations bet ween 
regions 1 and 2 coined with the contact discontinnity between 3 and 2 y ie!d UQU 



T 2 = {n^/mf^T 112 /^ ■ n 2 = 4r 2 ni ; e = e 2 = 4T 2 2 n imp c 2 , (54) 
Similar relations hold for the reverse shock: 

F 3 = {n A / ni )- x ^ 2 /V2 ; n 3 = 4T 3 n 4 . (55) 

Additinally 

e 3 = e ; f 3 = (T/r 2 + T 2 /T)/2 , (56) 

which follow from the equality of pressures and velocity on the contact discontinuity. Com- 
parable amounts of energy are converted to thermal energy in both shocks when both shocks 
are relativistic. 

The interaction between a relativistic flow and an external medium depends on the Sedov 
length that is defined generally as: 

rl 

E = rripC 2 / 4:Tm(r)r 2 dr . (57) 

J 

The rest mass energy within the Sedov sphere equals the energy of the explosion. For a 
homogeneous ISM: 

I = {— -)V3 „ I0 18 cm4 2 /3 nl /3 • (58) 

Note that in this section E stands for the isotropic equivalent energy. Because of the very 
large Lorentz factor angular structure on a scale larger than T _1 does not influence the 
evolution of the system and it behaves as if it is a part of a spherical system. A second 
length scale that appears in the problem is A, the width of the relativistic shell in the 
observer's rest frame. 

Initially the reverse shocks is Newtonian and only a negligible amount of energy is ex- 
tracted from the shell. At this stage the whole shell acts "together". Half of the shell's 
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kinetic energy is converted to therm al energy when the collected external mass is M/Y, 



where M is the shell's mass 



184. 



333J. This takes place at a distance: 



Rr = ^ = ( ^V /3 = 5.4 x 10" cm E^n, ^of, (59) 



T 2/3 \n ism m p c 2 T 2 

where E 52 is the equivalent isotropic energy in 10 52 ergs, n\ = ni sm /l particle/cm 3 . 

However, the reverse shock might become relativistic before Rr- Now energy extraction 
from the shell is efficient and one passage of the reverse shock through the shell is sufficient 
for complete conversion of the shell's energy to thermal energy. The energy of the shell 
will be extracted during a single passage of the reverse shock across the shell. Using the 
expression for the velocity of the reverse shock into th e sh ell (Eq. one finds that the 
reverse shock reaches the inner edge of the shell at R& \36^ : 

R A = /^A 1 / 4 « 10 15 cm/ 3 / 4 A}/ 4 . (60) 

The reverse shock becomes relativistic at R^, where 714/ni = T 2 : 

R N = l^/A^T 2 (61) 

Clearly, if R^ > Rr then the energy of the shell is dissipated while the shocks are still 
"Newtonain". If R^ < Rr the reverse shock becomes relativistic. In this case Rr looses its 
meaning as the radius where the energy is dissipated. The energy of the shell is dissipated in 
this "relativistic" ca se at r&. The question which of the two conditions is relevant depends 
on the parameter £ j.36 jj ] : 

£ = (//A) 1/2 r~ 4/3 = 2(/ 18 /A 12 ) 1 / 2 r i0 4 /3 . (62) 

I have used a canonical value for A as 10 12 cm. It will be shown later that within the internal- 
external scenario A/c corresponds to the duration of the bursts and 10 12 cm corresponds to 
a typical burst of 30sec. 

Using £ one can express the different radii as: 

Rintic = ra/^ /2 = = R N /e . (63) 

For completeness I have added to this equation Ri n t, where internal shocks take place (see 

Eq. 14*2*1) . The dimensionless quantity ( : ( = 5/A. Thus: 

( Ra < R-r < Rn £ > 1 (Newtonian reverse shock) 
I Rn < Rr < Ra £ < 1 (Relativistic reverse shock). 
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I have marked in bold face the location where the effective energy extraction does take place. 
With typical values for I, A and r £ is around unity. The radius where energy extraction 
takes place is marked in bold face! 

Expanding s hell : A physical shell is expected to expand during as it propagates with 
A = A + RT 2 315| . This will lead to a monotoneously decreasing £. As the value of 
Rr is independent of A it does not vary. However, Ra and Rn decrease from their initial 
values. If A < RrT 2 (corresponding to £0 > 1) then £ = 1 at R& = Rr = Rn and all three 
radii coincide. Given the fact that with typical parameters £ is of order unity this seems to 
be the "typical" case. The reverse shocks becomes mildly relativistic just when the energy 
extraction becomes efficient. However, if £ ^ 1 then the shell won't expand enough and 
still there will be a relativistic reverse shock operating at Ra- It is useful to note that in 
this case the effective energy extraction takes place at Ra for all initial values of £0. in the 
following I denote by £ the value of £ at Ra'- £~£oif£o<l and otherwise £ ~ 1. 

Overall the external shocks take place at: 

!max(//r 2//3 , / 3//4 A 1//4 ), Non spreading shell, 

I/T 2 / 3 « I 3 / 4 A 1 / 4 » 5 x lO^cmK^ 3 n} /3 rYoo /3 , Spreading shell. ' ° " ' 

Usually I will use the second relation (the spreading shell one) in the following discussion. 
Note that in the case of non spreading shell one uses the maximum of the two possible radii. 
For example in the Newtonian case where the extraction is at I/T 2 / 3 the shocks pass the 
shall many times and hence I /Gamma 2 / 3 > Z^A 1 / 4 . 



2. Synchrotron Spectrum from External Shocks 

The bulk of the kinetic energy of the shell is converted to thermal energy via the two 
shocks at around the time the shell has expanded to the radius Ra (this would be the case 
in either a thick shell with £ < 1 or with an expanding shell that begins with £ > 1 but 
reaches £ ~ 1 due to expansion of the shell around the time when R^ = Ra and efficient 
dissipation takes place . At this radius, the conditions at the forward shock are: 

r 2 = T£ 3/4 , n 2 = 4r 2 m, e 2 = \Y\n x m v c 2 , (66) 

while at the reverse shock: 

r 3 = r 3/4 , r 3 = r£ 3 / 4 , n 3 = ^ 4 rV, e 3 = e 2 . (67) 
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Substitution of T sh = T 2 = T£ 3 / 4 in Eq. [T3] yields, for the equipartition magnetic field: 

B = Vs^ce^Te^m 1 / 2 ^ 2 = (40 G) ef £ 3/4 r 100 ^ /2 . (68) 

If the magnetic field in region 2 behind the forward shock is obtained purely by shock com- 
pression of the ISM field, the field would be very weak, with e# <C 1. Such low fields are 
incompatible with observations of GRBs. I consider, therefore, the possibility that there 
may be some kind of a turbulent instability which brings the magnetic field to approximate 
equipartition In the case of the reverse shock, i.e. in region 3, magnetic fields 

of considerable strength might be present in the pre-shock shell material if the original ex- 
ploding fireball was magnetic. The exact nature of ma gnetic field evolution during fireball 
expansion depends on several assumptions. Thompson |406| found that the magnetic field 
will r emai n in equipartition if it started off originally in equipartition. Meszaros, Laguna & 
Rees [2 57j | on the other hand, estimated that if the magnetic field was initially in equiparti- 
tion then it would be below equipartition by a factor of 10 -5 by the time the shell expands 
to i?A- It is uncertain which, if either one, of the estimates is right. As in the forward shock, 
an instability could boost the field back to equipartition. Thus, while both shocks may have 

< 1 with pure flux freezing, both could achieve e# — > 1 in the presence of instabilities. 
In principle, e# could be different for the two shocks. For simplicity I will consider the same 
value in the following discussions. 

Following the discussion in §VB| I assume that in both regions 2 and 3 the electrons 
have a power law distribution with a minimal Lorentz factor 7e im j n given by Eq. E21 with 
the corresponding Lorentz factors for the forward and the reverse shocks. 

Forward shock: The typical energy of synchrotron photons as well as the synchrotron 
cooling time depend on the Lorentz factor 7 e of the relativistic electrons under consideration 
and on the strength of the magnetic field. Using Eq. E21 for je.min and Eq. EH for the 
characteristic synchrotron energy for the forward shock: 

(hv syn ) obshe , min = 160 keV eTelT^n 1 / 2 = 0.5 keV (e B /0.1) 1/2 (e B /0.1) 2 gr Vl /2 , (69) 

and 

= °- 085 sec ^V^Vr 1 = 0-085 sec e^Vr'rToX 1 (70) 

The characteristic frequency and the corresponding cooling time for the "typical" electron 
are larger and shorter by a factor of [(p — 2)/{p— l)] 2 , correspondingly. 
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The electrons at the forward shock are fast cooling and w the typical cooling frequency 



373j: 



v c = 6A;eV r (e B /0.1)- 3/2 (r 2 /100)-V ' % 2 , (71) 

where t s is the time in seconds. The photons from the early forward shock are in the low 
7-rays to X-ray range, but this depends strongly on the various parameters (note the strong 
T2 dependence in equation IfiTJj). For this set of canonical parameters v m < v c . However, the 
ratio of these two frequencies depends on T 8 ! For T slightly larger then 100 the inequality 
will reverse and the system will be in the fast cooling regime. 

Reverse Shock: The Lorentz factor of the reverse shock, T 3 is smaller by a factor of 
£ 3 / 2 r than the Lorentz factor of the forward shock T 2 . Similarly the Lorentz factor of a 
"typical electron" in the reverse shock is lower by a factor £ 3//2 r. Therefore the observed 
energy is lower by a factor £ 3 r 2 . The typical synchrotron frequency of the reverse shock is 

^reverse shock = 1.3 X 10 13 Hz (e B /0.1) 1 / 2 (e e /0.1) 2 r 2 00 . (72) 

This is in the IR regions but note again the strong dependence on the Lorentz factor and on 
e e , which could easily bring this frequency up to the optical regime. The cooling frequency 
in the reverse shock region is th e same as the cooling frequency of the forward shock (if both 
regions have the same eg) jsT j | hence: 

^reverse shock = 8 X 10 18 Hz(e B /0. l)- 3 / 2 (r 2 /100)- V = 

8.8 x 10 15 Hz(e B /0.1)- 3/2 E 52 1/2 n- 1 t; 1 / 2 . (73) 

In the forward shock u m is comparable or larger than v c . In the reverse shock v m < v c 
and it is usually in the slow cooling regime. The reverse shocks exists for a short time 
until it reaches the back of the relativistic shell. Then it turns into a rarefraction wave that 
propagates forwards. After some back and forth bounces of these wave all the matter behind 



the forward shock organizes itself in the form of the Blan dford - McKee se l 



discussed latter in WH Al This above estimates suggest 



236, 
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similar solution 
373^ that during 



the short phase in which the reverse shock exists it should produce a po werfu l optical flash. 
This flash should coincide with the late part of the GRB. Kobayashi jl92| calculates the 
light curves and typical frequencies of the reverse shock for a variety of conditions. 
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D. The Transition from Internal Shocks to External Shocks 

The internal shocks take place at a distance Ri nt ~ c5tT 2 ~ (£t/0.3sec)r2l0 14 cm. These 
shocks last as long as the inner engine is active. The typical observed time scale for this 
activity ~ 50sec (for long bursts) and ~ 0.5sec (for short ones). External shocks begin at 
Rext ~ 10 16 cm. If R ext /T 2 < T = A/c, namely if the burst is long, the afterglows begins 
while internal sho cks a re still going on and the initial part of the afterglow overlaps the late 
part of the GRB (36^. At the early time the afterglow emission (from the forward shock) 
peaks in the high X-rays contributing also to the observed 7-ray flux. One can expect, 
therefore, a transition within the GRB from hard (pure GRB) to softer and smoother (GRB 
and afterglow) signal. Some observational evidence for this transition was presented in 



E. Prompt Polarization 



In flll A 41 1 discussed the detection of very high linear polarization from GRB 021206 
While the data analysis is uncertain several papers claimed that this detection has 
strong implications. First Coburn and Boggs suggest t hat t his polarizat ion i ndicates 
that the emission mechanism is synchrotron. Lyutikov et al. 2281 ] and Granot 1331 ] suggest 
further that it implies uniform magnetic fields within the emitting regions and the first even 
conclude that this implies that the relativistic flow is Poynting flux d omin ated and that the 



dissi pation is in the form of external plasma instability. Waxman 



4321 ] and Nakar et al. 



275] show however that: (i) Random magnetic field in shock's plane could produce almost 



as high polarization as a uniform field (provided that the emitting jet is narrow and one is 
looking along the edge of the jet), (ii) Even if the magnetic field is uniform the flow does not 
have to be Poynting flux dominated. They also stress that while in the uniform field case 
we expect high polarization in almost every burst in the random field one we can expect 
iigh polarization only in very few bursts. The different time dependence of the polarization 
27q could also e nable us to distinguish between the two possibilities. 



Lazzati et al. 



21 lj and Dar and De Rujula |7l| suggest that this polarization implies 



IC (which can have in principle higher intrinsic polarization). This shows that even the 
simplest conclusion (that the polarization confirms synchrotron as the emission mechanism) 
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is uncertain. My overall conclusion is that without further data on other bursts (which is, 
unfortunately, quite unlikely in the nearby future) not much can be learnt from this tentative 
detection. 



VII. THE AFTERGLOW 



It is generally accepted that the afterglow is produced when the relativistic ejecta is 
slowed down by the surrounding matter [236]. The afterglow begins at R ex t where most of 
the energy of the ejecta is transferred to the shocked exter nal medium. For a long burst 



363] and WIDJ) . Initially the 



this takes place while the burst is still going on (see Sari 
process might be radiative, namely a significant fraction of the kinetic energy is dissipated 
and the radiation process effects the hydrodynamics of the shock. I discuss this phase in 
§VII CI Later the radiation processes become less efficient and an adiabatic phase begins 
during which the radiation losses are minor and do not influence the hydrodynamics. The 
hydrodynamic evolution at this stage is adiabatic. If the ejecta is in the form of a jet with 
an opening angle 9 then a jet transition will take place when V reaches 9 . A transition 
into the Newtonian regime takes place when r — 1 « 0.5. I begin the discussion of the 
afterglow with the hydrodynamics of the adiabatic phase and with the resulting synchrotron 
light curve. I continue with a discussion of the possible early radiative evolution. Then I 
turn to the jet break and to the Newtonian transition. I continue with various complications 
and variations on these themes. 



A. Relativistic Blast Waves and the Blandford-McKee solution 

The theory of relativistic blast waves has been worked out in a classical paper by Bland- 
ford & McKee (BM) already in 1976. The BM model is a self-similar spherical solution 
describing an adiabatic ultra relativistic blast wave in the limit T ^> 1. This solution is the 
relativistic analogue of the well known Newtonian Sedov- Taylor solution. Blandford and 
McKee also describe in the same paper a generalization for varying ambient mass density, 
p = po(R/ Ro)~ k , R being the distance from the center. The latter case would be particularly 
relevant for k = 2, as expected in the case of wind from a progenitor, prior to the GRB 
explosion. 
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The BM solution describes a narrow shell of width ~ R/T 2 , in which the shocked material 
is concentrated. For simplicity I approximate the solution with a thin homogenous shell. 
Then the adiabatic energy conservation yields: 

E = ^(poR k )R 3 ~ k T 2 c 2 , (74) 

where E is the energy of the blast wave and Q is the solid angle of the afterglow. For a full 
sphere Q = Att, but it can be smaller if the expansion is conical with an opening angle 9: 
Q = 47r(l — cosO) ~ 2ti6 2 (assuming a double sided jet). This expression can be simplified 
using a generalized Sedov scale: 



I 



(3 - k)E/p R k c 2 V(3 " fc) . (75) 



If Q does not change with time then the blast wave collects ambient rest mass that equals 
its initial energy at R = /.If we take into account sideway expansion (after the jet break) we 
find that r « 1 and the blast wave becomes Newtonian at: 

R = 1(^1/ Art) 1 ' {z ~ k \ (76) 

Using the approximate (the numerical factor in this equation assumes that the shell is 
moving at a constant velocity) time - radius relation Eq. |H]one can invert Eq. [73] (using the 
definition of I , Eq. ITo"]) and obtain R and T as a function of time: 

R = [^!]V(4-*) t V(4-*) (77) 

/3— k 

The time in these expressions is the observer time - namely the time that photons emitted at 
R arrive to the observer (relative to the time that photons emitted at R = 0). For spherical 
(or spherical like) evolution Q in these expressions is a constant. In general it is possible 
that Q varies with R or with T (as is the case in a sideways expansion of a jet). This will 
produce, of course, a different dependence of R and T on t. 

The values of R and T from Eq. [78] can be plugged now into the typical frequencies z/ c , 
v m and v sa as well into the different expression for F vmax to obtain the light curve of the 
afterglow. 

Alternatively, one can calculate the light curve using a more detailed integration over the 
BM density and energy profiles. To perform such integration recall that the radius of the 

72 



front of the shock is: 

R = t{l-[2(4-k)T 2 ]- 1 }, (78) 

where T(t) is the shock's Lorentz factor and t is the time since the explosion in its rest frame. 
The different hydrodynamic parameters behind the shock can be expressed as functions of 
a dimensionless parameter x' 

x =[l + 2(i-k)r 2 ](l-R/t), (79) 

as: 



n 



2vW;T {10 - 3fc)/[2(4 - fe)] , 



7 2 = \T\- 

p = 2 ir 2 x -(17-4fc)/(12-3fc) > (8Q) 
o 

where n\ and W\ are the number density and enthalpy density of the undisturbed circumburst 
material and n and p are measured in the fluid's rest frame. 

The BM solution is self-similar and assumes r » 1. Obviously, it breaks down when 
R ~ I. This Relativistic-Newtonian transition should take place around 

t NR = l/c « 1.2 yr(£ is0i52 /m) 1/3 , (81) 

where the scaling is for k = 0, E 52 is the isotropic equivalent energy, Ei SO = 4tiE/Q, in units 
of 10 52 ergs and n\ is the external density in cm™ 3 . After this transition the solution will 
turn into the Newtonian Sedov- Taylor solution with: 

R = RN R {t/t NR ) 2/ \ (82) 
v = v NR (t/t NR )~ 3/5 , 

e = e NR {t/t NR )- & l b . (83) 

The adiabatic approximation is valid for most of the duration of the afterglow. However, 
during the first hour or so (or even for the first day, for k — 2), the system could be radiative 
(provided that e e ss 1) or partially radiative. During a radiative phase the evolution can be 
approximated as: 

E = —^—AR 3 ~ k TT c 2 , (84) 
3 — k 
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where To is the initial Lorentz factor. Cohen et al. |58j derived an analytic self-similar 
solution describing this phase. 

Cohen and Piran [stJ describe a solution for the case when energy is continuously added 
to the blast wave by the central engine, even during the afterglow phase. A self-similar 
solution arises if the additional energy deposition behaves l ike a power law. This would 
arise naturally in some models, e.g. in the pulsar like model j^uj. 



B. Light Curves for the "Standard" Adiabatic Synchrotron Model 

In WC2I I discussed the instantaneous synchrotron spectrum. The light curve that 
corresponds to this spectrum depends simply on the varia tion of the F v ^ max and the break 



236, 1376|| . This in turn depends on the 



frequencies clS 8b function of the observer time 
variation of the physical quantities along the shock front. For simplicity I approximate 
here the BM solution as a spherical homogeneous shell in which the physical conditions are 
determined by the shock jump between the shell and the surrounding matter. Like in §V C 21 
the calculatio n is divided to two cases: fast cooling and slow cooling. 

Sari et al. \s7(\ estimate the observed emission as a series of power law segments in time 
and in frequency 6 : 

F v oc r a v~P , (85) 

that are separated by break frequencies, across which the exponents of these power laws 
change: the cooling frequency, u c , the typical synchrotron frequency v m and the self absorp- 
tion frequency v sa . To estimate the rates one plugs the expressions for T and R as a function 
of the observer time (Eq. ITSj) . using for a homogenous external matter k = 0: 



R{t) = (17 Et/4irm p nc 



r(t) = (17E/10247rnm p c¥) 1/8 , (86) 

to the expressions of the cooling frequency, u c , the typical synchrotron frequency v m and the 
self absorption frequency v sa (Eqs. l2T)j) and to the expression of the maximal flux (Eq. EH1 



,1/4 



6 The following notation appeared in the astro-ph version of |376| . Later during the proofs that author 
reali zed that a is used often in astrophysics to denote a spectral index and in the Ap. J. version of 
|376| the notations have been changed to F u oc t~@v~ a . However, in the meantime the astro-ph notation 
became generally accepted. I use these notations here. 
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for slow cooling and Eq. |2H]for fast cooling). Note that the numerical factors in the above 
expressions arise from an exact integration over the BM profile. This procedure results in: 

u c = 0.85 x 10 14 Hz(e B /0.1)" 3/2 J E 52 1/2 nr 1 t d 1/2 , 
v m = 1.8 x 10 12 E Z (e B /0.iy/ 2 (e e /0.1) 2 EU 2 t- d 3/2 , 
F u , max = 0.35 x !0 5 ^(e B /OA) 1 / 2 E 52 n\ /2 D W 2 . (87) 

A nic e fea ture of this light curve is that the peak flux is constant and does not vary with 



time [259| | as it moves to lower and lower frequencies. 

At sufficiently early times v c < u m , i.e. fast cooling, while at late times v c > z/ m , i.e., 
slow cooling. The transition between the two occurs when v c = v m . This corresponds (for 
adiabatic evolution) to: 

to = 0.5 hours( ei j/0.1) 2 (e e /0.1) 2 E5 2 ni . (88) 

Additionally one can translate Eqs. EH to the time in which a given break frequency 
passes a given band. Consider a fixed frequency v = z/ 15 10 15 Hz. There are two critical 
times, t c and t m , when the break frequencies, v c and v mi cross the observed frequency v: 

t c = 0.2 hours(e B /0.1)" 3 J B 52 1 nf 2 z/ 1 - 5 2 , 
t m = 0.2 hours(e B /0.1) 1 / 3 (e e /0.1) 4 / 3 J E 5 1 2 /3 z/r 5 2/3 • (89) 

In the Ray leigh- Jeans part of the black body radiation I u = kT(2u 2 /c 2 ) so that 
F v oc kTv 2 . Therefore, in the part of the synchrotron spectrum that is optically thick to syn- 
chrotron self absorption, we have F v oc kT e ffV 2 . For slow cooling kT e ff ~ 7 m m e c 2 = const. 
throughout the whole shell of shocked fluid behind the shock, and therefore F v oc v 2 below 
u sa where the optical depth to synchrotron self absorption equals one, r Uag = 1. For fast 
cooling, as we go down in frequency, the optical depth to synchrotron self absorption first 
equals unity due to absorption over the whole shell of shocked fluid behind the shock, most 
of which is at the back of the shell and has kT e ff ~ 7 C . The observer is located in front of 
the shock, and the radiation that escapes and reaches the observer is from r v ~ 1. As v 
decreases below v sa the location where t v ~ 1 moves from the back of the shell toward the 
front of the shell, where the electrons suffered less cooling so that kT e ff(r u = 1) oc z/ -5 / 8 . 
Consequently F n u oc z/ 11 / 8 . At a certain frequency r v ~ 1 at the location behind the shock 
where electrons with 7 m start to cool significantly. Below this frequency, (zv), even though 
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t v ~ 1 closer and closer to the shock with decreasing u, the effective temperature at that 
location is constant: kT e ff ~ / -f m m e c 2 = const., and therefore F v oc v 2 for v < u ac , while 
F v oc z/ 11//8 for z/ ac < v < v sa . Overall the expression for the self absorption frequency de- 
pends on the cooling regime. It divides to two cases, d enot ed v sa and u ac , for fast cooling 
and both expression are different from the slow cooling jl43| |. For fast cooling: 

IV = 1.7 x 10 9 Hz (tB/O.iy^iejO.iy^E^nl^it/WOsec) 3 / 10 , (90) 

v sa = 1.8 x 10 10 Hz (eB/O.lp^V^^t/lOOsec)- 1 / 2 . (91) 

For slow cooling: 

v sa = 1.24 x 10 9 Hz ^'^^ l + zY%-hfntEl£ (92) 

For a given frequency either to > t m > t c (which is typical for high frequencies) or 
to < t m < t c (which is typical for low frequencies). The results are summarized in two tables 
HI and |H] describing a and (3 for fast and slow cooling. The different light curves are depicted 
in Fig. El 

These results are valid only for p > 2 (and for 7 max , the maximal electron energy, much 
higher than 7 m j n ). If p < 2 then ^ max plays a critical role. The resulting temporal and 
spectral indices for slow cooling with 1 < p < 2 are given by Dai and Cheng 64] and by 
Bhattacharya [3] and summarized in table IIHI below. For completeness I include in this 
table also the cases of propagation into a wind (see WH E|) and a jet break (see WH H|) . 

The simple solution, that is based on a homogeneous shell approximation, can be modified 



by using t he fu ll BM solution and integrating over the entire volume of shocked fluid 140]. 
Following j27l| I discuss £IV11 G II a simple way to perform this integration. The detailed 
integration yields a smoother spectrum and light curve near the break frequencies, but the 
asymptotic slopes, away from the break freque ncies and the transition times, remain the 



same as in the simpler theory. Granot and Sari |144j describe a detailed numerical analysis 
of the smooth afterglow spectr um in cluding a smooth approximation for the spectrum over 
the transition regions (see also jl48| ). They also describe additional cases of ordering of the 
typical frequencies which were not considered earlier. 

A final note on this "standard" model is that it assumes adiabaticity. However, in reality a 
fraction of the energy is lost and this influences over a long run the hydrodynamic behavior. 
This could be easily corrected by an integration of the energy losses and an addit ion a 
variable energy to Eq. |7U followed by the rest of the procedure described above 
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C. Light Curves for the early radiative phase 



If the electrons' energy is large (namely if e e is not far from unity), then early on during 
the first few hours of the afterglow there will be a radiative phase in which a significant 
fraction of the kinetic energy is lost via the radiative processes. One can generalize the 



BM solution to this radiative stage (see Cohen et al. 58j and WH A|) . The essence of 
the radiative phase is that in this case the energy varies as E oc T, where T = (R/L)~ 3 . 
Note that L is calculated in terms of M and the initial energy of the explosion, E$, via 
M = E /r c 2 , where T is the initial Lorentz factor of the ejecta: 

R(t) = (Act/L) 1/7 L, 
r(t) £ (Act/L)- 3/7 (93) 

The transition time from the radiative to the adiabatic phase takes place when the radiation 
losses become negligible. This happens at: 

t rad = 0.17 hours (e B /0.1) 7 / 5 (e e /0.1) 7 / 5 i; 5 4 2 /5 (r/100)- 4 / 5 n? /5 . (94) 



Following Sari et al. [376j one can use the above expressions to express the different 
typical frequencies and fluxes as: 

u c = 4.1 x 10 14 Hz (e B /0.1)- 3 / 2 E 5 - 2 4/7 (r/100) 4 / 7 nr 13/14 t d 2/7 , 
v m = 3.8 x 10 11 Hz (e B /0.1) 1 / 2 (e e /0.1) 2 J E; 4 2 /7 (r/100)- 4 / 7 nr 1/14 t d 12/7 , 
F^max = 1-4 x 10 3 /iJ efE 5 8 2 /7 (r/100)' 8 / 7 ^ /14 J D 28 2 t d 3/7 . (95) 

Like in the adiabatic case this can be translated to the times of passage of the break fre- 
quencies at a given observed frequency: 



t c = 0.05 x 10~ 7 days(e B /0.1)^ 21/4 E 52 2 r 2 n 1 



13/4 -7/2 



15 5 



t m = 0.01 days( eB /0.1) 7 / 24 (e e /0.1) 7 / 6 E 1 2 /3 r 2 1/ Vr 5 7/12 ^r 1/24 ■ (96) 

Unlike the adiabatic case, here v c must be below v m . Otherwise the bulk of the electrons do 
not cool and the system won't be radiative. Indeed at t ra d (given by Eq. above) v c = v m . 



D. Light Curve During the Newtonian transition 
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At t ~ t]\ffi (see Eq. ETJ the afterglow reaches the Newtonian Sedov- Taylor phas e. During 



this phase the adiabatic hydrodynamic is described by Eq. |H21 Frail et al. [107J calculate 
the synchrotron spectrum and light curve of the afterglow in this stage. The energy scaling 
implies that B oc t~ 3 ^ 5 and ^ e ,min oc t _6//5 . Combined together this yields v m oc t~ 3 . Using 
the standard assumptions of equipartition and of a power law electron's distribution they 
find : 

v m = lGHz(e B /0.3) 1/2 (e e /0.3)V 1/2 , 
F Vm<v<Vc = lmJ( eB /0.3) 3 / 4 (e e /0.3)r^ 4 £ 5 ^ . (97) 

This late time light curve provides a simple "calorimetric" estimate of the afterglow energy 
at this stage (see ^11 Dj) . Additionally as the radio flux is rather large and as it varies on 
a scal e of several month it can be used for search of orphan radio afterglows (see Levinson 
et al. Q and avi I K2jl . 

E. Generalizations: I. Winds 

The simplest generalization of the previous models is to allow a variable circuburst density 
with n(R) oc R~ k . The hydrodynamic evolution of a relativistic blast wave in such a medium 
was considered already in the original paper of Bla ndfor d and McKee jscj. The synchrotron 
light curve was considered first by Meszaros et al. \l62\ and by Dai and Lu 

Chevalier and Li j^, |^ stressed the importance of the n(R) oc R~ 2 case which arises 
whenever there is a stellar wind ejected by the GRB's progenitor prior to the burst. This 
arises naturally in the Collapsar model that is based on the collapse of a massive star. The 
calculations follow those outlines in the previous sections, with the only difference that the 
relations determining R(t) and T(t) for homogeneous circumburst medium, Eqs. IHHI should 
be replaced by Eqs. EElwith k = 2 

The high initial densities in a wind density profile implies a low initial cooling frequency. 
Unlike the constant density case the cooling frequency here increase with time 3|. This 
leads to a different temporal relations between the different frequencies and cooling regimes. 
For example it is possible that the cooling frequency will be initially below the synchrotron 
self absorption frequency. Chevalier and Li |50| consider five different evolution of the light 
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curves for different conditions and observed frequencies. We list below the two most relevant 
cases, the first fits the X-ray and optical afterglows while the second is typical for the lower 
radio frequencies. 

Note that for v m ,v c < v both the spectral slop and the temporal evolution are similar for a 
wind and for a constant density profile. This poses, of course, a problem in the interpretation 
of afterglow light curves. 



F. Generalizations: II. Energy injection and refreshed shocks 

The simple adiabatic model assumes that the energy of the GRB is constant. However, 
the energy could change if additional slower material is ejected behind the initial matter. 
This would be expected generically in the internal shock model. In this model the burst is 
produced by a series of collisions between shells moving at different velocities. One naturally 
expect here also slower moving matter that does not catch up initially with the faster moving 
one. However, as the initially faster moving matter is slowed down by th e circum-burs t 



207 



335 



matter this slower matter eventually catch up and produces refreshed shocks 

There are two implications for the refreshed shocks. F irst the additional energy injection 
will influence the dynamics of the blast wave 



335, 



366]. This effect can be modelled by 



modifying E in Eq. d| but the effect of additional mass carrying the slower energy must 
be included in some cases. This would change the decay slope from the canonical one 
and produce a slower decay in the light curve. In the following section WH Gl I describe 
a scheme for calculating the light curve resulting from a variable blast wave energy. If 
the additional matter is emitted sporadically then the shel l col lision could produce initial 
temporal variability in the early afterglow signal. Fox et al. , for example, suggest that 
refreshed shocks are the origin of the variability in the early afterglow of GRB 021004. 

A second effect is the production of a reve rse shock propagating into the slower material 
when it catches up with the faster one |207 |. This is of course in addition to the forward 
shock that propagates into the outer shell. This reverse shock could be episodal or long 
lasting depending on the profile of the additional matter. Kumar and Piran |207 | consider 
two shells with energies E\ and E2 in the outer and the inner shells respectively. The outer 
shell is moving with a bulk Lorentz factor Tq c ~ 5{t/day) 3 ^ 8 at the (observed) time, t, of the 
collision. As the inner shell catches up with the outer one when both shells have comparable 
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Lorentz factors the reverse shocks is always mildly relativistic. The calculation of the shock 
is slightly different than the calculation of a shell propagating into a cold material (another 
shell or the ISM) discussed earlier. Here the outer shell has already collided with the ISM. 
Hence it is hot with internal energy exceeding the rest mass energy. The reverse shock 
produces emission at a characteristic frequency that is typically much lower than the peak 
of the emission from the outer shell by a factor of ~ 7Tq c (E2/ Ei) 1 ' 1 , and the observed flux 
at this frequency from the reverse shock is larger compared to the flux from the outer shell 
by a factor of ~ %{Y Q rEJE^ ) 5/3 . This emission is typically m the radio or the FIR range. 

Kumar and Piran 2071 ] suggest that due to angular spreading the refre shed shocks produce 
an enhancement with a typical time scale 5t ~ t. Granot et al. jl38| stress that the fact 
that energy necessarily increases in refreshed shocks, the overall light curve must have a 
step-wise shape (above the continues power-law decline) with a break at the corresponding 
shocks. This behavior was s een in GRB 030329. However there the transitions are fast with 
5t < t. Granot et al. jlSsj ] point out that if the refreshed shocks take place after the jet 
break (as is likely the case in GRB 030329) then if the later shells remain cold and do not 
spread sideways we would have 5t ~ tj et < t. This explains nicely the fast transitions seen 
in this burst. 



G. Generalizations: III. Inhomogeneous density profiles 



An interesting possibility that arose with the observation of the variable light curve 
of the afterglow of GR B 021004 is th at the ejecta encounters surrounding matter with an 



irregular density profile 



165 



209 



274j . To explore this situation one can resor t to n umerical 



simulation of the propagation of the blast wave into a selected densi ty pr ofile |209j |. Instead 



one can attempt to model this analytically or almost analytically [2711 ] . The key for this 
analytic model is the approximation of the light curve from an inhomogeneous density profile 
as a series of emission from instantaneous BM solutions, each with its own external density. 
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1. The light curve of a BM solution 



The observed fl ux, a t an observer time t, from an arbitrary spherically symmetric emitting 



region is given by 



14oj: 



oo 



F„ (t ) = ^ jf it' l rHrf^ d(co S e) n ' (r)P 'f: K ' r) S(f - t - I£2=?), (98) 

where n' is the emitters density and P' u is the emitted spectral power per emitter, both 
are measured in the fluid frame; 9 is the angle relative to the line of sight, and A -1 = 



1/7(1 — v cos9/c) (v is t 



Nakar and Piran 



le emitting matter bulk velocity) is the blue-shift factor. 



27l\ show 7 that using the self-similar nature of the BM profile (with 



an external density oc r k ) one can reduce Eq. 1381 to: 

! rRmax(t) „ 

F v (t) = ^J Q A v (R)g p (t, k)dR . (99) 

The integration over R is over the shock front of the BM solution. The upper limit R m ax 
corresponds to the shock position from where photons leaving along the line of sight reach 
the observer at t. The factor D is the distance to the source (neglecting cosmological factors). 
(3 is the local spectral index. 

The factor g@ is a dimensionless factor that describes the observed pulse shape of an 
instantaneous emission from a BM profile. The instantaneous emission from a thin shell 
produces a finite pulse (see fllV Bl and Fig. HI)) . This is generalized now to a pulse from an 
instantaneous emission from a BM profile. Note that even though the BM profile extends 
from to R most of the emission arise from a narrow regions of width ~ R/T 2 behind the 
shock front. g@ is obtained by integration Eq. 19*81 over cosO and r, i.e. over the volume of the 
BM profile. It depends only on the radial and angular structure of the shell. The self-similar 
profile of the shell enables us to express g@ as a general function that depends only on the 
dimensionless parameter t = (t — ti os (R))/t ang (R), with ti os (R) is the time in which a photon 
emitted at R along the line of sight to the center reaches the observer and t ang = R/2cT 2 . 
The second function, A u , depends only on the conditions of the shock front along the line- 
of-sight. It includes only numerical parameters that remain after the integration over the 
volume of the shell. 

7 See |l4f)j for an alternative method for integrating Eq. [~jEl 
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When all the significant emission from the shell at radius R is within the same power-law 
segment, j3, (i.e v is far from the break frequencies) then A v and g@ are given by: 



A V (R) 



d2 4/3 pl/3 A/f 



1/3 



v>v c . 



erg 



sec ■ cm ■ Hz ' 



(100) 



where R is the radius of the shock front, n ext (R) is the external density, E is the energy in 
the blast-wave, M(R) the total collected mass up to radius R and H u is a numerical factor 
which depends on the observed power law segment (see |27l| for the numerical values. 



g(t,/3, k) 



l+2(4-fc)t - M (/3,fc) 



~ \ -(2-/3) 
2(4-fc)t+l \ ^ ' 



2(4-fc) 
l + F)-(2-/3) 



2(4-fc) X J 



where 



fi((3, k) = 3 • (71 - 17Jfe)/(72 - 18fc) - /3 • (37 + fe)/(24 - 6Jfe). 



(101) 



(102) 



This set of equations is completed with the relevant relations between the different variables 
of the blast wave, the observer time and the break frequencies. 

These equations describe the light curve within one pow er law segment of the light curve. 



Matching between different power laws can be easily done 



2711 ]. The overall formalism can 



be used to calculate the complete light curve of a BM blast wave. 



2. The light curve with a variable density or energy 

The results of the previous section can be applied to study the effect of variations in 
the external density or in the energy of the blast-wave by approximating the solution as a 
series of instantaneous BM solutions whose parameters are determined by the instantaneous 
external density and the energy. Both can vary with time. This would be valid, of course, 
if the variations are not too rapid. The light curve can be expressed as an integral over the 
emission from a series of instantaneous BM solutions. 

When a blast wave at radius R propagates into the circumburst medium, the emitting 
matter behind the shock is replenished within AR « R[2 l ^ 4 ~ k ^ — 1). This is the length scale 
over which an external density variation relaxes to the BM solution. This approximation 
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is valid as long as the density variations are on a larger length scales than AR. It fails 
when there is a sharp density increase over a range of AR. However, the contribution to the 
integral from the region on which the solution breaks is small (AR/R <C 1) and the overall 
light curve approximation is acceptable. Additionally the density variation must be mild 
enough so that it does not give rise to a strong reverse shock that destroys the BM profile. 

A sharp density decrease is more complicated. Here the length scale in which the emitting 
matter behind the shock is replenished could be of the order of R. As an example we consider 
a sharp drop at some radius Rj and a constant density for R> R^. In this case the external 
density is negligible at first, and the hot shell cools by adiabatic expansion. Later the forward 
shock becomes dominant again. Kumar and Panaitescu 204j show that immediately after 
the drop the light curve is dominated by the emission during the adiabatic cooling. Later the 
the observed flux is dominated by emission from R « R^, and at the end the new forward 
shock becomes dominant. Our approximation includes the emission before the density drop 
and the new forward shock after the drop, but it ignores the emission during the adiabatic 
cooling phase. 

As an example for this method Fig 1241 depicts the v m < v < v c light curve for a Gaussian 
(AR/R = 0.1) over-dense region in the ISM. Such a density profile may occur in a clumpy 
environment. The emission from a clump is similar to the emission from a spherically over- 
dense region as long as the clump's angular size is much larger than 1/T. Even a mild short 
length-scale, over-dense region (with a maximal over-density of 2) influences the light curve 
for a long duration (mainly due to the angular spreading). This duration depends strongly 
on the magnitude of the over-density. 

The calculations presented so far do not account, however, for the reverse shock resulting 
from density enhancement and its effect on the blast-wave. Thus the above models are 
limited to slowly varying and low contrast density profiles. Now, the observed flux depends 
on the external density, n, roughly as n 1 ^ 2 . Thus, a large contrast is needed to produce 
a significant re-brightening. Such a large contrast will, however, produce a strong reverse 
shock which will sharply decrease the Lorentz factor of the emitting matter behind the shock, 
T s h, causing a sharp drop in the emission below v c and a long delay in the arrival time of 
the emitted photons (the observer time is oc Tj 2 ). Both factors combine to suppresses the 
flux and to set a strong limit on the steepness of the re-brightening events caused by density 
variations. 



S3 



The method can be applied also to variations in the blast wave's energy. Spherically 
symmetric energy variations are most likely to occur due to refre s hed shoc ks, when new 



inner shells arrive from the source and refresh the blast wave 



207 



335 



this approximation misses the effect of the reverse shock that arise in this case 



366] . Once more, 



2071 ] . However 



it enables a simple calculation of the observed light curve for a given energy profile. 
H. Generalizations: IV. Jets 

The afterglow theory becomes much more complicated if the relativistic ejecta is not 
spherical. The commonly called "jets" corresponds to relativistic matter ejected into a cone 
of opening angle 9. I stress that unlike other astrophysical jets this ejecta is non steady 
state and generally its width (in the direction parallel to the motion) is orders of magnitude 
smaller than the radius where the jet is. A "flying pancake" is a better description for these 
jets. 

The simplest implication of a jet geometry, that exists regardless of the hydrodynamic 
evolution, is that once V ~ 9" 1 relativistic beaming of light will become less effective. The 
radiation was initially beamed locally into a cone with an opening angle T _1 remained inside 
the cone of the original jet. Now with > 9 the emission is radiated outside of the initial 
jet. This has two effects: (i) An "on axis" observer, one that sees the original jet, will detect 
a jet break due to the faster spreading of the emitted radiation, (ii) An "off axis" observer, 
that could not detect the original emission will be able to see now an "orphan afterglow", 
an afterglow without a preceding GRB (see WIIKj) . The time of this transition is always 
give by Eq. 11041 below with C-i = 1. 

lydr odynamic evolution of the source changes when T ~ 9~ x . Initially, 



Additionally the 

as long as T ^ 9~ 1 |305j ] the motion would be almost conical. There isn't enough time, in 
the blast wave's rest frame, for the matter to be affected by the non spherical geometry, and 
the blast wave will behave as if it was a part of a sphere. When T = C^" 1 , namely at 8 : 



The exact values of the uncertain constants Ci and C\ are extremely important as they determine the jet 
opening angle (and hence the total energy of the GRB) from the observed breaks, interpreted as tj et , in 
the afterglow light curves. 
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sideways propagation begins. The constant C\ expresses the uncertainty at relation between 
the Lorentz factor and the observing time and it depends on the history of the evolution 
of the fireball. The constant C% reflects the uncertainty in the value of T, when the jet 
break begins vs. the value of opening angle of the jet 9. For the important case of constant 
external density k = this transition takes place at: 

The sideways expansion continues with 9 ~ r -1 . Plugging this relations in Eq. 1751 and 
letting Q vary like T~ 2 one finds that: 

R « const ; (105) 
T « (R/2t) 1/2 . 



A more detailed analysis 204, 309, 343, 344] reveals that according to the simple one di- 
mensional analytic models T decreases exponentially with R on a very short l engt h scale. 9 

Table 1X1 describes the parameters a and (3 for a post jet break evolution 374j . The jet 
break usually takes place rather late, after the radiative transition. Therefore, I include in 
this table only the slow cooling parameters. 

An important feature of the post jet-break evolution is that v c , the cooling frequency 
becomes constant in time . This means that the high frequency (optical and X-ray ) optical 
spectrum does not vary after the jet-break took place. On the other hand the radio spectrum 
varies (see Fig. l23j) . giving an additional structure that confir ms the interpretation of break 
as arising due to a sideway s exp ansion of a jet (see e.g. jl5^|). 

Panaitescu and Kumar [290| find that the jet break transition in a wind profile will be 
very long (up to four decades in time) and thus it will be hard to observe a jet break in 
such a case. On the other hand it is interesting to note that for typical values of a seen 
after a jet break (a ~ —2) the high frequency spectral index, j3 = a/2 « — 1, is similar 
to the one inferred from a spherically symmetric wind (3 = (2a + l)/3 ~ —1 [157 |. Note 
however, that the wind interpretation requires a high (^ 3) p value (which may or may 



Note that the exponential behavior is obtained after converting Eq. to a differential equation and 
integrating over it. Different appr oxim ations used in deriving the differential equation lead to slightly 
different exponential behavior, see [309| . 
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not be reasonable). Still from the optical observations alone it is difficult to distinguish 
between these two interpretati ons. Here the radio observations play a crucial role as the 
radio behavior is very different \l0(\ . 

The sideways expansion causes a change in the hydrodynamic behavior and hence a break 
in the light curve. The beaming outside of the original jet opening angle also causes a break. 
If the sideways expa nsion is at the speed of light than both transitions would take place at 
the same time |374j . If the sideways expansion is at the sound speed then the beaming 
Tansition would take place first and only later the hydrodynamic transition would occur 
292^ . This would cause a slower and wider transition with two distinct breaks, first a steep 
break when the edge of the jet becomes visible and later a shallower break when sideways 
expansion becomes important. 

T he analytic or s emi- analytic calculations of synchrotron radiation from jetted afterglows 



204, 



264, 



292 



3441 . 374j have led to different estimates of the jet break time tj et and of the 

curves assuming emission from 
one representative point, and obtained a smooth 'jet break', extending ~ 3 — 4 decades in 



time, after which F v>Vm oc t p . Sari et al. 374] assume that the sideways expansion is at the 
speed of light, and not at the speed of s ound (c/y3) as others assume, and find a smaller 



value for ij et . Panaitescu and Meszaros 



292] included the effects of geometrical curvature 



and finite width of the emitting shell, along with electron cooling, and obtained a relatively 
sharp break, extending ~ 1 — 2 decades in time, in the optical light curve. Moderski et al. 
[264] used a slightly different dynamical model, and a different formalism for the evolution of 
the electron distribution, and obtained that the change in the temporal index a (F v oc t~ a ) 
across the break is smaller than in analytic estimates (a = 2 after the break for v > v m , 
p = 2.4), while the break extends over two decades in time. 

The different analytic or semi-analytic models have different predictions for the sharpness 
of the 'jet break', the change in the temporal decay index a across t he br eak and its asymp- 
totic value after the break, or even the very existence a 'jet break' |l7l| . All these models 
rely on some common basic assumptions, which have a significant effect on the dynamics 
of the jet: (i) the shocked matter is homogeneous (ii) the shock front is spherical (within a 
finite opening angle) even at t > tj Ct (iii) the velocity vector is almost radial even after the 
jet break. 



However, recent 2D hydrodynamic simulations 



1371 ] show that these assumptions are not 
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a good approximation of a realistic jet. Using a very different approach Cannizzo et al. 
48j find in another set of numerical simulations a similar result - the jet does not spread 
sidew ays a s much. Figure shows the jet at the last time step of the simulation of Granot 
et al. |l37 |. The matter at the sides of the jet is propagating sideways (rather than in the 
radial direction) and is slower and much less luminous compared to the front of the jet. 
The shock front is egg-shaped, and quite far from being spherical. Figure EH shows the 
radius R, Lorentz factor T, and opening angle 9 of the jet, as a function of the lab frame 



time. The rate of increase o 
predicted by simple models 



9 w ith R ex ctiah, is much lower than the exponential behavior 



204, 



309 



343 



)44j |. The value of 9 averaged over the emissivity 



is practically constant, and most of the radiation is emitted within the initial opening angle 
of the jet. The radius R weighed over the emissivity is very close to the maximal value of R 
within the jet, indicating that most of the emission originates at the front of the jet 10 , where 
the radius is largest, while R averaged over the density is significantly lower, indicating that 
a large fraction of the shocked matter resides at the sides of the jet, where the radius is 
smaller. The Lorentz factor V averaged over the emissivity is close to its maximal value, 
(again since most of the emission occurs near the jet axis where V is the largest) while V 
averaged over the density is significantly lower, since the matter at the sides of the jet has a 
much lower V than at the front of the jet. The large differences between the assumptions of 
simple dynamical models of a jet and the results of 2D simulations, suggest that great care 
should be taken when using these models for predicting the light curves of jetted afterglows. 
Since the light curves depend strongly on the hydrodynamics of the jet, it is very important 
to use a realistic hyd ro dynamic model when calculating the light curves. 

Granot et al. 1371 ] used 2D numerical simulations of a jet running into a constant density 
medium to calculate the resulting light curves, taking into account the emission from the 
volume of the shocked fluid with the appropriate time delay in the arrival of photons to 
different observers. They obtained an achromatic jet break for v > f m (tj e t) (which typically 
includes the optical and near IR), while at lower frequencies (which typically include the 
radio) there is a more moderate and gradual increase in the temporal index a at tj et , and 
a much more prominent steepening in the light curve at a latter time when v m sweeps past 



10 This may imply that the expected rate of orphan afterglows should be smaller than estimated assuming 
significant sideways expansion! 
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the observed frequency. The jet break appears sharper and occurs at a slightly earlier time 
for an observer along the jet axis, compared to an observer off the jet axis (but within the 
initial opening angle of the jet). The value of a after the jet break, for v > u m , is found to 
be slightly larger than p (a = 2.85 for p = 2.5). Due to the fact that a significant fraction 
of the jet break occurs due to the relativistic beaming effect (that do not depend on the 
hydrodynamics) in spite of the different hydrodynamic behavior the numerical simulations 
show a jet break at roughly the same time as the analytic estimates. 



I. Generalizations: V. Angular Dependent Jets and the Structured Jet Model 



In a realistic jet one can expect either a random or regular angular dependent structure. 
Here there are two dominant effects. As the ejecta slows down its Lorentz factor decreases 
and an observer will detect radiation from an angular region of size T _1 (see §IVC)) . At 
the same time the mixing within the ejecta will lead to an intrinsic averaging of the angular 
structure. Thus, b oth effects lead to an averaging over the angular structure at later times. 



Several authors 



219 



347 



44q suggested independently a different interpretation to the 



observed achromatic breaks in the afterglow light curves. This interpretation is based on a 
jet with a regular angular structure. According to this model all GRBs are produced by a 
jet wi th a fixed angular structure and the break corresponds to the viewing angle. Lipunov 



et al. 



219j | considered a "u niversal" jet with a 3-step pro file: a spherical one, a 20° one and 



a 3° one. Rossi et al. |347j and Zhang and Meszaros 446] considered a special profile where 
the energy per solid angle, e(9), and the Lorentz factor T(t = 0, 6) are: 



< 9 < 6 r 



<0<0j 



(106) 



and 



(107) 



r c o < e < e c 

T c (l)- b , b>0 6 c <6<6j, 
where 9j is a maximal angle and core angle, 8 C , is introduced to avoid a divergence at 9 = 
and the parameters a and b here define the energy and Lorentz factor angular dependence. 
This core angle can be taken to be smaller than any other angle of interest. The power law 
index of T, b, is not important for the dynamics of the fireball and the computation of the 
light curve as long as T{t = 0,6) = T (6) > and F (6) > 1. 



angular structure with a = 2. Rossi et al. 347] approximate the evolution assuming that at 
every angle the matter behaves as if it is a part of a regular BM profile (with the local e and 
F(t, 9)) until r(i, 9) = 9~ l . Then the matter begins to expand sideways. The resulting light 
curve is calculated by averaging the detected light resulting from the different angles. They 
find that an observer at an angle o will detect a break in the light curve that appears around 
the time that T(t,9 ) = 9^ 1 (see Fig. l2*Hj) . A simple explanation of break is the following: 
As the evolution proceeds and the Lorentz factor decreases an observer will detect emission 
from a larger and larger angular regions. Initially the additional higher energy at small 
angles, 9 < 9 Q compensates over the lower energies at larger angles 9 > 9 Q . Hence the 
observer detects a roughly constant energy per solid angle and the resulting light curve is 
comparable to the regular pre-jet break light curve. This goes on until r _1 (0) = 0( o ). After 
this stage an further increase in the viewing angle will result in a decrease of the energy 
per unit solid angle within the viewing cone and will lead to a break in the light curve. 

This interpretation of the breaks in the light curves in terms of the viewing angles of a 
standard structured jets implies a different understanding of total energy within GRB jets 
and of the rate of GRBs. The total energy in this model is also a constant but now it is larger 
as it is the integral of Eq. I106l over all viewing angles. The distribution of GRB luminosities, 
which is interpreted in the uniform jet interpretation as a distribution of jet opening angles 
is interpreted here as a distribution of viewing angles. As such this distribution is fixed by 
geometrical reasoning with P(9 )d9 oc sin9 d9 (up to the maximal observing angle 9j). 
This leads to an implied isotropic energy distribution of 

P{log{E iso )) oc E£ . (108) 

Guetta et al. and Nakar et al. [266] find that these two distributions are somewhat 

inconsistent with current observations. However the present data that suffers from numerous 
observational biases is insufficient to reach a definite conclusions. 

In order to estimate better the role of the hydrodynamics on the light curves of a struc- 
tnred jet Granot and Knn lar Q, Granot e, a,,. Q considered two snnp.e .node.s tot the 
hydrodynamics. In the first (model 1) there is no mixing among matter moving at different 
angles i.e. e(9, t) = e(9, to)- In the second (model 2) e is a function of time and it is averaged 
over the region to which a sound wave can propagate (this simulates the maximal lateral 



To fit the constant energy result |105l |291. 310] Rossi et al. [347] consider a specific 
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energy transfer that is consistent with causality) . They consider various energy and Lorentz 
factors profiles and cal culat e the resulting light curves (see Fig. 



Granot and Kumar |135| | find that the light curves of models 1 and 2 are rather similar in 
spite of the different physical assumptions. This suggests that the widening of the viewing 
angle has a more dominant effect than the physical averaging. For models with a constant 
energy and a variable Lorentz factor ((a, b) = (0,2)) the light curve initially rises and 
there is no jet break, which is quite different from observations for most afterglows. For 
(a, b) = (2,2), (2,0) they find a jet break at tj when Y{6 ) ~ 9~ l . For (a, b) = (2,2) the 
value, ai, of the temporal decay slope at t < tj increases with 9 a , while a 2 = oc{t > tj) 
decreases with o . This effect is more prominent in model 1, and appears to a lesser extent 
in model 2. This suggests that 5a = a.\ — a 2 should increase with tj, which is not supported 
by observations. For (a, b) = (2, 0), ther e is a flattening of the light curve just before the jet 
break (also noticed by Rossi et al. |347 |). for 6 Q > 39 c . Again, this effect is larger in model 
1, compared to model 2 and again this flattening is not seen in the observed data. 

Clearly a full sol utio n of an angular dependent jet requires full numerical simulations. 
Kumar and Granot I202I present a simple 1-D model for the hydrodynamics that is obtained 
by assuming axial symmetry and integrating over the radial profile of the flow, thus con- 
siderably reducing the computation time. The light curves that they find resemble those of 
models 1 and 2 above indicating that these crude approximations are useful. Furthermore 
they find relatively little change in e(6) within the first few days, suggesting that model 
1 is an especially useful approximation for the jet dynamics at early times, while model 2 
provides a better approximation at late times. 

J. Afterglow Polarization - a tool that distinguished between the different jet 
models 

Synchrotron emission from a jet (in which the spherical symmetry is broken) would 



naturally produce polarized emission 



126, 



147. 13651 . Moreover, the level and the direction 



of the polarization are expected to vary with time and to give observational clues on the 
geometrical structure of the emitting jet and our observing angle with respect to it. 

The key feature in the determination of the polarization during the afterglow is the 
varying Lorentz factor and (after jet break) varying jet width. This changes changes the 
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overall geometry (see Fig. ITo^ and hence the observer sees different geometries |l77t \3Ql 
Initially, the relativistic beaming angle 1/T is narrower than the physical size of the jet 60, 
and the observer see a full ring and therefore the radial polarization averages out (the first 
frame, with T8o = 4 of the left plot in Fig. I3UJ) . As the flow decelerates, the relativistic 
beaming angle 1/T becomes comparable to #0 an d only a fraction of the ring is visible; net 
polarization is then observed. Assuming, for simplicity, that the magnetic field is along 
the shock then the synchrotron polarization will be radially outwards. Due to the radial 
direction of the polarization from each fluid element, the total polarization is maximal when 
a quarter (T6q = 2 in Figure IHU|) or when three quarters (T8q = 1 in Figure IHUJ) of the ring are 
missing (or radiate less efficiently) and vanishes for a full and a half ring. The polarization, 
when more than half of the ring is missing, is perpendicular to the polarization direction 
when less than half of it is missing. 

At late stages the jet expands sideways and since the offset of the observer from the 
physical center of the jet is constant, spherical symmetry is regained. The vanishing and 
re-occurrence of significant parts of the ring results in a unique prediction: there should 
be three peaks of polarization, with the polarization position angle during the central peak 
rotated by 90° with respect to the other two peaks. In case the observer is very close to the 
center, more than half of the ring is always observed, and therefore only a single direction of 
polarization is expected. A few possible polarization light curves are presented in Fig. EOJ 

The predicted polarization from a structured jet is drastically diffe rent from the one from a 
uniform jet, providing an excellent test between the two models 1.348^ . Within the structured 
jet model the polarization arises due to the gradient in the emissivity. This gradient has 
a clear orientation. The emissivity is maximal at the center of center of the jet and is 
decreases monotonously outwards. The polarization will be maximal when the variation in 
the emissivity within the emitting beam are maximal. This happens around the jet break 
when 9 Q bs ~ and the obser ved beam just reaches the center. The polarization expected 



in this case is around 20% 



348] and it is slightly larger than the maximal polarization from a 



uniform jet. As the direction of the gradient is always the same (relative to a given observer) 
there should be no jumps in the directi on o f the polarization. 



According to the patchy shell model 



2061 ] the jet can includes variable emitting hot spots. 



This could lead to a fluctuation in the light curve (as h ot spots e nter the observed beam) and 



also to corresponding fluctuations in the polarization 



133, 



There is a clear prediction 
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267 



2741 ] that if the fluctuations are angular fluctuations have a typical angular scale 9j 



then the first bump in the light curve should take place on the time when ~ 9f (the 
whole hot spot will be within the observed beam). The following bumps in the light curve 
should decrease in amplitude (due to statistical fluctuations). Nakar and Oren (267 1 show 
analytically and numerically that the jumps in the polarization direction should be random, 
sharp and accompanied by jumps in the amount of polarization. 



K. Orphan Afterglows 

Orphan afterglows arise as a natural prediction of GRB jets. The realization that GRBs 
are collimated with rather narrow opening angles, while the following afterglow could be 
observed over a wider angular range, led immediately to the search for orphan afterglows: 
afterglows which are not associated with observed prompt GRB emission. While the GRB 
and the early afterglow are collimated to within the original opening angle, 9j, the afterglow 
can be observed, after the jet break, from a viewing angle of The Lorentz factor, T, is 
a rapidly decreasing function of time. This means that an observer at 9 Q b s > 9j couldn't see 
the burst but could detect an afterglow once = 9 t, s . As the typical emission frequency 
and the flux decrease with time, while the jet opening angle 9 increases, this implies that 
observers at larger viewing angles will detect weaker and softer afterglows. X-ray orphan 
afterglows can be observed several hours or at most a few days after the burst (depending of 
course on the sensitivity of the detector). Optical afterglows (brighter than 25th mag ) can 
be detected in R band for a week from small (~ 10°) angles away from the GRB jet axis. 
On the other hand, at very late times, after the Newtonian break, radio afterglows could be 
detected by observers at all viewing angles. 

The search for orphan afterglows is an observational challenge. One has to search for a 
10~ 12 ergs/sec/cm 2 signal in the X-ray , a 23th or higher magnitude in the optical or a mJy 
in radio (at GHz) transients. Unlike afterglow searches that are triggered by a well located 
GRB for the orphan afterglow itself there is no information where to search and confusion 
with other transients is rather easy. So far there was no detection of any orphan afterglow 
in any wav eleng th. 



Rhoads 



343] was the first to suggest that observations of orphan afterglows would enable 



us to estimate the opening angles and the true rate of GRBs. Dalai et al. |£9j have pointed 
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out that as the post jet-break afterglow light curves decay quickly, most orphan afterglows 
will be dim and hence undetectable. They point out that if the maximal observing angle, 
dmaxi of an orphan afterglow will be a constant factor times 9j the ratio of observed orphan 
afterglows, R^t^-, to that of GRBs, Rqrbi w ^ n °t tell us much about the opening angles 
of GRBs and the true rate of GRBs, -R^rs = hR-GRB- However as we see below this 
assumption is inconsistent with the constant energy of GRBs that suggests that all GRBs 
will be detected to up to a fixed angle which is independent of their jet opening angle. 



1. Optical Orphan Afterglow 



Optical orphan afterglow is emitted at a stage when the ou tflow is st i 



105, 291 



1 re lativistic. The 
3101 and that the 



observation that GRBs have a roughly constant total energy 
observed variability in the apparent luminosity arises mostly from variation in the jet opening 
angl es leads to a remarkable result: The post jet-break afterglow light curve is universal 
jl39j|. Fig. depicts this universal light curve. This implies that for a given redshift, z, 
and a given limiting magnitude, m, there will be a fixed 9 max (z,m) (independent of 9j, for 
9j < 9 max ) from within which orphan afterglow can be detected 



This universal post jet-break li ght c urve can be estimated from the observations |409j or 



alternatively from first principles 273] . An observer at 9 obs > 9j will (practically) observe 
the afterglow emission only at tg when T = 9~ bs . Using Eq. 11041 and the fact that T oc t~ 1//2 
after the jet break (Eq. HOfij) one can estimate the time, tg when a emission from a jet would 
be detected at 9 obs : 

tg = A(9 0bs /9 j ) 2 t jet , (109) 

where A is a factor of order unity, and tj et is the time of the jet break (given by Eq. Il()4j) . 
The flux at this time is estim ated by substitution of this value into the post-jet-break light 
curve (see Nakar et al. 273] for details): 

F(9 obs )=F f(z)92 P , (HO) 

where F is a constant and f(z) = (1 + z) 1+/3 .D^f 8 includes all the cosmological effects and 
Dl28 is the luminosity distance in units of 10 28 cm. One notices here a very strong dependence 
on 9 obs- The peak flux drops quickly when the observer moves away from the axis. Note also 
that this maximal flux is independent of the opening angle of the jet, 9j. The observations 
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of the afterglows with a clear jet break (GRB990510 [159j,|39J, and GRB000926 |l60j ) can 
be used to calibrate F . 

Now, using Eq. I11UI one can estimate max (z, m) and more generally the time, 
(t b s (z,6,m) that a burst at a redshift, z, can be seen from an angle 9 above a limiting 
magnitude, m: 

t obs (z,e,m) « ^(e 2 max - e 2 obs ) . (in) 

One can then proceed and integrate over the cosmological distribution of bursts (assuming 
that this follows the star formation rate) and obtain an estimate of the number of orphan 
afterglows that would appear in a single snapshot of a given survey with a limiting sensitivity 

Flim- 

f°° n(z) dV(z) rdmax(z,m) 

N orph = / n \\ , { U z x / t obs {z, 9, m)6de oc (F F lim f /p • 112 

Jo (1 + z) dz JOj 

where n(z) is the rate of GRBs per unit volume and unit proper time and dV(z) is the 

differential volume element at redshift z. Note that m odifi c atio ns of this simple model may 
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arise with more refined models of the jet propagation 

The results of the intergration of Eq. 11121 are depicted in Fig. EH Clearly the rate 
of a single detection with a given limiting magnitude increases with a larger magnitude. 
However, one should ask what will be the optimal strategy for a given observational facility: 
short and shallow exposures that cover a larger solid angle or long and deep ones over a 
smaller area. The exposure time that is required in order to reach a given limiting flux, Fu m , 
is proportional to F^. Dividing the number density of observed orphan afterglows (shown 
in Fig. I3*2"j) by this time factor results in the rate per square degree per hour of observational 
facility. This rate increases for a shallow surveys that cover a large area. This result can be 
understood as follows. Multiplying Eq. 11121 bv Ff im shows that the rate per square degree 
per hour of observational facility oc Ff~^^ p . For p > 1 the exponent is positive and a shallow 
survey is preferred. The limiting magnitude should not be, however, lower than ~ 23rd as 
in this case more transients from on-axis GRBs will be discovered than orphan afterglows. 

Using these estimates Nakar et al. |273| find that with their most optimistic parameters 
15 orphan afterglows will be recorded in the Sloan Digital Sky Survey (SDSS) (that covers 
10 4 square degrees at 23rd mag) and 35 transients will be recorded in a dedicated 2m class 
;eles cope operating full time for a year in an orphan afterglow search. Totani and Panaitescu 
40s| find a somewhat higher rate (a factor ~ 10 above the optimistic rate). About 15% of 
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the transients could be discovered with a second exposure of the same area provided that 
it follows after 3, 4 and 8 days for mu m = 23, 25 and 27. This estimate does not tackle 
he challenging problem of identifying the afterglows within the collected data. Rhoads 
34s| suggested to identify afterglow candidates by comparing the multi-color SDSS data 
to an after glow template. One orphan afterglow candidate was indeed ide ntifie d using this 
technique \i2(\ . However, it turned out that it has been a variable AGN This event 

demonstrates the remarkable observational challenge involved in this project. 



2. Radio Orphan Afterglow 

After the Newtonian transition the afterglow is expanding spherical. The velocities are at 
most mildly relativistic so there are no relativistic beaming effects and the afterglow will be 
observed from all viewing angles. This implies that observations of the rate of orphan GRB 
afterglows at this stage will give a direct measure of the beaming factor of GRBs. Uppe r 
limits on the rate of orphan afterglows will provide a limit on the beaming of GRBs 298]. 
However, as I discuss shortly, somewhat surprisingly, upper limits on the rate of orphan 
radio afterglow (no det ectio n of orphan radio afterglow) provide a lower (and not upper) 



limit on GRB beam ing 



215 1 . 



Frail et al. |l07| estimate the radio emission at this stage using the Sedov- Taylor solution 
for the hydrodynamics (see §VII D|l . They find that the radio emission at GHz will be around 
1 mJy at the time of the Newtonian transition (typically three month after the burst) and 
it will decrease like t~ 3 (p _1 )/ 2 + 3 / 5 (see Eq. l9"Tj) . Using this limit one can estimate the rate 
of observed orphan radio afterglow within a given limiting flux. The beaming factor f b ~ l 
arises in two places in this calculations. First, the overall rate of GRBs: RcriB = fbR°GRBi 
increases with f b . Second the total energy is proportional to f^ 1 hence the flux will decrease 
when Jb increases. The first factor implies that the rate of orphan radi o af terglows will 



increase like f b . To estimate the effect of the second factor Levinson et al. |215| use the fact 
that (for a fixed observed energy) the time that a radio afterglow is above a given flux is 
proportional to E 10 / 9 in units of the NR transition time which itself is proportional to E 1 ^ 3 . 
Overall this is propo rtional to E 13 ^ 9 and hence to f b 13 ^ 9 . To obtain the overall effect of f b 
Levinson et al. [215] integrate over the redshift distribution and obtain the total number of 
orphan radio afterglow as a function of f b . For a simple limit of a shallow survey (which is 
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applicable to current surveys) typical distances are rather "small", i.e. less than lGpc and 
cosmological corrections can be neglected. In this case it is straight forwards to carry the 
integration ana lytic ally and obtain the number of radio orphan afterglows in the sky at any 



given moment 



215|: 



N R * 10VfW0.5) (|^)" /2 (^|) 3/2 (^) 9/ ^T 4 < 5 V 9 - 3/4 (V3t^)-/ 2 o. 

m V ' ' (113) 

where R is the observed rate of GRBs per Gpc 3 per year, and £j is the time in which the 
radio afterglow be come s isotropic. 



Levinson et al. 



2151 ] search the FIRST and NVSS surveys for point-like radio transients 



with flux densities greater than 6 mJy They find 9 orphan candidates. However, they argue 
that the possibility that most of these candidates are radio loud AGNs cannot be ruled out 
without further observations. This analysis sets an upper limit for the all sky number of 
radio orphans, which corresponds to a lower limit f^ 1 > 10 o n the beaming factor. Rejection 
of all candidates found in this search would imply > 40 



L. Generalizations: VI. Additional Physical Processes 

With the development of the theory of GRB afterglow it was realized that several ad- 
ditional physical ingredients may influence the observed afterglow light emission. In this 
section I will review two such processes: (i) Pre acceleration of teh surrounding matter by 
the prompt 7-rays emission and (ii) Decay of neutrons within the outflow. 



1. Pre- acceleration 



The surrounding regular ISM or even stellar wind is optically thin to the initial 7- 
rays pulse. Still the inte racti on of the pulse and the surrounding matter may not be trivial. 



Thompson and Madau 



4071 ] pointed out that a small fraction of the 7-rays radiation will 



be Compton scattered on the surrounding electrons. The backscattered photons could now 
interact with the outwards going 7-rays flux and produce pairs. The pairs will increase the 
rate of backscattering and this could lead to an instability. When sufficient number of pairs 
will be produced the surrounding matter will feel a significant drag by the 7-rays flux and 
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it will be accelerated outwards 2441. These pre-acceleration of the ambient medium could 

nn 

have several implications to the early afterglow [2JJ, [233J . 

The key issue is that while the optical depth of the surrounding medium (as "seen" 
the 7-rays photons) is very small, the mean free path of an ambient electron within the 
7-rays photons is large (at small enough radius) and each electron scatters many photons. 
While the medium absorbs only a small fraction of the prompt 7-rays energy, the effect 
of this energy can be significant. Beloborodov j^j characterizes the interaction of the 7- 
rays radiation front with the surrounding medium by a dimensionless parameter 11 : 

V = , Zf lS ° 2 = 6.5£ 52 i^ 2 , (H4) 



the energy that a single electron scatters relative to its rest mass. Beloborodov [2l| calculates 
the Lorentz factor of the ambient medium and the number of pairs per initial electron as 
functions of 77. where r]i oaa < = 20 — 30, depending on the spectrum of the gamma-rays, 
Vacc = Svioad = 100 - 150, and f acc = [exp(r] acc /r]i oad ) + exp(-r] acc /r]i oad )]/2 = 74. 

If 77 < rji oa d m 20 — 30, depending on the spectrum of the gamma-rays, the medium 
remains static and e ± -free. When the front has rj > r)i oad , a runaway loading occurs. The 
number of loaded pairs depends exponentially on 77 as long as 77 < r\ acc = 5r]i oad = 100 — 150. 
The medium is accelerated if 77 > r\ acc . r\ acc is around 100 because the electrons are coupled 
to the ambient ions, and and the other hand the loaded e ± increase the number of scatters 
per ion. At 77 = r\ gav m 3 x 10 3 , the matter is accelerated to a Lorentz factor T am b ie nt 
that exceeds the Lorentz factor of the ejecta. It implies that the radiation front pushes the 
medium away from the ejecta and opens a gap. 

As the GRB radiation front expands, the energy flux and hence 77 decreases oc R~ 2 . 
77 passes through 77^, r] acc , and r]i oad at R gap , R acc , and R toad , respectively. These three 
characteristic radii define four stages: 

i. R < Rga P ~ -Racc/3: The ejecta moves in a cavity produced by the radiation front with 

^ambient ^ejecta- 

II. R gap < R < R acc ~3x 10 15 cmE^ 2 cm: The ejecta sweeps the e ± -rich medium that 
has been preaccelerated to 1 <C F ambient < T e j ecta . 

III. R acc < R < Rioad ~ 2.3-R acc . The ejecta sweeps the "static" medium {T a mbient ~ 1) 
which is still dominated by loaded e ± . 



11 Note that Beloborodov uses the notation £ for this parameter 
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IV. R > Rioad- The ejecta sweeps the static pair-free medium. 

This influence of the 7-rays on the surrounding matter may modify the standard picture 
of interaction of external shocks with the surrounding medium (see £IV1 C II This depends 
mostly on the relation between R ext and R gap ~ lO 15 -^^ 2 cm. If R ext > R gap this effect 
won't be important. However, if R ext < R gap then effective decceleration will begin only at 
Rga P - At R < Rg a p the ejecta freely moves in a cavity cleared by the radiation front and 
only at R = R gap the blast wave gently begins to sweep the preaccelerated medium with 
a small relative Lorentz factor. With increasing R > R gap , ^ ambient falls off quickly, and it 
approaches T am bient = 1 at R = R acc ~ 3R gap as Y ambient = {R/R aC c)~ 6 - Thus, after a delay, 
the ejecta suddenly "learns" that there is a substantial amount of ambient material on its 
way. This resembles a collision with a wall and results in a sharp pulse (see Fig. . 

While Rgap does not depend on the external density R ext does (see Eq. ESJ). The condition 
R ext < Rg a p implies: 

Ell 6 n\ /3 T%1 > 0.02 (115) 

Thus it requires a dense external medium and large initial Lorentz factor. Otherwise R gap 
is too large and the deceleration takes place after the gap is closed. Hence the conditions 



for pre-acceleration will generally occur i 



the burst takes place in a d ense circumburst 



205] elaborate on this 



regions, like in a Wolf-Rayet progenitor |21|. Kumar and Panaitescu 
model and find that the observational limits by LOTIS and ROTSE on prompt emission 
from various burst limit the ambient ISM density (within 10 16 cm to less than I0 3 rmcm~ 3 . 
Similarly the find that in case of a wind the progenitors mass loss to wind's velocity ratio is 
below 10- 6 M o /yr/(10 3 km/sec). 



2. Neutron decoupling and decay 

Derishev et al. [3, Illl pointed out that neutrons that are included initially in the fireball 
will change its dynamics and modify the standard afterglow evolution. While the protons 
slow down due to the interaction with the surrounding matter the neutrons will coast freely 
after they decouple with T n , which equals to the Lorentz factor while decoupling took place. 

At 

Rdecay ~ 0.3 X 10 16 C m(r„/100) (116) 
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the neutrons decay. A new baryonic shell forms ahead of the original fireball shell, with 
energy comparable to the initial energy of the protons' shell (this depends, of course, on 
the initial ratio of neutrons to protons). At this stage the neutrons front that is not slowed 
down like the rest of the fireball is at a distance: 

AR = R(1/2T 2 - l/2Gamma 2 n ) , (117) 

from the fireball front, where T is the current Lorentz factor of the fireball. 

Once more the situation depends on whether Rdecay is smaller or larger than R ex t, the 
original deceleration radius. If Rdecay < Rext- 

E^n^T^iTn/Zm)- 1 < 0.06 , (118) 

the decaying neutron products will mix with the original protons and won't influence the 
evolution significantly (apart from adding their energy to the adiabatic fireball energy). 
Otherwise, they sit uatio n depends on Y n the Lorentz factor at decoupling. 

Pruet and Dalai |325^ consider a situation in which the neutron decouple with a low T n . 
In this case one will get a delayed shock scenario when the neutronic decay produce will 
eventually catch up with the slowing down protons (when their Lorentz factor is of order T n . 
Along the same line of thought Dalai et al. j3] suggest that a large neutronic component 
that ma y exi st within the initial fireball material may help to eliminate the baryon load 
problem |385|. 

n 

Beloborodov [22| considers a situation when T n « Tq, the initial Lorentz factor of the 
protons. In this case the decaying neutrons' products will be ahead of the shell of the 
protons. The decaying products will interact with the surrounding matter and will begin to 
slow down. There will be a triple interaction between the two shells and the surrounding 
ambient medium (resembling to some extend the pre-acceleration scenario described earlier) 
. This will take place at radii of a few times Rdecay and at an observed time of a few 
x Rdecay /2cT 2 « afewseconds/(r n /300) , i.e. extremely early. This will produce brightening 
when the fronts pass Rdecay 

The neutrons could also influence the behavior of the relativistic flow during the prompt 
(internal shocks) phase. Specifically inelastic collisions between differentially streaming pro- 
tons and neutrons can produce pions and eventually of 10 GeV as well as v e of 5 GeV 
I2I \l2>\ . These neutrino fluxes could produce ~ 7 events/year in km3 neutrino detectors. 
GeV photons will also be produced but it is unlikely that they could be detected. 
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VIII. ADDITIONAL EMISSION FROM GRBS 



A. TeV 7-rays 
Hurley 



172| reported on detection of 18 GeV photons from GRB 940217. Milagrito - A 
TeV detector - discovered a possible TeV signal coincident with GRB 970417 kjl. Gonzalez 

n 

et al. [130] discovered a high energy tail that extended up to 200 MeV from GRB 941017. 
A natural source for high energy 7-rays is the SSC (Synchro tron s elf C ompton) component 



produced by IC from the burst itself or from the afterglow |235l |260J. The SSC photons 
energy should be 7^ higher than the synchrotron photons. Typical random Lorentz factors of 
electrons, j e , within internal shocks are of order a thousand (in the fluid's rest frame). This 
implies that if the observed 7-rays emission is produced by synchrotron in internal shocks 
then the IC emission would produce a second peak around a few hundred GeV. This would 
be the analogue of the high energy component observed in Blazars. Note that emission abov e 
~ 10 — lOOGeV might be self absorbed by pair production with the source 29J, 3od ] . 

The SSC component would be even higher from the early afterglow. The synchrotron 
emission from the forward shock is expected to be around 10 keV (if the observed early 
afterglow is indeed produced by the external shocks). With a Lorentz factor of a typical 
electron around 10 5 the expected SSC component should be arou nd lOOTeV. Finally the 
reverse shock emission is expected to produce 100 eV photons [37^. With typical electrons 
Lorentz factor of a few thousand this should correspond to SSC photons with typical en- 
ergy of 100 MeV. Depending on the relevant Y parameter the fluxes of these high energy 
components should be comparable or even larger than the prompt GRB 7-rays fluxes. This 
emission should be simultaneous with the GRB emission. It is also possible that the forward 
shock electrons will IC reverse shocks pho t ons. It is likely that this is the cause of the high 



energy emission seen in GRB 941017 2951 1312 1. 



Other mechanisms can produce high energy emission as well. Vietri [423] suggested that 
as GRB can accelerate protons up to 10 20 eV (see §VIIICI below). These protons can emit 
0.01 of the GRB energy as high energy 7-rays with energies up to 300 GeV. Bottcher and 
Dermer J41J considered the synchrotron spectrum resulting from high energy protons and 
leptons produced in cascade initiated by photo-pion production. They predict a significant 
flux of lOMev-lOOGeV photons. 



100 



While the high energy photons flux could be significant these photons might not be 
detectable on earth. The high energy photon flux above 1 TeV would be attenuated sig- 
nific antly due to pair production of such high energy photons with the intergalactic NIR 
flux [l32j. Dai and Lu 66] suggest that secondary emission produced via these interaction 
(upscattering of the CMB by the produced pairs) would still point towards the initial di- 
rection and hence might be detectable as a delayed GeV emission. However, even a tiny 
intergalactic mag netic field (> 1CT 22 G would be sufficient to deflect the electrons and dilute 
these signal [150]. 



B. Neutrinos 



Neutrinos can be produced in several regions within GRB sources. First some models, 
like the Collapsar model or the neutron star merger model predict ample (~ 10 53 ergs) 
production of low (MeV) neutrinos. However, no existing or planned detector could see 
those from cosmological distances. Furthermore, this signal will be swamped in rate by the 
much more frequent SN neutrino signals which would typically appear closer. 

However, GRBs could be detectable sources of high energy neutrinos, with energies rang- 
ing from 10 14 eV to 10 17 eV. These neutrinos are produced by internal or external shocks of 
the GRB process itself and hence are independent of the nature of the progenitor. 

To understand the process of neutrino emission recall that neutrinos are "best" produced 
in nature following pions production in proton-photon or proton-proton collisions. The 
proton-photon process requires that the photon's energy is around the A resonance in the 
proton's energy frame: namely at ~ 200MeV. The resulting pion decays emitting neutrinos 
with a typical energy of ~ 50 MeV in the proton's rest frame. If the proton is moving 
relativistically, with a Lorentz factor 7 p within the laboratory frame the required photon 
energy in the lab frame is smaller by a factor of 7 P and the resulting neutrino energy is 
larger by a factor of 7 p . Depending on the surrounding environment very energetic pions 
may lose some of their energy before decaying producing a "cooling break" in the neutrino 
spectrum. In this case the resulting neutrinos' energy will be lower than this naive upper 
limit. 

Within GRBs protons are accelerated up to 10 20 eV 



422 



429|. The relevant Lorentz 



factors of these protons range from T up to 10 (at the very high energy tail of the protons 
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distribution). Thus we expect neutrinos up to 10 19 eV provided that there is a sufficient 
flux of photons at the relevant energies so that the pions can be produced and there are no 
energy loses to the pions. 



Paczynski and Xu 2891 ] and Waxman and Bahcall |433j ] calculated the flux of VHE 
neutrinos from internal shocks. They found that a significant flux of ~ 10 14 eV neutrinos 
can be produced by in terac tion of the accelerated internal shocks protons with the GRB 
photons. Guetta et al. jl55^ estimate that on average each GRB produces a flux of ~ 10~ 9 
GeV/cm 2 sec sr corresponding to 0.01 events in a km 3 detector. Calculations of specific 
fluxes from indiv idual bursts (that factor in t he ob served 7-rays spectrum) were performed 



by Guetta et al. 



1511 ] . Waxman and Bahcall 434j suggest that protons accelerated at the 



reverse shock (that arises at the beginning of the afterglow) would interact with the optical 
- uv flux of the afterglow and produce 10 18 eV neutrinos. 



24l|, Razzaque et al. \33l\ suggested 



Within the Collapsar model Meszaros and Waxman 

that as the jet punches the stellar shell it can p rodu ce a flux of TeV neutrinos. Within 

I 

the Supranova model the internal shock protons 1491 ] or external shocks protons [77] can 
also interact with external, pulsar wind bubble, photons producing 10 16 eV neutrinos with a 



comparable detection rate to the one obtained form interaction of the internal shock protons 
with 7-rays photons. If the external magnetic field is sufficiently large (as in the pulsar wind 
bubble) external shocks can also accelerate protons to high energy 4251 ] . In this cas e th e 



protons can interact with afterglow photons and can produce neutrinos up to 10 eV 



217|. 



C. Cosmic Rays and Ultra High Energy Cosmic Rays 



Already in 1990 Shemi and Piran 



385] noticed that a fireball may produce cosmic rays. 



However the flux of "low" energy (up to 10 14 eV) that they considered was smaller by 
several orders of magnitude than the observed flux of cosmic rays that are accelerated in 
SNRs. Henc e thi s component isn't important. 



4291 ] and independently Vietri [4221 ] noticed that protons can be accelerated up 



Waxman 

to 10 20 eV within the relativistic shocks that take place in GRBs. Namely internal shocks 
or the reverse shock in GRBs are among the few locations in the Universe where the shock 
acceleration condition (Eq. IT41 needed to accelerate protons up to 10 20 eV, the Hillas criterion 
can be satisfied. Moreover to within an order of magnitude the flux of 7-rays reaching earth 
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from GRBs is comparable to the observed flux of UHECRs (Ultra High Energy Cosmic 
Rays) \i2 9^ . Thus, if GRBs produce a comparable energy in 7-rays and in UHECRs they 
could be th e sou rce of the highest energy Cosm ic rays. 



Greisen 146] and Zatsepin and Kuzmin 445| (GZK) pointed out that the highest energy 
CR (above 10 19,5 eV) are attenuated as they propagate via the Cosmic Microwave back- 
ground (CMBR). This happens because at this high energies the protons can interact with 
the CMBR photons and produce pions. The typical mean free path of a ultra high energy 
proton in the CMBR decreases rapidly with energy and for a 10 20 eV proton it is only several 
tens Mpc. Thus, the observed UHECRs at energies above the GZK energy (~ 10 19 5 eV must 
arrive from relatively nearby (on cosmological scale) sources. However, there are no known 
steady state sources within this distance (see however Farrar and Piran GRBs as a 

transient phenomenon could be a "hidden" source of UHECRs. There won't be a direct 
association between GRBs and arrival of UHECRs as the later are deflected by the inter- 
galactic magnetic field. This leads to an angular deflection as well as a long time delay. If 
GRBs are sources of UHECRs then we expect a break in the UHECR spectrum at the GZK 
energy - below the GZK energy we will detect UHECRs from the whole universe. Above 
the GZK energy we will detect only "local" UHECRs from within the nearest several dozen 
Mpc. Bahcall and Waxman Q| suggested that recent observations imply that such a break 
has been seen. However, the observational situation is not clear as yet and a final resolution 
would most likely require the Auger UHECR detector. 



D. Gravitational Radiation 



Like GRBs, typical sources of gravitational radiation involve the formation of compact 
objects. Hence it is reasonable to expect that gravitational waves will accompany GRBs. 
This association is indirect: the gravitational waves are not directly related to the GRB. 
Additinionally, GRBs have their own, albeit weak, gravitational radiation pulse which arises 
during the acceleration of the jets to relativistic velocities. Unfortunately this signal is weak 
and moreover it is perpendicular to the GRB signal. 

To estimate the rates of observed gravitational radiation events associated with GRB we 
use the rate of long GRBs. The nearest (long) GRB detected within a year would be at 
lGpc. As GRBs are beamed the nearest (long) event would be at would be much nearer, 
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38l| estimates that the rate of short 



at 135#o,iMpc. However, this burst would be directed away from us. Still a GRB that is 
beamed away from us is expected to produce an "orp han" afterglow. 

The rate of short bursts is less certain. Schmidt 
GRBs is smaller by a factor of two than the rate of long ones. In this case the distances 
mentioned above should be revised up by a factor of 1.25. However, if the rate of short 
GRBs is larger by a factor 10 than the rate of long ones then the corresponding distances 
should be revised downwards by a factor of lO -1 / 3 This would put one event per year at 
~ 80#q j^Mpc, but once again this burst won't be pointing towards us. The nearest event 
with a burst in our direction would be at ~ 450Mpc. 



1. Gravitational Radiation from Neutron Star Mergers 

Binary neutron star mergers are the "canonical" sources of gravitational radiation emis- 
sion. LIGO and VIRGO both aim in detecting these sources. Specifically the goal of these 
detectors is to detect the characteristic "chirping" signals arising from the in-spiraling phase 
of these events. The possibility of detection of such signals has been extensively discussed 
(see e.g. jfi^)- Such events could be detected up to a distance of several tens of Mpc with 
LIGO I and up to ~ lOOMpc with LIGO II. 

Comparing with GRB rates we find that if, as some expect, neutron star mergers are 
associated with short GRBs and if the rate of short GRBs is indeed large, then we have 
one event per year within the sensitivity of LIGO II and marginally detectable by LIGO I. 
However, this burst will be pointing away from us. 

The detection of the chirping merger signal is based o n fitt ing the gravitational radiation 
signal to pre-calculated templets. Kochanek and Piran |l 9€j| suggested that the detection 
of a merger gravitational radiation signal would require a lower S/N ratio if this signal 
coincides with a GRB. This would increase somewhat the effective sensitivity of LIGO and 
VIRGO to such events. Finn et al. [97] suggest using the association of GRBs with sources of 
gravitational waves in a statistical manner and propose to search for enhanced gravitational 
radiation activity towards the direction of a GRB during the short periods when GRBs are 
detected. Given the huge distances of observed GRBs it is not clear if any of these techniques 
will be useful. 
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2. Gravitational Radiation from Collapsars 



The Collapsar model |242tl287l l443] is based on the collapse of the core of a massive star to 
a black hole surrounded by a thick massive accretion disk. As far as gravitational radiation is 
concerned this system is very similar to a regu lar supernova. Rotating gravitational collapse 



has been analyzed by Stark and Piran 



397]. They find that the gravitational radiation 



emission emitted in a rotating collapse to a black hole is dominated by the black hole's 
lowest normal modes, with a typical frequency of 20c 3 /GM. The total energy emitted is: 

AE GW = eMc 2 = min(1.4 • 10" 3 a 4 , e max )Mc 2 , (119) 

where a is the dimensionaless specific angular momentum and e max is a maximal efficiency 
which is of the order afew x 10~ 4 . The expected amplitude of the gravitational radiation 
signal, h, would be of the order of \feGM / 'c 2 d where d is the distance to the source. Even 
LIGO II won't be sensitive enough to detect such a signal from a distance of lGpc or even 
from 100 Mpc. 

3. Gravitational Radiation from Supranova 

According to the Supranova model a GRB arises after a neurton star collapse to a black 
hole. This collapse takes PLACE several weeks or months after the Supernova that formed 
he neutron star fsee IPX E*| . The expected gravitational waves signal from a Supranova 



426] includes two components. First the signal from the initial supernova is similar to 
the gravtitational waves from the collapsar model. However, here the first collapse (the 
Supernova) takes place several weeks or months before the GRB. Thus, there won't be any 
correlation between the gravitational waves emitted by the first collapse and the GRB. A 
second component may arise from the second collapse from the supramassive neutron star 
to a black hole. This signal should coincide with the GRB. 

4- Gravitational Radiation from the GRB 

The most efficient generation of gravitational radiation could take place here during the 
acceleration phase, in whic h th e mass is accelerated to a Lorentz factor T. To estimate this 



emission I follow Weinberg 



437j analysis of gravitational radiation emitted from a relativistic 
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collision between two particles. Consider the following simple toy model: two particles at 
rest with a mass M that are accelerated instantly at t = to a Lorentz factor T and energy 
E. Conservation of energy requires that some (actually most) of the rest mass is converted 
to kinetic energy during the acceleration and the rest mass of the accelerated particle is 
m = E/Y = M/Y. The energy emitted per unit frequency per unit solid angle in the 
direction at an angle a relative to f3 is: 

dE GM 2 (3 2 T 2 ((3 2 -cos 2 a) cos 2 a . 



dttdu cvr 2 L (l -/5 2 cos 2 a) 2 r 2 (l - (3 2 cos 2 a) 21 

The result is independent of the frequency, implying that the integral over all frequency 
will diverge. This nonphysical divergence arises from the nonphysical assumption that the 
acceleration is instantaneous. In reality this acceleration takes place over a time St, which 
is of order O.Olsec. This would produce a cutoff uj max ~ 27r/5t above which Eq. 11201 is 
not valid. The angular distribution found in Eq. 11201 is disappointing. The EM emission 
from the ultrarelativistic source is beamed forwards into a small angle 1/r, enhancing the 
emission in the forwards direction by a large factor (r 2 ). The gravitational radiation from 
this relativistic ejecta is spread rather uniformly in almost all 47r steradians. Instead of 
beaming there is "anti-beaming" with no radiation at all emitted within the forward angle 
1/r along the direction of the relativistic motion. 

Integration of the energy flux over different directions yields: 

As expected the total energy emitted is proportional to m 2 Y 2 . Further integration over 
frequencies up to the cutoff 2it/St yields: 

„ 2GM 2 r 2 , , 

E » — . (122) 

CTTOi 

In reality the situation is much more complicated than the one presented here. First, the 
angular width of the emitted blobs is larger than 1/r. The superposition of emission from 
different directions washes out the no emission effect in the forward direction. Additionally 
according to the internal shocks model the acceleration of different blobs go on independently. 
Emission from different blobs should be combined to get the actual emission. Both effects 
reduce the effective emission of gravitational radiation and makes the above estimate an 
upper limit to the actual emission. 
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The gravitational signal is spread in all directions (apart from a narrow beam along the di- 
rection of the relativistic motion of the GRB). It ranges in frequency from to f ma x ~ 100Hz. 
The amplitude of the gravitational radiation signal at the maximal frequency, f ma x ~ 100Hz, 
would be: h « (GMT 2 /c 2 d). For typical values of E = MT = 10 51 ergs, 5t = 0.01 sec and a 
distance of 500 Mpc, h « .5- 10~ 25 , it is far below the sensitivity of planned gravitational ra- 
diation detectors. Even if we a burst is ten times nearer this " direct" gravitational radiation 
signal would still be undetectable . 

Some specific mo dels for GRBs' inner engine p redic t additional amount of energy. For 
example van Putten |-417| . van Putten and Levinson 418 ] suggest a model of a black hole - ac- 
cretion torus in which a large fraction of the emitted energy of the black hole - accretion torus 
system escapes as gravitational radiation. The radiation arises due to instabilities within the 
torus that break down the axial symmetry. They estimate that as much as 10 53 ergs would be 
emitted as gravitational radiation which will have a characteristic signature corresponding to 
the normal mode of the black hole - accretion torus system with typical frequencies around 
few hundred Hz, conveniently within the frequency range of LIGO/VIRGO. If cor rect than 
GRBs are the most powerful burst-sources of gravitational waves in the Universe |417j |. 



IX. MODELS OF INNER ENGINES 



The Fireball model tells us how GRBs operate. However, it does not answer 
the most interesting astrophysical question: what produces them? which as- 
trophysical process generates the energetic ultrarelativistic flows needed for the 
Fireball model? Several observational clues help us answer these questions. 

• Energy: The total energy involved is large ~ 10 51 ergs, a signifi- 
cant fraction of the binding energy of a stellar compact object. The "in- 
ner engine" must be able to generate this energy and accelerate ~ 1O _5 M 
(or the equivalent in terms of Poynting flux) to relativistic velocities. 

• Collimation: Most GRBs are collimated with typical opening angles 
1° < 9 < 20°. The "inner engine" must be able to collimate the relativistic flow. 

• Long and Short Bursts: The bursts are divided to two groups according to their 
overall duration. Long bursts with T > 2sec and short ones with T < 2sec. As the duration 
is determined by the inner engine this may imply that there are two different inner engines. 
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• Rates: GRBs take place once per 3 • 10 5 yr per galaxy. 
GRBs are very rare at about 1/3000 the rate of supernovae. 

• Time Scales: The variability time scale, St, is at times as short as 1ms. The 
overall duration (of long GRBs), T, is of the order of 50sec. According to the internal 
shocks model these time scales are determined by the activity of the "inner engine". St ~ 1 
ms suggests a compact object. T ~ 50sec is much longer than the d ynamical time scale, 
suggesting a prolonged activity. 12 . This requires two (or possibly three \270\. \33t\l ) different 
time scales operating within the "inner engine" . This rules out any "explosive" model that 
release the energy in a single explosion. 

These clues, most specifically the last one suggest that GRBs arise due to accretion of a 
massive (~ O.lm ) disk onto a compact object, most likely a newborn black hole. A compact 
object is required because of the short time scales. Accretion is needed to produce the two 
different time scales, and in particular the prolonged activity. A massive (~ O.lm ) disk is 
required because of the energetics. Such a massive disk can form only simultaneously with 
the formation of the compact object. This leads to the conclusions that GRBs accompany 
the formation of black holes. This model is supported by the observations of relativistic (but 
not as relativistic as in GRBs) jets in AGNs, which are powered by accretion onto black 
holes. 

An important alternative to accretion is Usov's model 



412 



4131 ] in which the relativistic 



flow is mostly Poynting flux and it is driven by the magnetic and rotational energies of a 
newborn rapidly rotating neutron star. 



A. Black hole accretion 

Several scenarios could lead to a black 
include mergers (NS-NS binaries 0, 



role - massive accre tion disk system. This could 
, NS-BH binaries [3 WD-BH binaries Irbl 
BH- He-st ar binaries 112 1) and models based on "failed supernovae" or "Collapsars" [242, 
287l I443I . Narayan et al. 277 1 have recently shown that accretion theory suggests that 
from all the above scenarios only Collapsars could produce long bursts and only NS-NS (or 
NS-BH) mergers could produce short bursts. The basic idea is that the duration of the 



12 



The ratio 5t/T < 1 for short bursts as well |26 
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accretion depends on the size of the disk. So short burst must be produced by small disks 
and those are natureally produced in mergers. On the other hand long burst require large 
disks. However, those are inefficient. One can overcome this if we have a small disk that 
is fed continuously. In this case the efficiency can be large and the duration long. This 
happens naturally within the collapsar model. 



B. The Pulsar Model 



Several "inner engine" models involve pulsar like activity of the inner engine which is 
directly connected to a Poynting flux decimated relativistic flow (in a contrast to a baryonic 
flux dominated flow). Energy considerations require an extremely large magnetic fields of 
the order of 1 15 G within such sources. 

Usov l^l^ suggested that GRB arise during the formation of rapidly rotating highly 
magnetized neutron stars. Such objects could form by the gravitational collapse of accreting 
white dwarfs with anomalously high magnetic fields in binaries, as in magnetic cataclysmic 
binaries. The rapidly rotating and strongly magnetized neutron stars would lose their ro- 
tational kinetic energy on a timescale of seconds or less In a pulsar like mechanism. The 
energy available in this case is the rotational and magnetic energies of the neutron star that 
are of the order of a few xl0 51 ergs for a neutron star rotating near breakup. The rotation 
of the magnetic field creates a strong electric field and an electron-positron plasma which is 
initially optically thick and in quasi-thermodynamic equilibrium. Additionally a very strong 
magnetic field forms. The pulsar produces a relativistic Poynting flux dominated flow. 

W hile a Poynting flux domi nated flow may be dissipated in a regular internal shocks. 



Usov 413j | and Thompson [406J discuss a scheme in which the energy is dissipated from the 
magnetic field to the plasma and then via plasma instability to the observed 7-rays outside 
the 7-rays photosphere, which is at around 10 13 cm. At this distance the MHD approximation 
of the pulsar wind breaks down and intense electromagnetic waves are generated. The 
particles are accelerated by these electromagnetic waves to Lorentz factors of 10 6 and produce 



387, 388] and Drenkhahn and Spruit 



the non ther mal s pectrum. Smolsky and Usov 

Spruit et al. jsQ.^j ] discuss various aspects of the conversion of the Poynting flux energy to 
7-rays but these issues are more related to the nature of the emitting regions and only 
indirectly to the nature of the inner engine. 
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Usov's model is based on rotating highly magnetized neutron star and from this point of 
view it indeed resembles to a large ex tend a regular pulsar. Other authors consider pulsar 
like activities in other contexts. Katz jl85^ . for example, considers a black hole - thick disk 
model in which the electromagnetic proces s tur n rotational energy to particle energy in a 



pulsar like mechanism. Meszaros and Rees 



discuss related idea on the formation of a 



Poynting flux dominated flow within a black hole accretion disk model. 



C. Rotating black holes and the Blandford Znajek mechanism 



It is possible and even likely that the process of energy extraction involves the Blandford- 
Znajek mechanism 32] in which the black hole - torus system is engulfed in a magnetic 
field and the rotational energy of the black hole is extracted via this magnetic field. The 
exploration of the Blandford-Znajek mechanism involves relativistic MHD consideration 

which are beyond the scope of this review. I refer the reader to several recent extensive 

reviews on this subject (see e.g. Lee et al. |214j ). 



D. The Collapsar Model 

The evidence for the association of (long) GRBs with supernovae (s ee Bl oom et al. 
and flll C 4|) provides a strong support for the Collapsar model. Woosley 4431 ] proposed that 
GRB arise from the collapse of a sing le Wolf-Rayet star endowed with fast rotation ('failed' 
Type lb supernova). Paczynski |287| pointed out that there is tentative evidence that the 
GRBs 970228, 970508, and 970828 were close to star-forming regions and that this sug gests 
that GRBs are linked to catacly smic deat hs of massive stars. MacFadyen and Woosley 242] 



begun a series of calculations 



4. 



243, 



447j of a relativistic jet propagation through the stellar 



envelope of the collapsing star which is the most important ingredient unique to this model 
(other features like the accretion process onto the black hole, the corresponding particle 
acceleration and to some extend the collimation process are common to other models). The 
collimation of a jet by t he st ellar mantle was shown to occur analytically by Meszaros and 



Rees 



240j. Zhang et al. numerically confirmed and extended the basic features of this 



collimation process. 

According to the Collapsar model the massive iron core of a rapidly rotating massive 
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star, of mass M > 30M©, collapses to a black hole (either directly or during the accretion 
phase that follows the core collapse). An accretion disk form around this black hole and a 
funnel forms along the rotation axis, where the stellar material has relatively little rotational 
support. The mass of the accretion disk is around 0.1 M . Accretion of this disk onto the 
black hole takes place sev eral dozen seconds and powers the GRB. Energy can be extracted 



via neutrino annihilation 



242| or via the Bladford-Znajek mechanism. The energy deposited 



in the surrounding matter will preferably leak out along the rotation axis producing jets with 
opening angles of < 10°. If the jets are powerful enough they would penetrate the stellar 
envelope and pr oduc e the GRB. 



Zhang et al. 



44 7| find that relativistic jets are collimated by their passage through the 



stellar mantle. Starting with an initial half-angle of up to 20°, the jet emerges with half- 
angles that, though variable with time, are around 5°. The jet becomes very hot in this 
phase and it has only a moderate Lorentz factor, modulated by mixing, and a very large 
internal energy (more than 80% of the total energy). As the jet escapes, conversion of the 
remaining internal energy into kinetic energy gives terminal Lorentz factors along the axis 
of ~ 150 (depending, of course, on the initial conditions considered). Because of the large 
ratio of internal to kinetic energy in both the jet and its cocoon, the opening angle of the 
final jet is significantly greater than at breakout. A small amount of material emerges at 
large angles, but with a Lorentz factor still sufficiently large to make a weak GRB. When 
;he jet breaks out from the star it may produce a thermal precursor (seen in several GRBs) 
H Q. .nstabjlities in the aectetion ptoeess, or in the passage of the jet through 
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the stellar envelope |3j, l447| can produce the required variability in the Lorentz factor that 
is needed to produce internal shocks. 

The processes of core collapse, accretion along the polar column (which is essential in 
order to creat e the funnel) and the jet propagation through the stellar envelope take together 
~ lOsec I242I . The duration of the accretion onto the black hole is expected to take several 
dozen seconds. These argu ment s imply that Collapsars are expected to produce long GRBs 
(see however, Zhang et al. for a suggestion that the breakout of a relativistic jet and 
its collision with the stellar wind will produce a brief transient with properties similar to 
the class of "short-hard" GRBs.). 



Ill 



E. The Supranova Model 



Vietri and Stella 



4261 ] suggested that GRBs take place when a "supermassive" (or supra- 



massive as Vietri and Stella 



426( | call it) neutron star (namely a neutron star that is above 



the maximal cold nonrotating neutron star mass) collapses to a black hole. The collapse 
can take place because the neutron star losses angular momentum via a pulsar wind and 
it looses the extra support of the centrifugal force. Alternatively the supramassive neutron 
star can simply cool and become unstable if rotation alone is not enough to support it. The 
neutron star could also beco me ov er massive and collapse if it accretes slowly matter from 
a surrounding accretion disk 4271 ]. In this latter case the time delay from the SN could be 
very large and the SNR will not play any role in the GRB or its afterglow. 

The Supranova model is a two step event. First, there is a supernova, which may be 
more energetic than an average one, in which the supermassive neutron star forms. Then 
a few weeks or months later this neutron star collapses producing the GRB. While both 
the Supranova and the Collapsar (or hypernova) events are associated with Supernovae 
or Supernovae like event the details of the model are very different. First, while in the 
Collapsar model one expect a supernova bump on the afterglow light curve, such a bump 
is not expected in the Supranova model unless the time delay is a few days. On the other 
hand while it is not clear in the Collapsar model how does the Fe needed for the Fe X- 
ray lines reach the implied large distances form the center, it is obvious in this model, as the 
supernova shell was ejected to space several month before t he GRB. As mentioned earlier 
(see 311 C4I) the association of GRB 030329 with SN 2003dh [167] |395| is incompatible with 
the Supranova model. Proponents of this model, argue however, that there might be a 
distribution of delay times between the first and second collapses. 

The models are also very different in their physical content. First in the Supranova 
model the GRB jet does not have to punch a whole through the stellar envelop e. Instead 
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the ejecta propagates in almost free space polluted possibly by a pulsar wind 
In both models, like in many other models, the GRB is powered by accretion of a massive 
accretion disk surrounding the newborn black hole. This accretion disk forms, from the 
debris of the collapsing neutron star at the same time that the black hole is formed. Again, 
the t ime scale of the burst is determined by the accretion time of this disk. Narayan et al. 



277] (see also §IX A J) point however that long lived (50 sec) accretion disks must be large 
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and hence extremely i neffic ient. This may pose a prob lem for this model; 

Konigl and Granot jl83| . Guetta and Granot 149 1 and Inoue et al. |l78| considered the 
effects a strong pulsar wind (that may exist after the SN and before the second collapse) 
on this scenario. The pulsar wind can have several effects. First it would produce a denser 
highly magnetized medium into which the GRB jet propagates. The strong magnetic field 
will be amplified by the afterglow shock. This resolves the problem of the source of the 
strong magnetic field needed for the synchrotron afterglow model. This can also explain 
the high energy emission detected by EGRET in GRB 940217 (Hurley [172J and 311 A II) 
by Inverse Compton scattering on the pulsar wind bubble photons. On the other hand the 
density of this wind matter (~ 10 3 cm~ 3 ) might be too high for the spherical model. Note 
however, that this depends on the time delay as t~ 3 . However, the pulsar wind won't be 
spherical and one would expect that it will form an elongated supernova shell cavity within 
which the pulsar wind is bounded. If, as expected, the pulsar jet coincides with the GRB 
jet then the relativistic ejecta will move along the elongated direction of this shell. 



F. Merging neutron stars 



Neutron star binary mergers 



861 |276] or neutron star - black hole binary mergers [286] 
(hereafter called mergers) also produce a black hole - accretion disk system and are candi- 
dates for the inner engines of GRBs, specifically of short GRBs. These mergers take place 
because of the decay of the binary orbits due to gravitational radi ation emission as was 



405]. 



278, 



300 



414 1 . 



beautifully demonstrated in the famous binary pulsar PSR 1913+16 

These mergers take place at a rate of ~ 10~ 6 events per year per galaxy 
This is the rate of merger of binaries of the type of PSR 1913+16 whose life time is or 
order of several 10 8 years. Various population synthesis calcula tions suggest that there is 



also another population of short lived binaries |18L |2i 



410 



411]. These binaries form with 
very close orbits and hence with short lifetimes of the order of 10 5 yrs. Even thought the 
overall rate of such mergers could be comparable to those of the PSR 1913+16 type one 
cannot expect to catch a binary in our galaxy in such a stage. Similarly unlike the l ong lived 



276| this 



mergers that may be kicked from their host galaxy within their lon g li fe time (44 , 
short lived population remains within the galaxy when they merge |18j ]. 

Earlier simulations of mergers focused on the gravitational radiation from this system. 
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Davies et al. 



721 ] begun a ser ies o 



on GRB related aspects [lH, 1349 



nu meric al simulation of neutron star merger that focused 



350, 



351|- Using a SPH scheme they followed NS mergers 
under different assumptions (Newtonian with ad hoc addition of gravitational radiation 
back reaction or Post Newtonian), with different equations of state (adiabatic or realistic) 
and with different initial spin axis and mass rations and different estimates of the effects 



of ne utrino cooling. A 



para, 



180], Ruffert and Janka 



354, 



lei set o f simulations was carried out by Janka and Ruffert 



355 



357] who used particle in cell methods. 



3561 ] , Ruffert et al. 

Both kinds of simulations yield comparable results. The merger results in a black hole - 
accretion disk system. The mass of the accretion disk is of order 0.1M Q and it depends, of 
course somewhat on the orientation of the spins and the relative masses of the two neutron 
stars. 

A merger releases ~ 5 x 10 53 ergs but most of this energy is in the form of low energy 
neutrinos and gravitational waves. Still there is enough energy available to power a GRB 
but is not clear how the GRB is produced. A central question is, of course, how does a 



merger generate the relativistic wind required to power a GRB. Eichler et al. [8J| suggested 
that about one thousandth of these neutrinos annihilate and produce pairs that in turn 
produce gamma-rays via vv —>■ e + e~ — > 77. This idea was criticized on several grounds by 
different aut hors the main problem is that it does not produce enough energy. For example 



Jaroszynksi 



1811 ] pointed out that a large fraction of the neutrinos will be swallowed by 



the black hole that forms. An alternative source of energy within the merger model is the 
accretion power of a disk that forms around the black hole. This brings us back to the 
canonical black hole - accretion disk scenario. 



X. OPEN QUESTIONS AND FUTURE PROSPECTS 



I believe that overall we have a basic understanding of the GRB phenomenon. As usual 
some aspects are understood better than others. 

There is a very good understanding of the afterglow. Here there are numerous observa- 
tions over a wide range of wavelengths with which the theory can be confronted. The overall 
picture, of a slowing down relativistic flo w an d of synchrotron emis sion fit the data to a large 
extent (see e.g. Panaitescu and Kumar 2911 ] . Wijers and Galama 440] and many other fits 
of the observations to the model). We have already learned that the "cow is not spherical", 
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namely that the relativistic flow is collimated. New observations, like those of GRB 021004 
and GRb 030329, poses at times new puzzles and suggest that the basic simple picture has 
to be refined. It seems however, that the answers are within the scope of the current model, 
such as: refreshed shocks, patchy shells and variable external densities. All these phenom- 
ena are fairly reasonable in a realistic environment. Within the afterglow, I believe that the 
X-ray lines pose the greatest current puz zle, i n terms of their energy requirements and other 
implications on the source (see Lazzati 209]). Another interesting open question is what 
distinguished between GHOSTs and OTGRBs - an environmental (extinction or "improper 
conditions within the circum-burst matter) or an intrinsic mechanism? 

The main observational challenges concerning the afterglow are the determination 
whether short GRBs have afterglow. A wealth of information on long GRBs arises from the 
information on hosts, environments and redshifts, that are determined from the afterglow 
observations. All these are missing for short GRBs. If short GRBs don't have afterglows, 
then an immediate theoretical question is why? Is it possible that they are produced in 
a very different environment than long ones (such as outside their original galaxies) in a 
region with no circum-burst matter suitable for producing the afterglow? At the moment 
the observational situation is not clear. The coming Swift satellite may resolve this mystery. 

Another important observational question involves the search for orphan afterglows (ei- 
ther in radio or in optical). Their detection will establish the collimated jets picture. But 
even firm upper limits will set independent limits on the rates of GRBs. However, as 
mentioned in §VII Kl this is a very challenging observational task. This has important im- 
)lication f or th e nature of the jets - are GRB jets standard with a fixed angular structure 
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347, 446]? This question is related both to the overall energetics and to the rate of 



GRBs. 

Another interesting challenge will be the resolution of the afterglow image (see Granot 
et al. |l4l| ). This may be possible in radio for a nearby burst and the afterglow of GRB 
030329 provides an excellent candidate for that. Polarization measures could pave the way 
for understanding of the collimation geometry and for a confirmation of the synchrotron 
emission process. 

As we move backwards with time towards the burst we encounter the very early afterglow 
and the optical flash that coincides with the burst itself. Here a great progress was made 
with recent observations triggered by HETE II (e.g. the almost complete light curve of 
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GRB 021004 



10l(). SWIFT may contribute a lot to this issue. These observations could 



shed a light on issues like the role of pre-acceleration and neutrons that are unclear as 
yet. Here, I stress the importance of early and continuous radio observations, which could 
determine whet her t here are refreshed shocks during the early afterglow, that have a clear 
radio signature |207 |. 

The understanding of the 7-rays emitting regions is less clear. Here, within the internal 
shocks model there is a reasonable understanding of the temporal structure (in terms of 
activity of the inner engine). However, it is not clear how is the observed spectrum produces 
and it seems t hat the simple s ynchr otron spectrum has to be modified (see e.g. Lloyd 



and Petrosian 



222], Medvedev |253| for ideas on such modifications). Another possibly 



related puzzle is the origin of the n arrow E p distributions (see however, e.g. Daigne and 

Guetta et al. 0). Another set of open questions is what is the 



Mochkovitch 
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origin of the in trinsic correlation between luminosity (which in fact reflects the collimation 



angle |105ll291| ) discovere d by F enimore and Ramirez- Ruiz [2a] or the lag-luminosity relation 



discovered by Norris et al. [280] . Similarly or ev en more puzzling are the implied correlations 
)etw een redshift and intrinsic luminosity 223j] and between redshift and intrinsic hardness 
380] (note that thi s lat er correlation is essential in view of the narrow E p distribution of 

GRBs). Here pairs 1241 ] and IC can play an important role. Theoretical open basic physical 



tgrg low) deal with the pro cesse s 
, Nikto and Medvedev 



questions that arise here (as well as in the theory of the a 
of the behavior of collisionless shocks (see e.g. Medvedev 
particle acceleration (see §VB|) and the generation of strong magnetic field (see |25fi| ). Issues 
like relativistic turbulence and relativis tic p lasma instabilities might play an important role 



226|). 



here (see e.g. Lyutikov and Blandford 

^From an observational point, it will be a challenge to beat the statistical power of the 
BATSE data in terms of number of bursts. Somewhat surprisingly, the questions what is 
the luminosity function of GRBs and what is the rate of GRBs as a function of redshift and 
to what extend GRBs follow the star formation rate are still open. Detectors with better 
spectral resolutions could shed some additional light on the spectrum. Another hope for new 
data, or at least for upper limits, arises from observational windows in higher 7-rays b ands . 
On the low energy side it seems that there is a continuum between XRFs and GRBs [l^L Il89| . 
This result still has to be fully understood in the context of the narrow E p distribution. 
Looking far into the future one can hope to observe neutrinos or gravitational radiation 
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correlated to GRBs. UHE neutrinos (Fluxes of MeV neutrinos would be too weak to be 
detected from cosmological distances) could confirm that protons are accelerated to UHE 
energies within GRBs. In turn this would proof (or disprove) the possible role of GRBs as 
sources of UHECRs. Gravitational radiation could give a direct clue on the activity of the 
inner engine (see £IVlllD4l and identify, for example, merger events. 

There is a lot of observational evidence associating long GRBs with core collapse SNes. 
This gives a clear clue on what is the inner engine of long GRBs. There is no direct or 
indirect evidence on the progenitors of short GRBs. Even with this clue the situation is far 
from clear when we turn to the inner engine. Here most models assume some variant of a 
black-hole - torus system with various energy extraction mechanisms ranging from neutrino 
annihilation (which is less likely) to variants on the theme of electromagnetic extraction 
(magnetic turbulence within the accretion disk; the Blandford-Znajek mechanism which 
involves a disk-black hole- magnetic field interaction; pulsar like activity). Here there are 
open questions all around: What is the content of the ultrarelativistic flow - baryonic or 
Poynting flux? How is the flow accelerated and collimated? What determines the variability 
of the flow (required for internal shocks) and the different time scales? This part of the 
model seems to be in a rather poor shape - but this is understandable as we don't have any 
direct observations of this inner engine. One hope that arises is that there seem to be an 
emerging similarity between GRBs, galactic micro quasars and AGNs. All these systems 
accelerate collimated flows to relativistic velocities and they all seem to involve accretion 
onto black holes. Hopefully, this similarity could lead to a common resolution of how inner 
engines operate in all those systems. 
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Figures 




Break Ene-gy (keV) 



FIG. 1 : Dis tribution of the observed values of the break energy (a — (3)Eq in a sample of bright 
bursts 



. The solid line represents the whole sample while the dashed line represent a subset of the data. 



GRB920627 




Time since trigger (sec) 



FIG. 2: The light curve of GRB 920627. The total duration of the burst is 52sec, while typical 
pulses are 0.8sec wide. Two quiescent periods lasting ~10 seconds are marked by horizontal solid 
bold lines. 
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FIG. 3 : Th e pulse width distribution (right) and the distribution of intervals between pulses (left) 
(from gal). 
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Fig. 1 

FIG. 4: The hardness-duration correlation for BATSE bursts. HR is the ratio of fiuence between 
BATSE channels 3 and 2. Squares represent short burst and circle long bursts. The solid line is 
a regression line for the whole sampl e. T he two dotted lines are the regresions lines for the short 
and long samples respectively (from [s2f]].) 
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FIG. 5: The distribution of X-ray fluxes (2- 10 ke V) at t=ll hours after the GRB in 21 afterglows 
observed by BeppoSAX (from Piran et al. The sample includes all the fast observations 

performed by BeppoSAX on GRB from January 1997 to October 1999. No X-ray afterglow was 
detected in GB990217 to the limiting instrumental sensitivity of 10 -13 erg cm -2 s _1 , 6 hours after 
the burst. In the case of GB970111 a candidate was detected, but evidence of fading behaviour is 
marginal, so both cases are considered as upper limits (indicated by the arrow). 
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FIG. 6: The fraction of bursts with optical afterglow above three limiting magnitudes as a function 
of time (compared to the total number of bursts with optical afterglow). 




FIG. 7: The optical light curves of GRB990510. A fit for the observed optical light curves is 
obtained with a\ = 0.82 ± 0.02, a 2 = 2.18 ± 0.05 and U = 1.2 ± 0.08 days (from Harrison et al. 



141 



Radio 

FIG. 8: Venn diagram of the distribution of 47 Afterglows observed in different wavelength between 
1997-2001 (from a talk given by D. Frail at the Sackler GRB workshop, Harvard, May 2001 
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FIG. 9: Light curves of the radio afterglow of GRB970508 at 4.86GHz and 1.43GHz, compared 
with the predictions of the adiabatic fireball model (from Frail et al. 104]). 
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FIG. 10: Comparison of the peak magnitudes of GRB-SNe with those of local Type Ib/Ic SNe. 
The solid cumulative histogram to the left is for those GRBs with a believable detection of a SN 
bump; the brightness of all other claimed GRB-SN peaks or reported upper-limits are shown as a 
dotted histogram. The solid histogram to the right is for those local Ibc SNe for which the SN was 
observed before peak; all others were discovered after peak. The faintest GRB-SN non-detection 
(010921) only probes the top ~40th-percentile of local Type Ib/Ic SNe. It is clear that the current 
GRB-SNe population may have only revealed the tip of the iceberg; plausibly, then, SNe could 
accompany all long-duration GRBs, from 33] 



143 



d 

on 

d 
o 
U 
+ 

3 
in 




4000 5000 6000 7000 8000 
Observed Wavelength (A) 



FIG. 11: Evolution of the GRB 030329/SN2004dh spectrum, from April 1.13 UT (2.64 days after 
the burst), to April 8.13 UT (6.94 days after the burst). The early spectra consist of a power 
law continuum (F u oc v~ 9 ) with narrow emission lines originating from HII regions within the 
host galaxy at redshift z=0.168. Sp ectra taken after April 5 show the development of broad peaks 



characteristic of a supernova (from 
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FIG. 12: A histogram of GRB energies (E^) with three equal logarithmic spacings per decade. The 
histogram shows a narrow distribution of GRB energies about the standard energy E-y = 1.33- 10 51 
ergs , with an error of a = 0.07 dex. The observed burst-to-burst rms spread is 0.35 dex (a factor 
of 2.23) about this value. Bands of 1, 2, and 5 a about the standard energy are shown. There are 
five identifiable outliers, which lie more than 5 a from the mean, however, there is currently no 
basis other than discrepant energy 
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FIG. 13: GRB energy release versus redshift. Bands of 1, 2, and 5 a about the mean energy E^ — 
1.33 • 10 51 ergs are shown. Plotted are the trajectories of five GRBs with no known spectroscopic 
redshift (labeled with star symbols). While the energies of GRB 980329 and GRB 000630 could be 
consistent with the standard value at redshifts beyond z ~ 1.5, at no redshift could the energies of 
GRB 980326 and GRB 980519 be consistent citeBloomFrailKulkarni03. 
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FIG. 14: The observed pulse from an instantaneous flash from a spherical relativistic thin shell 
moving relativistically and emitting emitting with a power low 
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1341 ] ). The circle marks the and 



FIG. 15: Polarization from a uniform magentic field (following 
angle where the matter moves at an angle T _1 from the observer. Note that the polarization is 
maximal along the line perpendicular to the uniform field but it vanished in the other direction 
near the T _1 circle. 
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FIG. 16: Polarization from a random magnetic field in the plane of the shock [275(1 • The solid circle 
marks the and angle where the matter moves at an angle T _1 from the observer. Polarization can 
arise if we view a part of a jet (dashed line), or if the emission is dominated by hot spot (dashed 
- dotted regions) or if there is an overall gradient of the emissivity as would arise in the standard 
jet model. 




FIG. 17: The time integrated polarization (solid line) and the efficiency (dashed line) as a function 
of 6 f, s T for four different values of 9j for a random magnetic field. 
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FIG. 18: Different time scale from a relativistic expanding shell in terms of the arrival times (tj) 
of various photons: t ang = tp — tj±, Ir = tc — and t& = ie — 
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FIG. 19: The internal shocks model (from 36J]) Faster, shells collide with slower ones and produce 
the observed 7 rays. The variability time scale is L/c while the total duration of the burst is A/c . 
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FIG. 21: The shot-gun model (from 36J]). The inner engine emits narrow "bullets" that collide 
with the ISM. 




FIG. 22: The Lorentz factor T, the density p and the pressure p in the shocks. There are four 
regions: the ISM (region 1), the shocked ISM (region 2), the shocked shell (region 3) and the un- 
shocked shell (region 4), which are separ ated by the forward shock (FS), the contact discontinuity 
(CD) and the reverse shock (RS). From 194] . 
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FIG. 23: The different possible broad band synchrotron spectra from a relativistic blast wave, 
that accelerates the electrons to a power law distribution of energies. The thin solid line shows 
the asymptotic power law segments, and their points of intersection, where the break frequencies, 
z/&, and the corresponding flux densities, F Ubfixt , are defined. The different PLSs are labelled A 
through H, while the different break frequencies are labelled 1 through 11. The temporal scalings 
of the power law segments and the break frequencies, for an ISM (k = 0) or stellar wind (k = 2) 
environment, are indicated by the arrows. The different spectra are labelled 1 through 5, from top 
to bottom. The relevant spectrum is determined by the ordering of the break frequencies. The 
top two panels (spectra 1 and 2) correspond to slow cooling (y m < v c ). Spectrum 1 applies when 
v S a < "m, while spectrum 2 applies when v m < v sa < v c . The two bottom panels (spectra 4 and 5) 
correspond to fast cooling (v c < f m ). Spectrum 5 applies when v sa < u c , and spectrum 4 applies 
when v c < v sa < v m . Spectrum 3 (middle panel) applies when v sa > v m ,v c , where in this case the 
relative ordering of v r an d v m is unimportant (i.e. spectrum 3 may apply both to slow cooling or 
fast cooling). From 144]. 




FIG. 24: The light curves results from a Gaussian (AR/R = 0.1) over-dense region in the ISM. The 
different thick lines are for a maximal over-densities of 40 (dashed- dot), 5 (dashed) and 2 (solid). 
The thin line is the light curve for a constant ISM density. The inset depicts the ratio of the mass, 
M(R ) over the mass of an ISM (without the over-dense region), Mjsm(R) (from Nakar and Piran 
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FIG. 2 5: O bserved and predicted light curve at 8.6Ghz light curves of GRB990510 (from Harrison 



et al. 



1591 ] ). The different behavior of the optical and radio light curves after the jet break is 



clearly seen. 




FIG. 26: A relativistic jet at the last time step of the simulation [1371 ]. (left) A 3D view of the 
jet. The outer surface represents the shock front while the two inner faces show the proper number 
density (lower face) and proper emissivity (upper face) in a logarithmic color scale, (right) A 2D 
'slice' along the jet axis, showing the velocity field on top of a linear color-map of the lab frame 
density. 
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FIG. 27: The radius R (left frame), Lorentz factor T — 1 (middl e fra me) and opening angle 9 of 
the jet (right frame), as a function of the lab frame time in days |l37j ]. 
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FIG. 28: The light curves of an inhomogeneous jet observed from different angles (From Rossi 
et al. [347! ) . From the top 9 Q = 0.5, 1, 2, 4, 8, 16°. The break time is related only to the observer 
angle: tf, oc 9%. The dashed line is the on-axis light curve of an homogeneous jet with an opening 
angle 29 and an energy per unit solid angle e{0 o ). The blow up is the time range between 4 hours 
and 1 month, where most of the optical observations are performed. Comparing the solid and the 
dashed lines for a fixed 9 Q , it is apparent that one can hardly distinguish the two models by fitting 
the afterglow data. 
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FIG. 29: Light curves for structured jets (initially e oc 9 a and F oc 9 b ), for models 1 
and 2, in the optical (y = 5 x 10 14 Hz), for a jet core angle 9 C = 0.02, viewing angles 
9 ohs = 0.01,0.03,0.05,0.1,0.2,0.3,0.5, p = 2.5, e e = e B = 0.1, n = 1 cm" 3 , r = 10 3 , and e 
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was chosen so that the total energy of the jet would be 10 
in some of the panels, for comparison. From Granot et al. 
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(GK). A power law of t v is added 
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FIG. 30: Left: Shape of the emitting region. The dashed line marks the physical extent of the jet, 
and solid lines give the viewable region I/7. The observed radiation arises from the gray shaded 
region. In each frame, the percentage of polarization is given at the top right and the initial size of 
the jet relative to \/T is given on the left. The frames are scaled so that the size of the jet is unity. 
Right: Observed and theoretical polarization light curves for three possible offsets of the observer 
relative to the jet axis Observational data for GRB 990510 is marked by crosses (x), assuming 
tj e t = 1.2 d ays. The upper limit for GRB 990123 is given by a triangle, assuming tj et = 2.1 days 



FIG. 31: A schematic afterglow light curve. While the burst differ before the jet break (due to 
different opening angles, the light curve coincide after the break when the energy per unit solid 
angle is a constant. 



(from |3J5j). 




Log(t) 



154 



Canonical 
Optimistic 

Canonical + Z peaJ = 2 
Model B - Canonical 
Model B - Optimistic 
Canonical Uncertainly 



FIG. 32: The number of observed orphan afterglows per square degree (left vertical scale) and in 
the entire sky (right vertical scale), in a single exposure, as a function of the limiting magnitude 
for detection. The thick lines are for model A with three different sets of parameters: i) Our 
"canonical" normalization Fq = 0.003 /xJy, z pea k = 1, 0j = 0.1 (solid line). The gray area around 
this line corresponds an uncertainty by a factor of 5 in this normalization, ii) Our most optimistic 
model with a relatively small 0j = 0.05 and a large -Fo = 0.015^/Jy (dashed- dotted line), iii) 
The same as our "canonical" model, except for z pea k = 2 (dotted line). The thin lines are for 
model B, where the solid line is for our "optimistic" parameters, while the dashed line is for our 
"canonical" parameters. Both models are similar for the "optimistic" parameters while mo del B 
predicts slightly more orphan afterglows then model A for the "canonical" parameters. From 27^ 
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FIG. 33: Afterglow from a GRB ejecta decelerating in a wind of a Wolf-Rayet progenitor with 
M = 2 x 1O _5 M0 yr -1 and w = 10 3 km s _1 . The burst has an isotropic energy E = 10 53 erg) 
and a thin ejecta shell with kinetic energy E = 10 53 erg and a Lorentz factor V = 200. Dashed 
curves show the prediction of the standard model that neglects the impact of the radiation front, 
and solid curves show the actual behavior. Two extreme cases are displayed in the figure: 77 = 
(adiabatic blast wave) and rj = 1 (radiative blast wave). Four zones are marked: I - R < R gap (the 
gap is opened), II - R gap < R < Race (the gap is closed and the ejecta sweeps the relativistically 
preaccelerated e ± -loaded ambient medium), III - R acc < R < Ri oa d (e ± -loaded ambient medium 
with T arn t, ient ~ 1), and IV - R > Ri oa d (pair-free ambient medium with T am \, ient rj 1). The radie 
are measured in units of R aC c ~ 10 16 cm. Top panel: dissipation rate. Bottom panel: synchrotron 
peak frequency (assuming e# = 0.1) in units of m e c 2 /h. From Beloborodov 2C] 
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a 


P 


V < Va 


1 


2 


V a < V < U c 


1/6 


1/3 


V c < V < V m 


-1/4 


-1/2 


V m <V 


-(3p-2)/4 


-p/2 = (2a - l)/3 


TABLE I: a and /3 for fast coolin 




a 


P 


v < v a 


1/2 


2 


V a < V < V m 


1/2 


1/3 


v m <v <v c 


-3(p-l)/4 


-(p- l)/2 = 2a/3 


v c <v 


-(3p-2)/4 


-p/2 = (2a - l)/3 



TABLE II: a and P for slow cooling (v a <v m < u c ) into a constant density ISM 



ISM 



wind 



Jet 
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v < v a 


(17p-26)/16(p-l) 


(13p-18)/18(p-l) 


3(p-2)/4(p-l) 


v a <v <v m 


(p+l)/8(p-l) 


5(2-p)/12(p-l) 


(8-5p)/6(p-l) 


v m <v <v c 


-3(p+2)/16 


-(p+8)/8 


-(p+6)/4 


v c <v 


-(3p+10)/16 


-(p+6)/8 


-(p+6)/4 



TABLE III: a for slow cooling [y a < v m < v c ) into a constant density ISM, wind and jet for 
electron distribution with 1 < p < 2. 





a 


P 


V c < V < v m 


-1/4 


-1/2 


V m ,V c < V 


-(3p-2)/4 


-p/2 = (2a - l)/3 


U m < V < v c 


-(3p-l)/4 


-(p-l)/2 = (2a + l)/3 



TABLE IV: a and /? for X-ray and optical frequencies from a blast wave into a wind profile when 

n 

z^ a < v c ,v m ,v |50|]. Note that the order of the table is according to the evolution of the light curve 
at a fixed high observed frequency. 





a 


P 


v c <u <u a < v m 


7/4 


5/2 


V < v c < v a < v m 


2 


2 


V < V a < V m < U c 


1 


2 


V a < V < V m < V c 





1/3 


Va<V m <V <V C 


-(3p-l)/4 


-(p-l)/2 = (2a + l)/3 



TABLE V: a and P for radio frequencies from a blast wave into a wind profile 



50]. Note that the 



order of the table is according to the evolution of the light curve at a fixed low observed frequency. 





a 


P 


V <v a 





2 


v a <v < v m 


-1/3 


1/3 


V m < V < u c 


-P 


-(p-l)/2 = (a + l)/2 


v c <v 


-P 


-p/2 = a/2 



TABLE VI: a and P for slow cooling [y a <v m < u c ) after a jet break. 
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